-// Copyright (C) 2007-2015 CEA/DEN, EDF R&D, OPEN CASCADE
-//
-// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
-// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
-//
+// 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
#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();
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();
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;