Infrastructure Modeling DevBlog

Creating an AutoCAD Selection set from the Map 3D Platform Selection set (MgSelectionBase)

By Partha Sarkar

I am trying to use the _MAPExtractFeatureGeometry from SendStringToExecute() function and I want to pass in a selection set containing a Feature.

I have used Geospatial Platform API and MgSelectionBase to create a selection set of features which works fine as I could highlight the selected feature. But subsequently when I use  _MAPExtractFeatureGeometry from SendStringToExecute(), I get a response 'there is no selection set'. I am trying to avoid the user intervention to select a feature on screen. How do I create an AutoCAD selection set from the Map 3D Platform Selection set (MgSelectionBase) ?

 

Here is a C# .NET code snippet demonstrating the same –

// Select using API

 

AcMapMap currentMap = AcMapMap.GetCurrentMap();

MgLayerCollection layers = currentMap.GetLayers();

MgLayerBase layer = layers.GetItem("Parcels");

MgFeatureService fs = AcMapServiceFactory.GetService(MgServiceType.FeatureService) as MgFeatureService;

 

string fsId = layer.GetFeatureSourceId();

string className = layer.GetFeatureClassName();

MgFeatureQueryOptions query = new MgFeatureQueryOptions();

// Parcels.shp file is used

query.SetFilter("FeatId = 166");

MgResourceIdentifier resId = new MgResourceIdentifier(fsId);

MgFeatureReader featureReader = fs.SelectFeatures(resId, className, query);

 

MgSelectionBase selectionSet = new MgSelectionBase(currentMap);

 

//Add all features in the feature reader to the newly constructed selection set

selectionSet.AddFeatures(layer, featureReader, 0);

 

AcMapFeatureEntityService.HighlightFeatures(selectionSet);

 

// create an AutoCAD selection

SelectionSet acadSelSet = AcMapFeatureEntityService.AddFeaturesToSelectionSet(null, selectionSet);

 

// force this selection to be the implied selection

Autodesk.AutoCAD.ApplicationServices.Application.DocumentManager.MdiActiveDocument

.Editor.SetImpliedSelection(acadSelSet);  

Autodesk.AutoCAD.ApplicationServices.Application.DocumentManager.MdiActiveDocument

.SendStringToExecute("_MAPExtractFeatureGeometry ", true, false, true);

 

featureReader.Close();        

 

Hope this is useful to you!

Read FDO Geometry Area using Map 3D Geospatial Platform API

By Partha Sarkar

This blog post is in continuation of my previous post on “Read FDO Features using Map 3D Geospatial Platform API “.  If FDO data has an attribute like “AREA” as shown in the screenshot below, we can access the same using –

String ftrArea = ftrRdr.GetString("AREA");

However, if there is no such attribute or data field in the FDO data, in that case we can use the following approach to access the FDO Geometry and then read it’s Area property as shown in the code snippet below :

//this is just an example; change this according to your FDO data field                         

 

Int32 ftrID = ftrRdr.GetInt32("Autogenerated_SDF_ID");

 

String geomName = mapLayer.GetFeatureGeometryName();

MgByteReader reader = ftrRdr.GetGeometry(geomName);

 

// Convert the AGF binary data to MgGeometry.

MgAgfReaderWriter agfReaderWriter = new MgAgfReaderWriter();

MgGeometry ftrGeometry = agfReaderWriter.Read(reader);

ed.WriteMessage("\n Feature ID :" + ftrID.ToString() + "   " + "Area  : " + ftrGeometry.Area.ToString());  

 

Hope this is useful to you! 

Coordinate Transformation using Map LISP API

By Partha Sarkar

How to use Map LISP function ade_projptforward to achieve coordinate transformation?

Here is a LISP code snippet which demonstrates the same :

;;Coordinate transform with LISP.

(defun c:test (/ result pt)
(ade_projsetsrc "LL")
(ade_projsetdest "ITALY-W-ROME")
(setq result(ade_projptforward '(10.0 20.0 0.0)))
(if(null result)
(prompt "\nError in Transformation ")
)
(princ result)
);

How to apply Coordinate system Transformation using API ?

By Partha Sarkar

Coordinate transformation is converting a coordinate value in one coordinate system into its equivalent in another coordinate system.

 

Here is a code snippet which demonstrates how to apply Coordinate system Transformation :

string mapSRS = AcMapMap.GetCurrentMap().GetMapSRS();

// create MgCoordinateSystemFactory

 

MgCoordinateSystemFactory coordSysFactory = new MgCoordinateSystemFactory();

MgCoordinateSystem coordSys = coordSysFactory.Create(mapSRS);

MgCoordinateSystem wgs84Sys = coordSysFactory.Create("LL84");

 

MgCoordinateSystemTransform transform = null;

transform = coordSysFactory.GetTransform(wgs84Sys, coordSys);

MgCoordinate coord = null;

 

if (transform != null && !double.IsNaN(pnt.Z))

  coord = transform.Transform(pnt.Z, pnt.Y, pnt.Z);

else

  coord = transform.Transform(pnt.X, pnt.Y);   

How to extract all of the coordinate system names from the coordinate system dictionary and check their types?

By Partha Sarkar

The OSGeo.MapGuide namespace contains a set of classes which provide access to the CS-Map coordinate system functionality in both the AutoCAD Map 3D and AIMS environments.

MgCoordinateSystemFactory is used to get the coordinate system categories and the coordinate system catalog. MgCoordinateSystemCatalog is used to get the various dictionaries. MgCoordinateSystemDictionary contains definitions for earth-based and Euclidean coordinate systems.

 

The following code extracts all of the coordinate system names from the coordinate system dictionary, uses those names to get all of the coordinate system definitions from the dictionary and determines the type (projected, geographic, or arbitrary) of the coordinate system from its definition, and stores the name of the coordinate system in a list according to its type. 

//MgCoordinateSystemFactory coordSysFactory = new MgCoordinateSystemFactory();

//MgCoordinateSystemCatalog csCatalog = coordSysFactory.GetCatalog();

//MgCoordinateSystemDictionary csDict = csCatalog.GetCoordinateSystemDictionary();

 

List<string> arbitraryCs = new List<string>();

List<string> geographicCs = new List<string>();

List<string> projectedCs = new List<string>();

 

MgCoordinateSystemEnum csDictEnum = csDict.GetEnum();

 

// the following line gets all of the names in the dictionary

// uses those names to get the coordinate

 

int csCount = csDict.GetSize();

ed.WriteMessage("\nCoordinate System Count : " + csCount.ToString());

ed.WriteMessage("\n-------------------------------------------------");

MgStringCollection csNames = csDictEnum.NextName(csCount);

string csName = null;

MgCoordinateSystem cs = null;

int csType = 0;

for (int i = 0; i < csCount; i++)

{

 csName = csNames.GetItem(i);

 cs = csDict.GetCoordinateSystem(csName);

 csType = cs.GetType();

 if (csType == MgCoordinateSystemType.Arbitrary)

 {

  arbitraryCs.Add(csName);

 }

 else if (csType == MgCoordinateSystemType.Geographic)

 {

  geographicCs.Add(csName);

 }

 else if (csType == MgCoordinateSystemType.Projected)

 {

  projectedCs.Add(csName);

 }         

  ed.WriteMessage("\nCoordinate System Name : " + csName.ToString());

}

Hope this is useful to you !

How do I turn off the "Map Base" Layer using API ?

By Partha Sarkar

I want to switch off the "Map Base" Layer using Map 3D API.

 

 

Here is the C# code snippet to switch off the "Map Base" Layer using Map 3D managed API :

 

using (Transaction trans = transManager.StartTransaction())

{

  ObjectId managerId = DisplayManager.Create(proj).MapManagerId(proj, true);

  MapManager manager = trans.GetObject(managerId, OpenMode.ForRead) as MapManager;

 

  if (manager != null)

  {

    currentmapId = manager.CurrentMapId;

    Map currentMap = trans.GetObject(currentmapId, OpenMode.ForWrite) as Map;

 

    IEnumerator iterator = currentMap.NewIterator(true, true);

    Type groupType = typeof(Autodesk.Gis.Map.DisplayManagement.Group);

 

    while (iterator.MoveNext())

    {

      ObjectId objId = (ObjectId)iterator.Current;

      Object obj = trans.GetObject(objId, OpenMode.ForWrite);             

 

      if (obj.GetType().ToString() == "Autodesk.Gis.Map.DisplayManagement.BaseElement")

      {

        Autodesk.Gis.Map.DisplayManagement.BaseElement baseElement = obj as Autodesk.Gis.Map.DisplayManagement.BaseElement;

        baseElement.SetVisible(0, false);                 

      }

    }

  }

  trans.Commit();

}

Setting the center of map programmatically

By Daniel Du

Q:

I want to zoom a map to screen center without changing the scale, how can I do it?

A:

In Map 3D, if you have connected to FDO data sources and added them to map, you can call AcMapMap.ZoomToExtent() to set the new extent of map. It accepts a parameter of MgEnvelope type. If you do not have any FDO data source, you have to use AutoCAD APIs.

Here is the code snippet:

[CommandMethod("ZoomCenter")]

public void ZoomCenter()

{

 

  Document doc = Application.DocumentManager.MdiActiveDocument;

  Database db = doc.Database;

  Editor ed = doc.Editor;

 

  Point3d centerPt;

  PromptPointOptions ppo = new PromptPointOptions(

                "Click on map to zoom center:");

  PromptPointResult ppr = ed.GetPoint(ppo);

  if (ppr.Status == PromptStatus.OK)

  {

    centerPt = ppr.Value;

 

    AcMapMap map = AcMapMap.GetCurrentMap();

    MgEnvelope mapExtent = map.GetMapExtent();

 

    //There are some FDO feature source connected

    if (!mapExtent.IsNull())

    {

      double centerX = mapExtent.LowerLeftCoordinate.X

                                    + mapExtent.Width / 2;

      double centerY = mapExtent.LowerLeftCoordinate.Y

                                    + mapExtent.Height / 2;

 

      ed.WriteMessage("center:" + centerX.ToString() + ","

                                + centerY.ToString() + "\n");

 

      MgEnvelope newExtent = new MgEnvelope(

          centerPt.X - mapExtent.Width  / 2,

          centerPt.Y - mapExtent.Height / 2,

          centerPt.X + mapExtent.Width  / 2,

          centerPt.Y + mapExtent.Height / 2);

 

      map.ZoomToExtent(newExtent);

    }

    else

    {

      //no FDO data source, use AutoCAD API

 

      using (Transaction Tx = db.TransactionManager.StartTransaction())

      {

        ed.UpdateTiledViewportsInDatabase();

 

        ViewportTableRecord viewportTableRec = Tx.GetObject(

                ed.ActiveViewportId, OpenMode.ForWrite)

                as ViewportTableRecord;

 

        viewportTableRec.CenterPoint = new Point2d(centerPt.X, centerPt.Y);

 

        ed.UpdateTiledViewportsFromDatabase();

 

        Tx.Commit();

      }

 

    }

 

 

  }

 

}

Hope this helps.

Computing distance in Map 3D using .net API

By Daniel Du

Q:

I'd like to know how to compute the distance of two points using map 3D .net API

A:

Here is the sample code to do the same:

[CommandMethod("ComputeDistance")]

public void ComputeDistance()

{

  //I am using LL84 coordinate system for this sample

  //As an example, this is longtitue/Latitue of two points

 

  double x1 = -87.7104750022991;

  double y1 = 43.7017449116101;

  double x2 = -87.703061972587;

  double y2 = 43.7016702994388;

 

  Editor ed = Application.DocumentManager.MdiActiveDocument.Editor;

 

  //Get coordinate system of current map

  AcMapMap currentMap = AcMapMap.GetCurrentMap();

  string srsWkt = currentMap.GetMapSRS();

  //ed.WriteMessage("srs = " + srsWkt + "\n");

  MgCoordinateSystemFactory coordSysFactory = new MgCoordinateSystemFactory();

  MgCoordinateSystem coordSys = coordSysFactory.Create(srsWkt);

 

  //compute gread circle distance

  double distance = coordSys.MeasureGreatCircleDistance(x1, y1, x2, y2);

  distance = coordSys.ConvertCoordinateSystemUnitsToMeters(distance);

  ed.WriteMessage("gread circle dist = " + distance.ToString() + "\n");

 

  //compute Euclidean distance

  distance = coordSys.MeasureEuclideanDistance(x1, y1, x2, y2);

  distance = coordSys.ConvertCoordinateSystemUnitsToMeters(distance);

  ed.WriteMessage("Euclidean distance = " + distance.ToString() + "\n");

 

 

  //Another method, compute the distance from Newyork to Boston

  MgCoordinateSystemMeasure coordSysMeasure = coordSys.GetMeasure();

  double dist = coordSysMeasure.GetDistance(-74.806394, 40.714169, -71.061342, 42.355892);

  dist = coordSys.ConvertCoordinateSystemUnitsToMeters(dist);

  ed.WriteMessage("distance from Newyork to Boston is " + dist.ToString() + " meters\n");

 

}

Hope this helps.

Mimic OnCheckIn Event in Map 3D

By Daniel Du

Have you ever run into such scenario that value of one field is computed from other fields? For example we can calculate the total land price from land price per unit and area. Most RDBMS can do this with stored procedure or similar, but for file based data source like SDF or SHP file, we have to do this ourselves by code. A good chance to do the calculation is when the modification on features are checked in, it would be better if Map 3D API provides an event like OnCheckIn, unfortunately, there is no such events in Geospatial Platform API.

But it can be done by a combination of serials events, they are FeatureInserted, FeatureUpdated, and FeatureDeleted of AcMapFeatureService.

Following code snippet domos how to use FeatureUpdated event, you can implement the other two to mimic OnCheckIn event is Map 3D. I blogged this topic a few years ago on my personal blog, if you prefer to read Chinese, please go to this link.

// (C) Copyright 2012 by Autodesk

//

using System;

using Autodesk.AutoCAD.Runtime;

using Autodesk.AutoCAD.ApplicationServices;

using Autodesk.AutoCAD.DatabaseServices;

using Autodesk.AutoCAD.Geometry;

using Autodesk.AutoCAD.EditorInput;

using Autodesk.Gis.Map.Platform;

using OSGeo.MapGuide;

 

// This line is not mandatory, but improves loading performances

[assembly: CommandClass(typeof(MapOnCheckIn.MyCommands))]

 

namespace MapOnCheckIn

{

 

  // This class is instantiated by AutoCAD for each document when

  // a command is called by the user the first time in the context

  // of a given document. In other words, non static data in this class

  // is implicitly per-document!

  public class MyCommands

  {

 

    Editor ed = Application.DocumentManager.MdiActiveDocument.Editor;

 

    //Following code sample demos how to mimic the OnCheckIn event

    //when features are updated using fs.FeatureUpdated,

    //You will need to implement FeatureDeleted  and FeatureInserted yourself

 

    [CommandMethod("AddFsEvents")]

    public void AddFeatureServicesEvents()

    {

      AcMapFeatureService fs = AcMapServiceFactory.GetService(

          MgServiceType.FeatureService) as AcMapFeatureService;

      fs.FeatureUpdated += new FeatureUpdatedHandler(fs_FeatureUpdated);

 

      AcMapMap currentMap = AcMapMap.GetCurrentMap();

      currentMap.FeatureInstanceModified +=

        new FeatureInstanceModifiedHandler(currentMap_FeatureInstanceModified);

 

    }

 

    // Triggered at once when the feature is updated, even it is not checked,

    // This is not what we want.

    void currentMap_FeatureInstanceModified(object sender,

                  FeatureInstanceModifiedEventArgs args)

    {

      ed.WriteMessage("FeatureInstanceModified fired!\n");

    }

 

    //Triggered  when the feature is updated and checked in,

    //This is the one we are looking for, it only triggered when checked in

    //We can get the feature which is being checked in by Args.GetFeature()

    //if more than one feature are updated, it will be fired more times

    //so you can catch this even to handle each of them.

 

    void fs_FeatureUpdated(object sender, AcMapFeatureEventArgs args)

    {

      ed.WriteMessage("FeatureUpdated fired! \n");

      ed.WriteMessage("--------------------------------------\n");

      ed.WriteMessage("PropertyName \t Value \n");

      MgPropertyCollection props = args.GetFeature().GetProperties();

      int propCount = props.Count;

      for (int i = 0; i < propCount; i++)

      {

        string propName = props[i].Name;

        string valueString;

        switch (props[i].GetPropertyType())

        {

          case MgPropertyType.Int32:

            valueString = (props[i] as MgInt32Property).GetValue().ToString();

            break;

 

          case MgPropertyType.String:

            valueString = (props[i] as MgStringProperty).GetValue().ToString();

            break;

          default:

            valueString = "<***>";

            break;

        }

 

        ed.WriteMessage(propName + "\t" + valueString + "\n");

 

      }

      ed.WriteMessage("------------------------------------\n");

 

 

    }

  }

 

}

 

Hope this helps.

Run MapClean Profile Programmatically

By Daniel Du

I am asked about this question: I know how to run MapClean in Map 3D UI, and create/load profile to run the mapclean process repeatedly. but is it possible to run MapClean programmatically? The answer is yes, we can create a custom application to run MapClean programmatically. Here is the code snippet, please note that we need add reference to ManagedMapApi.dll, which can be found in Map 3D installation folder.

 

[CommandMethod("clean")]

public void DrawingCleanUp()

{

 

  Autodesk.Gis.Map.Topology.Variable cadAction

    = new Autodesk.Gis.Map.Topology.Variable();

  //Load a dpf script file, it can be generated by Map 3D UI

  cadAction.LoadProfile(@"C:\Work\mapcleandemo\myclean.dpf");

 

  //The instance of this class should be released by explicitly

  //calling Dispose() in order to avoid memory leak.

  //See also 'using' keyword in C# and VB.Net.

  using (Autodesk.Gis.Map.Topology.TopologyClean cadCleanobj

      = new Autodesk.Gis.Map.Topology.TopologyClean())

  {

    //Open a file first to avoid MapTopologyException

    cadCleanobj.Init(cadAction, null);

 

    cadCleanobj.Start();

    cadCleanobj.GroupNext();

 

    while (!cadCleanobj.Completed)

    {

      cadCleanobj.GroupFix();

      cadCleanobj.GroupNext();

    }

    //Commit the changes

    cadCleanobj.End();

  }

 

}

 

Hope this helps!