+// Copyright (C) 2014-2015 EDF-R&D
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
#include "HYDROData_Bathymetry.h"
#include "HYDROData_Document.h"
#include <TDataStd_AsciiString.hxx>
#include <TDataStd_Integer.hxx>
+#include <QColor>
#include <QFile>
#include <QFileInfo>
#include <QPointF>
aCoordsArray->SetValue( i * 3 + 2, aPoint.Z() );
}
- SetToUpdate( true );
+ Changed( Geom_Z );
}
-HYDROData_Bathymetry::AltitudePoints HYDROData_Bathymetry::GetAltitudePoints() const
+HYDROData_Bathymetry::AltitudePoints HYDROData_Bathymetry::GetAltitudePoints(bool IsConvertToGlobal) const
{
AltitudePoints aPoints;
if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), aCoordsArray ) )
return aPoints;
+ Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
for ( int i = aCoordsArray->Lower(), n = aCoordsArray->Upper(); i <= n; )
{
if ( i + 3 > n + 1 )
aPoint.SetY( aCoordsArray->Value( i++ ) );
aPoint.SetZ( aCoordsArray->Value( i++ ) );
+ if( IsConvertToGlobal )
+ aDoc->Transform( aPoint, false );
aPoints.Append( aPoint );
}
void HYDROData_Bathymetry::RemoveAltitudePoints()
{
- TDF_Label aLabel = myLab.FindChild(DataTag_AltitudePoints, false);
- if (!aLabel.IsNull())
- {
- aLabel.ForgetAllAttributes();
- SetToUpdate(true);
- }
+ TDF_Label aLabel = myLab.FindChild( DataTag_AltitudePoints, false );
+ if ( !aLabel.IsNull() )
+ {
+ aLabel.ForgetAllAttributes();
+ Changed( Geom_Z );
+ }
}
void interpolateAltitudeForPoints( const gp_XY& thePoint,
double HYDROData_Bathymetry::GetAltitudeForPoint(const gp_XY& thePoint) const
{
- //DEBTRACE("GetAltitudeForPoint p(" << thePoint.X() << ", " << thePoint.Y() << ")");
+ DEBTRACE("GetAltitudeForPoint p(" << thePoint.X() << ", " << thePoint.Y() << ")");
double anInvalidAltitude = GetInvalidAltitude();
double aResAltitude = anInvalidAltitude;
}
std::map<double, const gp_XYZ*> dist2nodes;
- aQuadtree->NodesAround(thePoint, dist2nodes, 1.0);
+ aQuadtree->NodesAround(thePoint, dist2nodes, aQuadtree->getPrecision());
+ while (dist2nodes.size() == 0)
+ {
+ aQuadtree->setPrecision(aQuadtree->getPrecision() *2);
+ DEBTRACE("adjust precision to: " << aQuadtree->getPrecision());
+ aQuadtree->NodesAround(thePoint, dist2nodes, aQuadtree->getPrecision());
+ }
+ aQuadtree->NodesAround(thePoint, dist2nodes, 5.0);
if (dist2nodes.size())
{
std::map<double, const gp_XYZ*>::const_iterator it = dist2nodes.begin();
aResAltitude = it->second->Z();
- //DEBTRACE(" number of points found: " << dist2nodes.size() << " nearest z: " << aResAltitude);
+ DEBTRACE(" number of points found: " << dist2nodes.size() << " nearest z: " << aResAltitude);
+ }
+ else
+ {
+ DEBTRACE(" number of points found: 0");
}
return aResAltitude;
TDataStd_Integer::Set( myLab.FindChild( DataTag_AltitudesInverted ), (Standard_Integer)theIsInverted );
- SetToUpdate( true );
+ Changed( Geom_Z );
if ( !theIsUpdate )
return;
// Try to import the file
if ( aFileSuf == "xyz" )
aRes = importFromXYZFile( aFile, aPoints );
+ else if ( aFileSuf == "asc" )
+ aRes = importFromASCFile( aFile, aPoints );
// Close the file
aFile.close();
return true;
}
+bool HYDROData_Bathymetry::importFromASCFile( QFile& theFile,
+ AltitudePoints& thePoints ) const
+{
+ if ( !theFile.isOpen() )
+ return false;
+
+ QString aLine;
+ QStringList aStrList;
+
+ int aNCols;
+ int aNRows;
+ double anXllCorner;
+ double anYllCorner;
+ double aCellSize;
+ double aNoDataValue;
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "ncols" )
+ return false;
+ aNCols = aStrList[1].toInt();
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "nrows" )
+ return false;
+ aNRows = aStrList[1].toInt();
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "xllcorner" )
+ return false;
+ anXllCorner = aStrList[1].toDouble();
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "yllcorner" )
+ return false;
+ anYllCorner = aStrList[1].toDouble();
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "cellsize" )
+ return false;
+ aCellSize = aStrList[1].toDouble();
+
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+ if ( aStrList.length() != 2 && aStrList[0].toLower() != "nodata_value" )
+ return false;
+ aNoDataValue = aStrList[1].toDouble();
+
+ bool anIsAltitudesInverted = IsAltitudesInverted();
+
+ int i = 0;
+ int aStrLength = 0;
+ while ( !theFile.atEnd() )
+ {
+ aLine = theFile.readLine().simplified();
+ aStrList = aLine.split( ' ', QString::SkipEmptyParts );
+
+ aStrLength = aStrList.length();
+ if ( aStrLength == 0 )
+ continue;
+
+ if ( aStrLength != aNRows )
+ return false;
+
+ for (int j = 0; j < aNCols; j++)
+ {
+ if (aStrList[j].toDouble() != aNoDataValue)
+ {
+ AltitudePoint aPoint;
+ aPoint.SetX(anXllCorner + aCellSize*(j + 0.5));
+ aPoint.SetY(anYllCorner + aCellSize*(aNRows - i + 0.5));
+ aPoint.SetZ(aStrList[j].toDouble());
+
+ if ( anIsAltitudesInverted )
+ aPoint.SetZ( -aPoint.Z() );
+
+ thePoints.Append(aPoint);
+ }
+ }
+ i++;
+
+ }
+
+ return true;
+
+}
+
Handle_HYDROData_PolylineXY HYDROData_Bathymetry::CreateBoundaryPolyline() const
{
aResult->AddPoint( 0, HYDROData_IPolyline::Point( Xmin, Ymax ) );
aResult->AddPoint( 0, HYDROData_IPolyline::Point( Xmax, Ymax ) );
aResult->AddPoint( 0, HYDROData_IPolyline::Point( Xmax, Ymin ) );
+
+ aResult->SetWireColor( HYDROData_PolylineXY::DefaultWireColor() );
+
aResult->Update();
return aResult;