string pointTransformAndWriteZ(string geom, MgMap map)
{
double[] x = new double[1];
double[] y = new double[1];
double[] z = new double[1];
MgWktReaderWriter wktrwdmr = new MgWktReaderWriter();
MgPoint cntr = wktrwdmr.Read(geom).Centroid;
x[0] = cntr.Coordinate.X;
y[0] = cntr.Coordinate.Y;
z[0] = 0;
StringBuilder output = new StringBuilder();
output.Append("<table width=\"100%\" class=\"results\">");
//WGS preko MG API
MgCoordinateSystemFactory fact = new MgCoordinateSystemFactory();
string wktFrom = map.GetMapSRS();
string wktTo = "GEOGCS[\"LL84\",DATUM[\"WGS84\",SPHEROID[\"WGS84\",6378137.000,298.25722293]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.01745329251994]]";
MgCoordinateSystem CSSource = fact.Create(wktFrom);
MgCoordinateSystem CSTarget = fact.Create(wktTo);
MgCoordinateSystemTransform coordTransform = fact.GetTransform(CSSource, CSTarget);
MgGeometryFactory geomFact = new MgGeometryFactory();
// GK
//samo prvo višino zaenkrat, dokler ni enačbe za Z
output.Append(String.Format("<tr><td class='header'><b>Map koordinates:</b></td><td class=\"results\"><table><tr><td><b>Y:</b></td><td>{0}</td></tr><tr><td><b>X:</b></td><td>{1}</td></tr></table></td></tr>", string.Format("{0:0.0}", x[0]), string.Format("{0:0.0}", y[0])));
int i = 0;
//transformacija preko MG koordinatnega sistema
foreach (double pointX in x)
{
MgCoordinate coord = geomFact.CreateCoordinateXY(x[i], y[i]);
coord = coordTransform.Transform(coord);
x[i] = coord.X;
y[i] = coord.Y;
i++;
}
double[] xwgs = x;
double[] ywgs = y;
output.Append(String.Format("<tr><td class='header'><b>WGS84:</b></td><td class=\"results\"><table><tr><td><b>Lon:</b></td><td>{0}</td></tr><tr><td><b>Lat:</b></td><td>{1}</td></tr></table></td></tr>", string.Format("{0:0.000000}", xwgs[0]), string.Format("{0:0.000000}", ywgs[0])));
output.Append("</table>");
return(output.ToString());
}
// Create a temporary Feature Source to store geocode results.
public void CreateLocationMarkerFeatureSource(MgFeatureService featureService, MgResourceIdentifier locationMarkerDataResId, MgMap map)
{
//String ll84Wkt = "GEOGCS[\"LL84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.25722293287],TOWGS84[0,0,0,0,0,0,0]],PRIMEM[\"Greenwich\",0],UNIT[\"Degrees\",1]]";
//string ll84Wkt = "LOCAL_CS[\"*XY-MT*\",LOCAL_DATUM[\"*X-Y*\",10000],UNIT[\"Meter\", 1],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]] ";
MgCoordinateSystemFactory sysCoordsFactory = new MgCoordinateSystemFactory();
string srsDefMap = map.GetMapSRS();
string mapSrsUnits = string.Empty;
srsDefMap = !string.IsNullOrEmpty(srsDefMap) ? srsDefMap : "LOCAL_CS[\"*XY-MT*\",LOCAL_DATUM[\"*X-Y*\",10000],UNIT[\"Meter\", 1],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]] ";
MgCoordinateSystem srsMap = sysCoordsFactory.Create(srsDefMap);
MgClassDefinition locationClass = new MgClassDefinition();
locationClass.SetName("LocationMarker");
MgPropertyDefinitionCollection properties = locationClass.GetProperties();
MgDataPropertyDefinition idProperty = new MgDataPropertyDefinition("ID");
idProperty.SetDataType(MgPropertyType.Int32);
idProperty.SetReadOnly(true);
idProperty.SetNullable(false);
idProperty.SetAutoGeneration(true);
properties.Add(idProperty);
MgDataPropertyDefinition addressProperty = new MgDataPropertyDefinition("Address");
addressProperty.SetDataType(MgPropertyType.String);
addressProperty.SetLength(512);
properties.Add(addressProperty);
MgGeometricPropertyDefinition locationProperty = new MgGeometricPropertyDefinition("Location");
locationProperty.SetGeometryTypes(MgGeometryType.Point);
locationProperty.SetHasElevation(false);
locationProperty.SetHasMeasure(false);
locationProperty.SetReadOnly(false);
locationProperty.SetSpatialContextAssociation(srsMap.GetCsCode());
properties.Add(locationProperty);
MgPropertyDefinitionCollection idProperties = locationClass.GetIdentityProperties();
idProperties.Add(idProperty);
locationClass.SetDefaultGeometryPropertyName("Location");
MgFeatureSchema locationSchema = new MgFeatureSchema();
locationSchema.SetName("LocationMarkerSchema");
locationSchema.GetClasses().Add(locationClass);
MgCreateSdfParams sdfParams = new MgCreateSdfParams(srsMap.GetCsCode(), srsDefMap, locationSchema);
featureService.CreateFeatureSource(locationMarkerDataResId, sdfParams);
}
public static string TransformGK2Wgs(double x, double y, string ss, string mpN)
{
MgMap map = new MgMap();
MgResourceService resourceSrvc = GetMgResurceService(ss);
map.Open(resourceSrvc, mpN);
//Create coordinate system factory
MgCoordinateSystemFactory fact = new MgCoordinateSystemFactory();
string wktFrom = map.GetMapSRS();
string wktTo = "GEOGCS[\"LL84\",DATUM[\"WGS84\",SPHEROID[\"WGS84\",6378137.000,298.25722293]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.01745329251994]]";
MgCoordinateSystem coordinateSystemSource = fact.Create(wktFrom);
MgCoordinateSystem coordinateSystemTarget = fact.Create(wktTo);
MgGeometryFactory geomFact = new MgGeometryFactory();
MgCoordinateSystemTransform coordTransform = fact.GetTransform(coordinateSystemSource, coordinateSystemTarget);
MgCoordinate coord = coordTransform.Transform(x, y);
return coord.X.ToString().Replace(',', '.') + ";" + coord.Y.ToString().Replace(',', '.');
}
public static string TransformGK2Wgs(double x, double y, string ss, string mpN) //Transform from Map to WGS 84 Coordtnate System
{
MgMap map = new MgMap();
MgResourceService resourceSrvc = GetMgResurceService(ss);
map.Open(resourceSrvc, mpN);
//Create coordinate system factory
MgCoordinateSystemFactory fact = new MgCoordinateSystemFactory();
string wktFrom = map.GetMapSRS();
string wktTo = "GEOGCS[\"LL84\",DATUM[\"WGS84\",SPHEROID[\"WGS84\",6378137.000,298.25722293]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.01745329251994]]";
MgCoordinateSystem coordinateSystemSource = fact.Create(wktFrom);
MgCoordinateSystem coordinateSystemTarget = fact.Create(wktTo);
MgGeometryFactory geomFact = new MgGeometryFactory();
MgCoordinateSystemTransform coordTransform = fact.GetTransform(coordinateSystemSource, coordinateSystemTarget);
MgCoordinate coord = coordTransform.Transform(x, y);
return(coord.X.ToString().Replace(',', '.') + ";" + coord.Y.ToString().Replace(',', '.'));
}
public bool ShowSpatialFilter()
{
bool result = true;
MgUserInformation userInfo = new MgUserInformation(Request["SESSION"]);
MgSiteConnection siteConnection = new MgSiteConnection();
siteConnection.Open(userInfo);
MgResourceIdentifier sdfResId = new MgResourceIdentifier("Session:" + Request["SESSION"] + "//Filter.FeatureSource");
MgResourceService resourceService = siteConnection.CreateService(MgServiceType.ResourceService) as MgResourceService;
MgFeatureService featureService = siteConnection.CreateService(MgServiceType.FeatureService) as MgFeatureService;
MgFeatureCommandCollection updateCommands = new MgFeatureCommandCollection();
MgMap map = new MgMap();
map.Open(resourceService, Request["MAPNAME"]);
MgLayer layer = null;
MgLayerCollection layers = map.GetLayers();
if (layers.Contains("_QuerySpatialFilter"))
{
layer = (MgLayer)layers.GetItem("_QuerySpatialFilter");
//updateCommands.Add(new MgDeleteFeatures("Filter", "ID > 0"));
}
else
{
// Create the Feature Source (SDF)
MgFeatureSchema sdfSchema = this.CreateFilterSchema();
MgCreateSdfParams sdfParams = new MgCreateSdfParams("MAPCS", map.GetMapSRS(), sdfSchema);
featureService.CreateFeatureSource(sdfResId, sdfParams);
// Create the Layer
MgResourceIdentifier layerResId = new MgResourceIdentifier("Session:" + Request["SESSION"] + "//Filter.LayerDefinition");
String layerDefinition = File.ReadAllText(GetQueryXmlTemplatePath());
layerDefinition = layerDefinition.Replace("%s", sdfResId.ToString());
MgByteReader reader = new MgByteReader(layerDefinition, "text/xml");
resourceService.SetResource(layerResId, reader, null);
layer = new MgLayer(layerResId, resourceService);
layer.SetName("_QuerySpatialFilter");
layer.SetLegendLabel("תחום זמני");
layer.SetDisplayInLegend(true);
layer.SetSelectable(true);
layer.ForceRefresh();
layer.NeedsRefresh();
layers.Insert(0, layer);
}
// Make the layer visible
layer.SetVisible(true);
map.Save(resourceService);
// Add the geometry to the filter feature source
MgPolygon polygon = this.CreatePolygonFromGeomText(Request["GEOMTEXT"].ToString());
MgAgfReaderWriter agfWriter = new MgAgfReaderWriter();
MgByteReader byteReader = agfWriter.Write(polygon);
MgPropertyCollection propertyValues = new MgPropertyCollection();
propertyValues.Add(new MgGeometryProperty("Geometry", byteReader));
try
{
updateCommands.Add(new MgInsertFeatures("Filter", propertyValues));
featureService.UpdateFeatures(sdfResId, updateCommands, false);
}
catch { }
return result;
}
//----------------------------------------------------------------------------------------
// �� �ܣ� ��ȡ��ͼ������ϵͳ
//
// �� �ߣ�
//
//
// �� �ڣ�2007.05.#
//
//-----------------------------------------------------------------------------------------
public String GetMapSrs(MgMap map)
{
String srs = map.GetMapSRS();
if (srs != "")
return srs;
return "LOCALCS[\"*XY-MT*\",LOCAL_DATUM[\"*X-Y*\",10000],UNIT[\"Meter\", 1],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]]";
}
public MgCoordinateSystemTransform GetTransform()
{
MgCoordinateSystemFactory coordSysFactory = new MgCoordinateSystemFactory();
_resourceService = (MgResourceService)this.site.CreateService(MgServiceType.ResourceService);
MgMap map = new MgMap();
map.Open(_resourceService, GetParameter(this.args, "MAPNAME"));
MgCoordinateSystem source = coordSysFactory.Create(map.GetMapSRS());
//MgCoordinateSystem target = coordSysFactory.Create(MarkupManager.LL84WKT);
MgCoordinateSystem target = coordSysFactory.Create(map.GetMapSRS());
return coordSysFactory.GetTransform(source, target);
}
string pointTransformAndWriteZ(string geom, MgMap map)
{
double[] x = new double[1];
double[] y = new double[1];
double[] z = new double[1];
MgWktReaderWriter wktrwdmr = new MgWktReaderWriter();
MgPoint cntr = wktrwdmr.Read(geom).Centroid;
x[0] = cntr.Coordinate.X;
y[0] = cntr.Coordinate.Y;
z[0] = 0;
StringBuilder output = new StringBuilder();
output.Append("<table width=\"100%\" class=\"results\">");
//WGS preko MG API
MgCoordinateSystemFactory fact = new MgCoordinateSystemFactory();
string wktFrom = map.GetMapSRS();
string wktTo = "GEOGCS[\"LL84\",DATUM[\"WGS84\",SPHEROID[\"WGS84\",6378137.000,298.25722293]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.01745329251994]]";
MgCoordinateSystem CSSource = fact.Create(wktFrom);
MgCoordinateSystem CSTarget = fact.Create(wktTo);
MgCoordinateSystemTransform coordTransform = fact.GetTransform(CSSource, CSTarget);
MgGeometryFactory geomFact = new MgGeometryFactory();
// GK
//samo prvo višino zaenkrat, dokler ni enačbe za Z
output.Append(String.Format("<tr><td class='header'><b>Map koordinates:</b></td><td class=\"results\"><table><tr><td><b>Y:</b></td><td>{0}</td></tr><tr><td><b>X:</b></td><td>{1}</td></tr></table></td></tr>", string.Format("{0:0.0}", x[0]), string.Format("{0:0.0}", y[0])));
int i = 0;
//transformacija preko MG koordinatnega sistema
foreach (double pointX in x)
{
MgCoordinate coord = geomFact.CreateCoordinateXY(x[i], y[i]);
coord = coordTransform.Transform(coord);
x[i] = coord.X;
y[i] = coord.Y;
i++;
}
double[] xwgs = x;
double[] ywgs = y;
output.Append(String.Format("<tr><td class='header'><b>WGS84:</b></td><td class=\"results\"><table><tr><td><b>Lon:</b></td><td>{0}</td></tr><tr><td><b>Lat:</b></td><td>{1}</td></tr></table></td></tr>", string.Format("{0:0.000000}", xwgs[0]), string.Format("{0:0.000000}", ywgs[0])));
output.Append("</table>");
return output.ToString();
}
public static string GetMapSrs(MgMap map)
{
string srs = map.GetMapSRS();
if (srs != "")
{
return srs;
}
return "LOCALCS[\"Non-Earth (Meter)\",LOCAL_DATUM[\"Local Datum\",0],UNIT[\"Meter\", 1],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]]";
}
