Coordinate ConventionalCoordinateSystem::coord2coord(const ICoordinateSystem &sourceCs, const Coordinate &crdSource) const
{
if (sourceCs->isEqual(*this))
return crdSource;
LatLon ll = sourceCs->coord2latlon(crdSource);
if ( ll.isValid()) {
return latlon2coord(ll);
}
return Coordinate();
}
Coordinate ConventionalCoordinateSystem::coord2coord(const ICoordinateSystem &sourceCs, const Coordinate &crdSource) const
{
//TODO isLatLon guard doesn't consider latlon cs other than WGS84!
if (sourceCs->id() == id()) // the 'real'isEqual test is too expensive here, as this method can be called 100000's of times (resample)
return crdSource;
LatLon ll = sourceCs->isLatLon() ? LatLon(crdSource.y,crdSource.x) : sourceCs->coord2latlon(crdSource);
if ( ll.isValid()) {
return isLatLon() ? ll : latlon2coord(ll);
}
return Coordinate();
}
Coordinate ConventionalCoordinateSystem::coord2coord(const ICoordinateSystem &sourceCs, const Coordinate &crdSource) const
{
//TODO isLatLon guard doesn't consider latlon cs other than WGS84!
if (sourceCs->id() == id()) // the 'real'isEqual test is too expensive here, as this method can be called 100000's of times (resample)
return crdSource;
LatLon ll = sourceCs->isLatLon() ? LatLon(crdSource.y,crdSource.x) : sourceCs->coord2latlon(crdSource);
if (ll.isValid()) {
if (hasType(sourceCs->ilwisType(), itCONVENTIONALCOORDSYSTEM)) {
const IConventionalCoordinateSystem & srcCs = sourceCs.as<ConventionalCoordinateSystem>();
if (srcCs->datum().get() && datum().get() && !srcCs->datum()->equal(*datum().get())) { // different datums given, datum shift needed
ll = srcCs->datum()->llToWGS84(ll, *srcCs->ellipsoid().ptr());
ll = datum()->llFromWGS84(ll, *ellipsoid().ptr());
}
}
return isLatLon() ? ll : latlon2coord(ll);
}
return Coordinate();
}
void SubGridLayer::calcEnvelope(Coordinate& cmin, Coordinate& cmax) const {
cmin = layerManager()->rootLayer()->zoomEnvelope().min_corner();
cmax = layerManager()->rootLayer()->zoomEnvelope().max_corner();
if (!cmin.isValid() || !cmax.isValid())
return;
ICoordinateSystem screenCsy = layerManager()->rootLayer()->screenCsy();
if (!screenCsy.isValid())
return;
if (screenCsy->isUnknown()) // cant do any extra calcs with unknown
return;
bool isLatLon = screenCsy->isLatLon();
if (isLatLon) {
if (cmin.y <= -89)
cmin.y = -89;
if (cmin.x < -180)
cmin.x = -180;
if (cmax.y >= 89)
cmax.y = 89;
if (cmax.x > 180)
cmax.x = 180;
}
else {
LatLon llmin = screenCsy->coord2latlon(cmin);
LatLon llmax = screenCsy->coord2latlon(cmax);
if (llmin.isValid() && llmax.isValid()) {
if (llmin.Lat() <= -85)
cmin.y = screenCsy->latlon2coord(LatLon(-85, llmin.Lon())).y;
if (llmin.Lon() < -180)
cmin.x = screenCsy->latlon2coord(LatLon(llmin.Lat(), -180)).x;
if (llmax.Lat() > 85)
cmax.y = screenCsy->latlon2coord(LatLon(85, llmax.Lon())).y;
if (llmax.Lon() > 180)
cmax.x = screenCsy->latlon2coord(LatLon(llmax.Lat(), 180)).x;
}
}
}