public void Dispose()
{
if (pjNative != IntPtr.Zero)
{
Proj6Native.proj_destroy(pjNative);
}
pjNative = IntPtr.Zero;
}
public CoordinateTransformation(ICRS source, ICRS target)
{
this.SourceCRS = source;
this.TargetCRS = target;
//shoud we create a threading context here instead of the default IntPtr.Zero?
pjNative = Proj6Native.proj_create_crs_to_crs(IntPtr.Zero, source.WKT, target.WKT, IntPtr.Zero);
if (pjNative == IntPtr.Zero)
{
throw new Exception("Could not create coordinate transformation");
}
IntPtr pn = Proj6Native.proj_normalize_for_visualization(IntPtr.Zero, pjNative);
Proj6Native.proj_destroy(pjNative);
pjNative = IntPtr.Zero;
if (pn == IntPtr.Zero)
{
throw new Exception("Could not create coordinate transformation - proj_normalize_for_visualization returned zero");
}
pjNative = pn;
}
public bool IsEquivalent(ICRS other)
{
if (other == null)
{
return(false);
}
IntPtr pjThis = Proj6Native.proj_create(IntPtr.Zero, this.WKT);
IntPtr pjOther = Proj6Native.proj_create(IntPtr.Zero, other.WKT);
try
{
if (pjThis != IntPtr.Zero && pjOther != IntPtr.Zero)
{
bool same = Proj6Native.proj_is_equivalent_to(pjThis, pjOther, Proj6Native.PJ_COMPARISON_CRITERION.PJ_COMP_EQUIVALENT_EXCEPT_AXIS_ORDER_GEOGCRS) != 0;
if (!same)
{
//proj_is_equivalent_to doesn't seem to compare different WKT representations
//convert both to ESRI WKT and compare
string wkt1 = Proj6Native.Proj_as_wkt(IntPtr.Zero, pjThis, Proj6Native.PJ_WKT_TYPE.PJ_WKT1_ESRI);
string wkt2 = Proj6Native.Proj_as_wkt(IntPtr.Zero, pjOther, Proj6Native.PJ_WKT_TYPE.PJ_WKT1_ESRI);
same = string.Compare(wkt1, wkt2, StringComparison.OrdinalIgnoreCase) == 0;
if (!same)
{
IntPtr pjWkt1 = Proj6Native.proj_create(IntPtr.Zero, wkt1);
IntPtr pjWkt2 = Proj6Native.proj_create(IntPtr.Zero, wkt2);
try
{
if (pjWkt1 != IntPtr.Zero && pjWkt2 != IntPtr.Zero)
{
same = Proj6Native.proj_is_equivalent_to(pjWkt1, pjWkt2, Proj6Native.PJ_COMPARISON_CRITERION.PJ_COMP_EQUIVALENT_EXCEPT_AXIS_ORDER_GEOGCRS) != 0;
}
}
finally
{
if (pjWkt1 != IntPtr.Zero)
{
Proj6Native.proj_destroy(pjWkt1);
}
if (pjWkt2 != IntPtr.Zero)
{
Proj6Native.proj_destroy(pjWkt2);
}
}
}
}
return(same);
}
}
finally
{
if (pjThis != IntPtr.Zero)
{
Proj6Native.proj_destroy(pjThis);
}
if (pjOther != IntPtr.Zero)
{
Proj6Native.proj_destroy(pjOther);
}
}
return(false);
}
public static CRS FromWKT(string wkt, bool identify = false)
{
//although using proj_context_create will make proj6 thread safe
//it is expensive (10x). Faster to use .net lock and just use the default context
IntPtr context = IntPtr.Zero;// Proj6Native.proj_context_create();
try
{
IntPtr p = Proj6Native.proj_create(context, wkt);
//this code returns null?
//IntPtr p = Proj6Native.Proj_create_from_wkt(IntPtr.Zero, wkt);
if (p != IntPtr.Zero && identify)
{
string name = Proj6Native.GetAuthName(p);
int confidence = 0;
IntPtr p2 = Proj6Native.Proj_identify(context, p, name, out confidence);
if (p2 != IntPtr.Zero && confidence > 85)
{
Proj6Native.proj_destroy(p);
p = p2;
}
else
{
if (p2 != IntPtr.Zero)
{
Proj6Native.proj_destroy(p2);
}
}
}
try
{
Proj6Native.PJ_TYPE pType = Proj6Native.proj_get_type(p);
//Console.Out.WriteLine("pType=" + pType);
string name = Proj6Native.GetName(p);
string authName = Proj6Native.GetAuthName(p);
string id = Proj6Native.ProjGetIdCode(p);
CRSBoundingBox areaOfUse = new CRSBoundingBox()
{
WestLongitudeDegrees = -1000,
NorthLatitudeDegrees = -1000,
EastLongitudeDegrees = -1000,
SouthLatitudeDegrees = -1000
};
Proj6Native.proj_get_area_of_use(context, p,
ref areaOfUse.WestLongitudeDegrees,
ref areaOfUse.SouthLatitudeDegrees,
ref areaOfUse.EastLongitudeDegrees,
ref areaOfUse.NorthLatitudeDegrees,
IntPtr.Zero);
if (identify)
{
string axisName;
string axisAbbrev;
string axisDirection;
double unit_conv_factor = 1;
string unit_name;
string unit_auth_name;
string unit_code;
//int axisCount = Proj6Native.proj_cs_get_axis_count(IntPtr.Zero, p);
//if (axisCount > 0)
//{
// if (Proj6Native.Proj_cs_get_axis_info(IntPtr.Zero, p, 0, out axisName, out axisAbbrev, out axisDirection, out unit_conv_factor, out unit_name, out unit_auth_name, out unit_code))
// {
// Console.Out.WriteLine(axisName);
// }
//}
IntPtr pCrs = Proj6Native.proj_crs_get_coordinate_system(context, p);
if (pCrs != IntPtr.Zero)
{
int axisCount = Proj6Native.proj_cs_get_axis_count(context, pCrs);
if (axisCount > 0)
{
if (Proj6Native.Proj_cs_get_axis_info(context, pCrs, 0, out axisName, out axisAbbrev, out axisDirection, out unit_conv_factor, out unit_name, out unit_auth_name, out unit_code))
{
//System.Diagnostics.Debug.WriteLine(axisName);
// System.Diagnostics.Debug.WriteLine("unit_conv_factor:" + unit_conv_factor);
}
}
Proj6Native.proj_destroy(pCrs);
}
}
if (pType == Proj6Native.PJ_TYPE.PJ_TYPE_GEOGRAPHIC_2D_CRS)
{
return(new GeographicCRS()
{
Id = id,
Name = name,
Authority = authName,
WKT = wkt,
AreaOfUse = areaOfUse
});
}
else if (pType == Proj6Native.PJ_TYPE.PJ_TYPE_PROJECTED_CRS)
{
return(new ProjectedCRS()
{
Id = id,
Name = name,
Authority = authName,
WKT = wkt,
UnitsToMeters = 1,
AreaOfUse = areaOfUse
});
}
else if (pType == Proj6Native.PJ_TYPE.PJ_TYPE_BOUND_CRS)
{
if (wkt.IndexOf("PROJECTION", StringComparison.OrdinalIgnoreCase) >= 0)
{
return(new ProjectedCRS()
{
Id = id,
Name = name,
Authority = authName,
WKT = wkt,
UnitsToMeters = 1,
AreaOfUse = areaOfUse
});
}
else
{
return(new GeographicCRS()
{
Id = id,
Name = name,
Authority = authName,
WKT = wkt,
AreaOfUse = areaOfUse
});
}
}
else
{
//Console.Out.WriteLine("pType = " + pType);
}
return(null);
}
finally
{
if (p != IntPtr.Zero)
{
Proj6Native.proj_destroy(p);
}
}
}
finally
{
if (context != IntPtr.Zero)
{
Proj6Native.proj_context_destroy(context);
}
}
}