+// 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 "HYDROData_Tool.h"
#include "HYDROData_PolylineXY.h"
-#include <boost/math/special_functions/fpclassify.hpp>
-
#include <gp_XY.hxx>
#include <gp_XYZ.hxx>
#include <TDataStd_AsciiString.hxx>
#include <TDataStd_Integer.hxx>
+#include <QColor>
#include <QFile>
#include <QFileInfo>
#include <QPointF>
{
}
-QStringList HYDROData_Bathymetry::DumpToPython( MapOfTreatedObjects& theTreatedObjects ) const
+QStringList HYDROData_Bathymetry::DumpToPython( const QString& thePyScriptPath,
+ MapOfTreatedObjects& theTreatedObjects ) const
{
QStringList aResList = dumpObjectCreation( theTreatedObjects );
QString aBathymetryName = GetObjPyName();
- aResList << QString( "%1.SetAltitudesInverted( %2 );" )
+ aResList << QString( "%1.SetAltitudesInverted( %2 )" )
.arg( aBathymetryName ).arg( IsAltitudesInverted() );
TCollection_AsciiString aFilePath = GetFilePath();
- aResList << QString( "%1.ImportFromFile( \"%2\" );" )
+ aResList << QString( "%1.ImportFromFile( \"%2\" )" )
.arg( aBathymetryName ).arg( aFilePath.ToCString() );
aResList << QString( "" );
- aResList << QString( "%1.Update();" ).arg( aBathymetryName );
+ aResList << QString( "%1.Update()" ).arg( aBathymetryName );
aResList << QString( "" );
return aResList;
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 );
}
if ( !aLabel.IsNull() )
{
aLabel.ForgetAllAttributes();
- SetToUpdate( true );
+ Changed( Geom_Z );
}
}
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();
if ( !isXOk || !isYOk || !isZOk )
return false;
- if ( boost::math::isnan( aPoint.X() ) || boost::math::isinf( aPoint.X() ) ||
- boost::math::isnan( aPoint.Y() ) || boost::math::isinf( aPoint.Y() ) ||
- boost::math::isnan( aPoint.Z() ) || boost::math::isinf( aPoint.Z() ) )
+ if ( HYDROData_Tool::IsNan( aPoint.X() ) || HYDROData_Tool::IsInf( aPoint.X() ) ||
+ HYDROData_Tool::IsNan( aPoint.Y() ) || HYDROData_Tool::IsInf( aPoint.Y() ) ||
+ HYDROData_Tool::IsNan( aPoint.Z() ) || HYDROData_Tool::IsInf( aPoint.Z() ) )
return false;
// Invert the z value if requested
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;
}
+
+void HYDROData_Bathymetry::UpdateLocalCS( double theDx, double theDy )
+{
+ gp_XYZ aDelta( theDx, theDy, 0 );
+ AltitudePoints aPoints = GetAltitudePoints();
+ AltitudePoints::Iterator anIter( aPoints );
+ for ( int i = 0 ; anIter.More(); ++i, anIter.Next() )
+ {
+ AltitudePoint& aPoint = anIter.ChangeValue();
+ aPoint += aDelta;
+ }
+ SetAltitudePoints( aPoints );
+}
+
