/// <summary>
/// Coordinate transformation between two world coordinate systems
/// This Version works with coordinates as X, Y doubles
/// </summary>
/// <param name="wkTextSource">The wk text source.</param>
/// <param name="wkTextTarget">The wk text target.</param>
/// <param name="inputPointX">The input point X.</param>
/// <param name="inputPointY">The input point Y.</param>
/// <param name="transformedPointX">The transformed point X.</param>
/// <param name="transformedPointY">The transformed point Y.</param>
/// <returns></returns>
public static bool TransformPoint2d(string wkTextSource, string wkTextTarget, double inputPointX, double inputPointY, ref double transformedPointX, ref double transformedPointY)
{
bool RetVal = false;
//_logger.Debug("Start TransformPoint2d");
try
{
//Creating coordinate system factory
MgCoordinateSystemFactory CSFactory = new MgCoordinateSystemFactory();
//Creating coordinate system objects
MgCoordinateSystem SourceCS = CSFactory.Create(wkTextSource);
MgCoordinateSystem TargetCS = CSFactory.Create(wkTextTarget);
//Creating geometry factory
MgGeometryFactory GeomFactory = new MgGeometryFactory();
//Populating geometry factory CreateCoordinateXY
MgCoordinate SourceCoord = GeomFactory.CreateCoordinateXY(inputPointX, inputPointY);
//Getting transformation definition
MgCoordinateSystemTransform CSTransform = CSFactory.GetTransform(SourceCS, TargetCS);
//Transforming coordinate
MgCoordinate TargetCoord = CSTransform.Transform(SourceCoord);
//Populating return coordinate objects
transformedPointX = TargetCoord.X;
transformedPointY = TargetCoord.Y;
RetVal = true;
}
catch (System.Exception ex)
{
//_logger.Error("Error in TransformPoint2d", ex);
throw;
}
//_logger.Debug("End TransformPoint2d");
return(RetVal);
}
internal LocalNativeSimpleTransform(string sourceWkt, string targetWkt)
{
var fact = new MgCoordinateSystemFactory();
var source = fact.Create(sourceWkt);
var target = fact.Create(targetWkt);
_trans = fact.GetTransform(source, target);
}
internal LocalNativeSimpleTransform(string sourceWkt, string targetWkt)
{
var fact = new MgCoordinateSystemFactory();
var source = fact.Create(sourceWkt);
var target = fact.Create(targetWkt);
_trans = fact.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 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(',', '.'));
}
private void btnCreate_Click(object sender, EventArgs e)
{
btnCreate.Enabled = false;
try
{
var layerName = txtBufferLayer.Text.Trim();
if (string.IsNullOrEmpty(layerName))
{
MessageBox.Show(Strings.MsgEnterNameForLayer);
return;
}
if (lstLayers.SelectedItems.Count == 0)
{
MessageBox.Show(Strings.MsgIncludeLayersToBuffer);
return;
}
var map = _viewer.GetMap();
var layers = map.GetLayers();
var provider = _viewer.GetProvider();
//From here, it's the same logic as buffer.aspx in .net MapGuide AJAX viewer
MgResourceIdentifier fsId = new MgResourceIdentifier("Session:" + _sessionId + "//" + txtBufferLayer.Text + "_Buffer.FeatureSource"); //NOXLATE
MgResourceIdentifier ldfId = new MgResourceIdentifier("Session:" + _sessionId + "//" + txtBufferLayer.Text + "_Buffer.LayerDefinition"); //NOXLATE
MgLayerBase layer = Util.FindLayer(layers, txtBufferLayer.Text);
string[] layerNames = GetLayerNames();
double distance = Convert.ToDouble(numBufferDistance.Value);
MeasurementUnit bUnits = (MeasurementUnit)cmbUnits.SelectedItem;
switch (bUnits)
{
case MeasurementUnit.Feet:
distance *= 0.30480;
break;
case MeasurementUnit.Kilometers:
distance *= 1000;
break;
case MeasurementUnit.Miles:
distance *= 1609.35;
break;
}
String srsDefMap = Util.GetMapSrs(map);
MgCoordinateSystem srsMap = provider.GetMapCoordinateSystem();
string mapSrsUnits = "";
bool arbitraryMapSrs = (srsMap.GetType() == MgCoordinateSystemType.Arbitrary);
if (arbitraryMapSrs)
mapSrsUnits = srsMap.GetUnits();
String xtrans = String.Format("{0:x2}", ((int)(255 * Convert.ToInt32(numFillTransparency.Value) / 100))); //NOXLATE
var lineColor = Util.ToHtmlColor(pnlBorderColor.BackColor);
var foreColor = Util.ToHtmlColor(pnlFillColor.BackColor);
var backColor = Util.ToHtmlColor(pnlFillBackColor.BackColor);
String layerTempl = string.Format(Properties.Resources.AreaLayerDef,
fsId.ToString(),
"BufferSchema:Buffer", //NOXLATE
"GEOM", //NOXLATE
cmbFillPattern.SelectedItem,
xtrans + foreColor,
((0 != 1/*transparent*/) ? "ff" : "00") + backColor, //NOXLATE
cmbBorderPattern.SelectedItem,
numLineThickness.Value.ToString(NumberFormatInfo.InvariantInfo),
lineColor
);
byte[] bytes = Encoding.UTF8.GetBytes(layerTempl);
MgByteSource src = new MgByteSource(bytes, bytes.Length);
MgByteReader layerDefContent = src.GetReader();
_resSvc.SetResource(ldfId, layerDefContent, null);
bool newBuffer = false;
if (layer == null)
{
newBuffer = true;
//Targetting a new layer. create a data source for it
//
MgClassDefinition classDef = new MgClassDefinition();
classDef.SetName("Buffer"); //NOXLATE
classDef.SetDescription("Feature class for buffer layer"); //NOXLATE
classDef.SetDefaultGeometryPropertyName("GEOM"); //NOXLATE
//Set KEY property
MgDataPropertyDefinition prop = new MgDataPropertyDefinition("KEY"); //NOXLATE
prop.SetDataType(MgPropertyType.Int32);
prop.SetAutoGeneration(true);
prop.SetReadOnly(true);
classDef.GetIdentityProperties().Add(prop);
classDef.GetProperties().Add(prop);
//Set ID property. Hold this segment ID
prop = new MgDataPropertyDefinition("ID"); //NOXLATE
prop.SetDataType(MgPropertyType.Int32);
classDef.GetProperties().Add(prop);
//Set geometry property
MgGeometricPropertyDefinition geomProp = new MgGeometricPropertyDefinition("GEOM"); //NOXLATE
//prop.SetGeometryTypes(MgFeatureGeometricType.mfgtSurface); //TODO use the constant when exposed
geomProp.SetGeometryTypes(4);
classDef.GetProperties().Add(geomProp);
//Create the schema
MgFeatureSchema schema = new MgFeatureSchema("BufferSchema", "Temporary buffer schema"); //NOXLATE
schema.GetClasses().Add(classDef);
//finally, creation of the feature source
MgFileFeatureSourceParams sdfParams = new MgFileFeatureSourceParams("OSGeo.SDF", "LatLong", map.GetMapSRS(), schema); //NOXLATE
_featSvc.CreateFeatureSource(fsId, sdfParams);
//Add layer to map
layer = provider.CreateLayer(ldfId);
layer.SetName(txtBufferLayer.Text);
layer.SetLegendLabel(txtBufferLayer.Text);
layer.SetDisplayInLegend(true);
layer.SetSelectable(true);
layers.Insert(0, layer);
}
else
{
//data source already exist. clear its content
//
Util.ClearDataSource(_featSvc, fsId, "BufferSchema:Buffer"); //NOXLATE
}
var sel = _viewer.GetSelection();
var selLayers = sel.GetLayers();
MgAgfReaderWriter agfRW = new MgAgfReaderWriter();
MgGeometryCollection bufferGeometries = new MgGeometryCollection();
MgGeometry geomBuffer;
MgFeatureCommandCollection commands = new MgFeatureCommandCollection();
int featId = 0;
MgBatchPropertyCollection propCollection = new MgBatchPropertyCollection();
int excludedLayers = 0;
MgCoordinateSystem srsDs = null;
MgGeometryCollection inputGeometries = new MgGeometryCollection();
int bufferFeatures = 0;
for (int li = 0; li < selLayers.GetCount(); li++)
{
MgLayerBase selLayer = selLayers.GetItem(li);
bool inputLayer = false;
String selLayerName = selLayer.GetName();
for (int il = 0; il < layerNames.Length; il++)
{
if (layerNames[il].Equals(selLayerName))
{
inputLayer = true;
break;
}
}
if (inputLayer == false)
{
continue;
}
// get the data source SRS
//
MgResourceIdentifier featSourceId = new MgResourceIdentifier(selLayer.GetFeatureSourceId());
MgSpatialContextReader ctxs = _featSvc.GetSpatialContexts(featSourceId, false);
String srsDefDs = string.Empty;
if (ctxs != null && ctxs.ReadNext())
srsDefDs = ctxs.GetCoordinateSystemWkt();
if (srsDefDs == null || srsDefDs.Length == 0)
{
excludedLayers++;
continue;
}
var srsFactory = new MgCoordinateSystemFactory();
srsDs = srsFactory.Create(srsDefDs);
bool arbitraryDsSrs = (srsDs.GetType() == MgCoordinateSystemType.Arbitrary);
String dsSrsUnits = string.Empty;
if (arbitraryDsSrs)
dsSrsUnits = srsDs.GetUnits();
// exclude layer if:
// the map is non-arbitrary and the layer is arbitrary or vice-versa
// or
// layer and map are both arbitrary but have different units
//
if ((arbitraryDsSrs != arbitraryMapSrs) || (arbitraryDsSrs && (dsSrsUnits != mapSrsUnits)))
{
excludedLayers++;
continue;
}
// calculate distance in the data source SRS units
//
double dist = srsDs.ConvertMetersToCoordinateSystemUnits(distance);
// calculate great circle unless data source srs is arbitrary
MgCoordinateSystemMeasure measure;
if (!arbitraryDsSrs)
measure = srsDs.GetMeasure();
else
measure = null;
// create a SRS transformer if necessary
MgCoordinateSystemTransform srsXform;
if (!srsDefDs.Equals(srsDefMap))
srsXform = srsFactory.GetTransform(srsDs, srsMap);
else
srsXform = null;
String featureClassName = selLayer.GetFeatureClassName();
String filter = sel.GenerateFilter(selLayer, featureClassName);
if (filter == null || filter.Length == 0)
continue;
MgFeatureQueryOptions query = new MgFeatureQueryOptions();
query.SetFilter(filter);
MgResourceIdentifier featureSource = new MgResourceIdentifier(selLayer.GetFeatureSourceId());
MgFeatureReader features = _featSvc.SelectFeatures(featureSource, featureClassName, query);
if (features.ReadNext())
{
MgClassDefinition classDef = features.GetClassDefinition();
String geomPropName = classDef.GetDefaultGeometryPropertyName();
do
{
MgByteReader geomReader = features.GetGeometry(geomPropName);
MgGeometry geom = agfRW.Read(geomReader);
if (!chkMergeBuffers.Checked)
{
geomBuffer = geom.Buffer(dist, measure);
if (geomBuffer != null)
{
if (srsXform != null)
geomBuffer = (MgGeometry)geomBuffer.Transform(srsXform);
Util.AddFeatureToCollection(propCollection, agfRW, featId++, geomBuffer);
bufferFeatures++;
}
}
else
{
if (srsXform != null)
geom = (MgGeometry)geom.Transform(srsXform);
inputGeometries.Add(geom);
}
}
while (features.ReadNext());
features.Close();
}
}
if (chkMergeBuffers.Checked)
{
if (inputGeometries.GetCount() > 0)
{
double dist = srsMap.ConvertMetersToCoordinateSystemUnits(distance);
MgCoordinateSystemMeasure measure;
if (!arbitraryMapSrs)
measure = srsMap.GetMeasure();
else
measure = null;
MgGeometryFactory geomFactory = new MgGeometryFactory();
geomBuffer = geomFactory.CreateMultiGeometry(inputGeometries).Buffer(dist, measure);
if (geomBuffer != null)
{
Util.AddFeatureToCollection(propCollection, agfRW, featId, geomBuffer);
bufferFeatures = 1;
}
}
}
if (propCollection.GetCount() > 0)
{
commands.Add(new MgInsertFeatures("BufferSchema:Buffer", propCollection)); //NOXLATE
//Insert the features in the temporary data source
//
Util.ReleaseReader(_featSvc.UpdateFeatures(fsId, commands, false), commands);
}
// Save the new map state
//
layer.ForceRefresh();
_viewer.RefreshMap();
//build report message
if (newBuffer)
MessageBox.Show(string.Format(Strings.MsgBufferLayerCreated, txtBufferLayer.Text));
else
MessageBox.Show(string.Format(Strings.MsgBufferLayerUpdated, txtBufferLayer.Text));
}
finally
{
btnCreate.Enabled = true;
}
}
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();
}
private void btnCreate_Click(object sender, EventArgs e)
{
btnCreate.Enabled = false;
try
{
var layerName = txtBufferLayer.Text.Trim();
if (string.IsNullOrEmpty(layerName))
{
MessageBox.Show(Strings.MsgEnterNameForLayer);
return;
}
if (lstLayers.SelectedItems.Count == 0)
{
MessageBox.Show(Strings.MsgIncludeLayersToBuffer);
return;
}
var map = _viewer.GetMap();
var layers = map.GetLayers();
var provider = _viewer.GetProvider();
//From here, it's the same logic as buffer.aspx in .net MapGuide AJAX viewer
MgResourceIdentifier fsId = new MgResourceIdentifier("Session:" + _sessionId + "//" + txtBufferLayer.Text + "_Buffer.FeatureSource"); //NOXLATE
MgResourceIdentifier ldfId = new MgResourceIdentifier("Session:" + _sessionId + "//" + txtBufferLayer.Text + "_Buffer.LayerDefinition"); //NOXLATE
MgLayerBase layer = Util.FindLayer(layers, txtBufferLayer.Text);
string[] layerNames = GetLayerNames();
double distance = Convert.ToDouble(numBufferDistance.Value);
MeasurementUnit bUnits = (MeasurementUnit)cmbUnits.SelectedItem;
switch (bUnits)
{
case MeasurementUnit.Feet:
distance *= 0.30480;
break;
case MeasurementUnit.Kilometers:
distance *= 1000;
break;
case MeasurementUnit.Miles:
distance *= 1609.35;
break;
}
String srsDefMap = Util.GetMapSrs(map);
MgCoordinateSystem srsMap = provider.GetMapCoordinateSystem();
string mapSrsUnits = "";
bool arbitraryMapSrs = (srsMap.GetType() == MgCoordinateSystemType.Arbitrary);
if (arbitraryMapSrs)
{
mapSrsUnits = srsMap.GetUnits();
}
String xtrans = String.Format("{0:x2}", ((int)(255 * Convert.ToInt32(numFillTransparency.Value) / 100))); //NOXLATE
var lineColor = Util.ToHtmlColor(pnlBorderColor.BackColor);
var foreColor = Util.ToHtmlColor(pnlFillColor.BackColor);
var backColor = Util.ToHtmlColor(pnlFillBackColor.BackColor);
String layerTempl = string.Format(Properties.Resources.AreaLayerDef,
fsId.ToString(),
"BufferSchema:Buffer", //NOXLATE
"GEOM", //NOXLATE
cmbFillPattern.SelectedItem,
xtrans + foreColor,
((0 != 1 /*transparent*/) ? "ff" : "00") + backColor, //NOXLATE
cmbBorderPattern.SelectedItem,
numLineThickness.Value.ToString(NumberFormatInfo.InvariantInfo),
lineColor
);
byte[] bytes = Encoding.UTF8.GetBytes(layerTempl);
MgByteSource src = new MgByteSource(bytes, bytes.Length);
MgByteReader layerDefContent = src.GetReader();
_resSvc.SetResource(ldfId, layerDefContent, null);
bool newBuffer = false;
if (layer == null)
{
newBuffer = true;
//Targetting a new layer. create a data source for it
//
MgClassDefinition classDef = new MgClassDefinition();
classDef.SetName("Buffer"); //NOXLATE
classDef.SetDescription("Feature class for buffer layer"); //NOXLATE
classDef.SetDefaultGeometryPropertyName("GEOM"); //NOXLATE
//Set KEY property
MgDataPropertyDefinition prop = new MgDataPropertyDefinition("KEY"); //NOXLATE
prop.SetDataType(MgPropertyType.Int32);
prop.SetAutoGeneration(true);
prop.SetReadOnly(true);
classDef.GetIdentityProperties().Add(prop);
classDef.GetProperties().Add(prop);
//Set ID property. Hold this segment ID
prop = new MgDataPropertyDefinition("ID"); //NOXLATE
prop.SetDataType(MgPropertyType.Int32);
classDef.GetProperties().Add(prop);
//Set geometry property
MgGeometricPropertyDefinition geomProp = new MgGeometricPropertyDefinition("GEOM"); //NOXLATE
//prop.SetGeometryTypes(MgFeatureGeometricType.mfgtSurface); //TODO use the constant when exposed
geomProp.SetGeometryTypes(4);
classDef.GetProperties().Add(geomProp);
//Create the schema
MgFeatureSchema schema = new MgFeatureSchema("BufferSchema", "Temporary buffer schema"); //NOXLATE
schema.GetClasses().Add(classDef);
//finally, creation of the feature source
MgFileFeatureSourceParams sdfParams = new MgFileFeatureSourceParams("OSGeo.SDF", "LatLong", map.GetMapSRS(), schema); //NOXLATE
_featSvc.CreateFeatureSource(fsId, sdfParams);
//Add layer to map
layer = provider.CreateLayer(ldfId);
layer.SetName(txtBufferLayer.Text);
layer.SetLegendLabel(txtBufferLayer.Text);
layer.SetDisplayInLegend(true);
layer.SetSelectable(true);
layers.Insert(0, layer);
}
else
{
//data source already exist. clear its content
//
Util.ClearDataSource(_featSvc, fsId, "BufferSchema:Buffer"); //NOXLATE
}
var sel = _viewer.GetSelection();
var selLayers = sel.GetLayers();
MgAgfReaderWriter agfRW = new MgAgfReaderWriter();
MgGeometryCollection bufferGeometries = new MgGeometryCollection();
MgGeometry geomBuffer;
MgFeatureCommandCollection commands = new MgFeatureCommandCollection();
int featId = 0;
MgBatchPropertyCollection propCollection = new MgBatchPropertyCollection();
int excludedLayers = 0;
MgCoordinateSystem srsDs = null;
MgGeometryCollection inputGeometries = new MgGeometryCollection();
int bufferFeatures = 0;
for (int li = 0; li < selLayers.GetCount(); li++)
{
MgLayerBase selLayer = selLayers.GetItem(li);
bool inputLayer = false;
String selLayerName = selLayer.GetName();
for (int il = 0; il < layerNames.Length; il++)
{
if (layerNames[il].Equals(selLayerName))
{
inputLayer = true;
break;
}
}
if (inputLayer == false)
{
continue;
}
// get the data source SRS
//
MgResourceIdentifier featSourceId = new MgResourceIdentifier(selLayer.GetFeatureSourceId());
MgSpatialContextReader ctxs = _featSvc.GetSpatialContexts(featSourceId, false);
String srsDefDs = string.Empty;
if (ctxs != null && ctxs.ReadNext())
{
srsDefDs = ctxs.GetCoordinateSystemWkt();
}
if (srsDefDs == null || srsDefDs.Length == 0)
{
excludedLayers++;
continue;
}
var srsFactory = new MgCoordinateSystemFactory();
srsDs = srsFactory.Create(srsDefDs);
bool arbitraryDsSrs = (srsDs.GetType() == MgCoordinateSystemType.Arbitrary);
String dsSrsUnits = string.Empty;
if (arbitraryDsSrs)
{
dsSrsUnits = srsDs.GetUnits();
}
// exclude layer if:
// the map is non-arbitrary and the layer is arbitrary or vice-versa
// or
// layer and map are both arbitrary but have different units
//
if ((arbitraryDsSrs != arbitraryMapSrs) || (arbitraryDsSrs && (dsSrsUnits != mapSrsUnits)))
{
excludedLayers++;
continue;
}
// calculate distance in the data source SRS units
//
double dist = srsDs.ConvertMetersToCoordinateSystemUnits(distance);
// calculate great circle unless data source srs is arbitrary
MgCoordinateSystemMeasure measure;
if (!arbitraryDsSrs)
{
measure = srsDs.GetMeasure();
}
else
{
measure = null;
}
// create a SRS transformer if necessary
MgCoordinateSystemTransform srsXform;
if (!srsDefDs.Equals(srsDefMap))
{
srsXform = srsFactory.GetTransform(srsDs, srsMap);
}
else
{
srsXform = null;
}
String featureClassName = selLayer.GetFeatureClassName();
String filter = sel.GenerateFilter(selLayer, featureClassName);
if (filter == null || filter.Length == 0)
{
continue;
}
MgFeatureQueryOptions query = new MgFeatureQueryOptions();
query.SetFilter(filter);
MgResourceIdentifier featureSource = new MgResourceIdentifier(selLayer.GetFeatureSourceId());
MgFeatureReader features = _featSvc.SelectFeatures(featureSource, featureClassName, query);
if (features.ReadNext())
{
MgClassDefinition classDef = features.GetClassDefinition();
String geomPropName = classDef.GetDefaultGeometryPropertyName();
do
{
MgByteReader geomReader = features.GetGeometry(geomPropName);
MgGeometry geom = agfRW.Read(geomReader);
if (!chkMergeBuffers.Checked)
{
geomBuffer = geom.Buffer(dist, measure);
if (geomBuffer != null)
{
if (srsXform != null)
{
geomBuffer = (MgGeometry)geomBuffer.Transform(srsXform);
}
Util.AddFeatureToCollection(propCollection, agfRW, featId++, geomBuffer);
bufferFeatures++;
}
}
else
{
if (srsXform != null)
{
geom = (MgGeometry)geom.Transform(srsXform);
}
inputGeometries.Add(geom);
}
}while (features.ReadNext());
features.Close();
}
}
if (chkMergeBuffers.Checked)
{
if (inputGeometries.GetCount() > 0)
{
double dist = srsMap.ConvertMetersToCoordinateSystemUnits(distance);
MgCoordinateSystemMeasure measure;
if (!arbitraryMapSrs)
{
measure = srsMap.GetMeasure();
}
else
{
measure = null;
}
MgGeometryFactory geomFactory = new MgGeometryFactory();
geomBuffer = geomFactory.CreateMultiGeometry(inputGeometries).Buffer(dist, measure);
if (geomBuffer != null)
{
Util.AddFeatureToCollection(propCollection, agfRW, featId, geomBuffer);
bufferFeatures = 1;
}
}
}
if (propCollection.GetCount() > 0)
{
commands.Add(new MgInsertFeatures("BufferSchema:Buffer", propCollection)); //NOXLATE
//Insert the features in the temporary data source
//
Util.ReleaseReader(_featSvc.UpdateFeatures(fsId, commands, false), commands);
}
// Save the new map state
//
layer.ForceRefresh();
_viewer.RefreshMap();
//build report message
if (newBuffer)
{
MessageBox.Show(string.Format(Strings.MsgBufferLayerCreated, txtBufferLayer.Text));
}
else
{
MessageBox.Show(string.Format(Strings.MsgBufferLayerUpdated, txtBufferLayer.Text));
}
}
finally
{
btnCreate.Enabled = true;
}
}
