Latest Contributions by MelitaKennedy

Melita,     Thank you for your response. Yes, you are correct that the original GEOTIFF is georeferenced, but comes without a spatial reference. Based on the Burned Area Manual, I know that it has been reprojected from Sinusoidal into Platte Carree and I have no issue defining the GEOTIFF as such. But, I need a projection that preserves accurate area, thus Albers Equal Area for the CONUS is ideal. this is what I've tried: 1) in ArcCatalog or in ArcMap (I've tried both) I go to the properties of the file and give it the Platte Carre projection (WGS 1984 by default) I've also now tried assigning it your MODIS sphere proj. 2) Next, I go to arctools > data management > projections and transformations > project raster. For the input I select my file, it shows up with the projection I just defined. For the Output I select projected coord system > continental > North America > USA contiguous Albers Equal Area Conic.prj. Then the only thing I modify is the latitude of origin, changing it to 0.here it is: Projection: Albers False_Easting: 0.000000 False_Northing: 0.000000 Central_Meridian: -96.000000 Standard_Parallel_1: 29.500000 Standard_Parallel_2: 45.500000 Latitude_Of_Origin: 0.000000 Linear Unit: Meter (1.000000) Geographic Coordinate System: GCS_North_American_1983 Angular Unit: Degree (0.017453292519943299) Prime Meridian: Greenwich (0.000000000000000000) Datum: D_North_American_1983   Spheroid: GRS_1980     Semimajor Axis: 6378137.000000000000000000     Semiminor Axis: 6356752.314140356100000000     Inverse Flattening: 298.257222101000020000 3) as soon as I click ok, the program shuts down, every time. It seems the issue is related to the difference in spheroids between WGS 1984 or MODIS sphere and the NAD_1983 sphere. And I know there are some transformations for them, but I can't get them to work. The other part of my problem is that I'm trying to get this burned area to overlay with a cropland data layer I have for the conUS. here is the projection specs for my cropland data layer. This is not the original projection, it was reprojected to this via ArcInfo workstation, which as I understand assigns the clarke 1866 as a default. I was given this layer this way and have to make the two overlay accurately to determine how much burned area is cropland. Projection: Albers False_Easting: 0.000000 False_Northing: 0.000000 Central_Meridian: -96.000000 Standard_Parallel_1: 29.500000 Standard_Parallel_2: 49.500000 Latitude_Of_Origin: 0.000000 Linear Unit: Meter (1.000000) Geographic Coordinate System: GCS_Clarke_1866 Angular Unit: Degree (0.017453292519943299) Prime Meridian: Greenwich (0.000000000000000000) Datum: D_Clarke_1866   Spheroid: Clarke_1866     Semimajor Axis: 6378206.400000000400000000     Semiminor Axis: 6356583.799998980900000000     Inverse Flattening: 294.978698200000000000 So, I've also been trying to get the Burned Area geotiffs into a Clarke_1866 albers, but I have no experience with ArcInfo Workstation and haven't been able to reproduce what my boss did with the cropland data layer. This must seem like such a mess, at least it does to me, hopefully not to you! The bottom line is that I need to get those GEOTIFF's into an alber's projection that will overlay with the cropland data. Is this possible? Can you think of any way to do so? Again, thank you, Melita. Carie Hi Carie, Oh yes, transformations. First, can you overlay the Burned Area data with your cropland data in ArcMap? Ignore any warning messages about different geographic coordinate systems/datums. If it doesn't overlay, try defining a custom transformation. Data frame properties, coordinate systems tab, Transformations, New button. In the New dialog, change the method to Geocentric Translation, leave the parameters set to zeroes. OK all and recheck the overlay. Melita ... View more

Hi, I have recently been compiling maps of Iraq on a UTM grid scale, so have been using the UTM grid projection system. However, what projection system and datum should one use when creating maps for the entire country? EPSG (http://www.epsg-registry.org) lists an ED 1950-based Iraq National Grid. We added it at version 10. It uses transverse Mercator. There's a newer IGRS geographic coordinate system, and 3 UTM zones based upon it, but no national grid. IGRS is close to WGS84, according to EPSG. Melita ... View more

Hello, I can't seem to get anywhere on this and am hoping that Melita might be able to respond since it seems she has been dealing with this a lot, but any help is appreciated. I have GEOTIFF's provided by the University of Maryland for MODIS burned area for years 2003-2009. I am trying to take it from a Platte Carree projection to Albers Equal Area. please help me! thanks to all, Carie Hi Carie, What problem are you encountering? Does the original, downloaded GeoTIFF have a defined coordinate system? Does the data appear to be georeferenced (the raster extents aren't equal to 'cell' counts)? If georeferenced but doesn't have a defined coordinate system, what happens when you define its coordinate system? The online documentation is less than helpful about the coordinate system details, isn't it? I'm assuming that the GCS/Datum is the same as listed in the report. Because I can't attached .prj files, here's the contents of the a .prj file: GEOGCS["GCS_MODIS_Burned_Area_Sphere",DATUM["D_MODIS_Burned_Area_Sphere",SPHEROID["MODIS_Burned_Area_Sphere",6371007.181,0.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] Place it in a text file (no word wrapping, single line) with a .prj extension and place it in coordinate systems, geographic coordinate systems. Now you can create a custom projected coordinate system with this as its geographic coordinate system. Just leave the Plate Carree parameters as-is. I assume the central meridian is zero. Melita ... View more

Again - thanks for your reply - I have the standard toolbar active and the scale is not on there as defaulted by the add data button. Apparently some predecessor of mine turned it off or something. 😞 I have looked in the customized options to add that 'button' (which I have not found) as well as going through the options on the data frame properties. This should be fairly simple can anyone please help - Thanks!!!! Hi, I'm not sure this will work, but click the arrow at the right of the standard toolbar and select "Customize". If you hover over it, you should see "Toolbar options". This will open the customize dialog. Select the Toolbars tab and highlight the Standard toolbar. Click the Reset button at right. This will get rid of any custom buttons on the toolbar--it removed my project wizard button, for instance. Hopefully, it will replace the missing control! I also saw a "Map Scale" control on the Commands tab, Pan/Zoom set, but wasn't able to drop it onto a toolbar, so I don't know if it's the same control or not. Melita ... View more

Argh.... Its things like this that make me despise working with ArcSDE and ESRI products.  How is there two different WGS 84s?  SRID 4326 (and Oracle's coordsys) are geospatial standard WGS 84 definitions. EPSG:4326 .prj definition: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]] ArcGIS WGS 84 .prj definition: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]] These are identical definitions.  How does ArcSDE get confused?  How does Oracle Spatial screw up 4326? Here's Oracle's 4326 from 11g: GEOGCS [ "WGS 84", DATUM ["World Geodetic System 1984 (EPSG ID 6326)", SPHEROID ["WGS 84 (EPSG ID 7030)", 6378137, 298.257223563]], PRIMEM [ "Greenwich", 0.000000 ], UNIT ["Decimal Degree", 0.01745329251994328]] Here's ours from ArcGIS 10: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984", 6378137.0, 298.257223563]],PRIMEM["Greenwich", 0.0],UNIT["Degree", 0.0174532925199433],AUTHORITY["EPSG",4326]] We've started including an AUTHORITY tag for the top level object. This would be one of the easier entries to try to "make equal", but it's complicated in code to do it reliably without just taking the SRIDs at face value. The values are all 'equal', but does that really mean that the WKT versions (if you didn't have the SRIDs) represent the same coordinate reference system? Melita ... View more

Hi there, I want to reproject a few raster datasets from GCS_WGS_1984 to NAD_1983.  The data is for the province of British Columbia in Canada.  Could someone help me choose the correct transformation?  My options range for NAD_1983_To_WGS_1984_(1 through 8).  There is no mention of British Columbia in this ESRI document (http://resources.arcgis.com/content/kbase?fa=articleShow&d=24159) describing each of these transformations. Thanks, Rob Hi Rob, What's the cell size on the rasters? If they're more than 1m, don't worry about it. If they're significantly less, try using NAD_1983_To_WGS_1984_5. Melita ... View more

Hi Melita, Thank you for your reply. The projection information is correct (NAD83 Albers USGS version). I created a new raster in UTM 18N using the Project Raster tool. Now the data aren't aligning (screen shot attached) with another feature with the same projection information (UTM 18N). The planar distance units for the raster and feature are both meters. I've had this problem in the past, but I can't remember what the issue was. [ATTACH]3120[/ATTACH] Why aren't the data lining up? Please let me know if you need more info. Again, thank you. Regards, Mike Hi Mike, Does the original raster line up? If it doesn't, then we have an issue with the original raster. If it does, then we can focus on what happened when the raster was projected. How far off is it? Are you certain that the other data is in the correct location? Melita ... View more

Hi there...                I am a student jus started on my college project.. I am thinking of finding the best projection suited for various regions of the  world taking certain factors into consideration....  well I am not clear enough like what image format should I prefer to use in my scenario here?!! where can I possibly get the  free data to implement this... i am jus confused like where/ how to start this... So any ideas or suggestions are truly welcome... Thank you... Hi, I can't help you on where to find a free world raster. See if you can get access to the Esri Data & Maps CDs/DVDs. Or check USGS website. For your project, see if you can find a copy of Canters, Frank. Small Scale Projection Design. Also useful might be: Snyder, John P. Flattening the Earth. which discusses the different projections created over the centuries and what they were used for. Melita ... View more

Hi, I attempted to find a thread that already exists for this problem, but I could not find one. I would really appreciate some guidance with this. I am having trouble projecting a raster grid. When I try to load the file, I get a message saying "error, inconsistent extent.." Also, when I try to transform the projection it fails with a general 999999 message, with no useful information. I defined the projection as: North American Datum of 1983. Then attempted to transform the projection to: NAD_1983_UTM_Zone_18N. The process fails and I get the inconsistent extent message. Extent information from the layer property source table is as follows: Top: 2584425 Left: 1605495 Right: 1932525 Bottom: 2123805 Here is the spatial ref information from the metadata: Spatial_Reference_Information: Horizontal_Coordinate_System_Definition: Planar: Map_Projection: Map_Projection_Name: Albers Conical Equal Area Albers_Conical_Equal_Area: Standard_Parallel: 29.500000 Standard_Parallel: 45.500000 Longitude_of_Central_Meridian: -96.000000 Latitude_of_Projection_Origin: 23.000000 False_Easting: 0.000000 False_Northing: 0.000000 Planar_Coordinate_Information: Planar_Coordinate_Encoding_Method: row and column Coordinate_Representation: Abscissa_Resolution: 30.000000 Ordinate_Resolution: 30.000000 Planar_Distance_Units: meters Geodetic_Model: Horizontal_Datum_Name: North American Datum of 1983 Ellipsoid_Name: Geodetic Reference System 80 Semi-major_Axis: 6378137.000000 Denominator_of_Flattening_Ratio: 298.257222 What is going wrong here? Hi, The metadata description looks like NAD83 Albers USGS version (in projected coordinate systems, continental, north america) except the latitude of origin parameter is missing. Based on the raster's extent, I think it is the USGS Albers because the extents place the raster around northern NJ, eastern NY up to southern VT. If the raster does not represent that part of the world, the Albers coordinate system information could be incorrect. I'm confused by your statement that you defined the coordinate system as NAD83--the geographic coordinate system? You need to leave it as Albers (hopefully, it's correct) and then use the Project Raster to create a new raster in UTM 18N. The 'inconsistent extent' is classic if a dataset is georeferenced to a projected coordinate system, but is assigned a geographic coordinate system. "Projecting" data via the Define Projection tool or the data's property page in ArcCatalog updates the metadata only. The data values are still relative to the previous coordinate system. Melita ... View more

I have similar problem when I tried to migrate my 9.3 C# code to 10.0. I got error message for ISpatialReferenceFactory spatialReferenceFactory2 = new SpatialReferenceEnvironmentClass(); the error message is "Retrieving the COM class factory for component with CLSID {7B5B7020-C4F6-11D1-BC92-0000F875BCCE} failed due to the following error: 80040111." Then I tried to use activator so I revised the code to following: // Set up the SpatialReferenceEnvironment.             // SpatialReferenceEnvironment is a singleton object and needs to use the Activator class.             Type factoryType = Type.GetTypeFromProgID("esriGeometry.SpatialReferenceEnvironment");             System.Object object1 = Activator.CreateInstance(factoryType);             ISpatialReferenceFactory2 spatialReferenceFactory2 = object1 as ISpatialReferenceFactory2; However, I still got the same error message. Can anyone help on this one? Thanks a lot. Wei Hi Wei, Another post with a similar issue said that he needed to add this: ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Engine); (or other licensing information). Melita ... View more

Hi Melita, Your understanding of the problem sounds correct but I may not be attempting to solve it properly.  If you have any samples source code that I can look at, I would appreciate it. Just to clarify what i m attempting:  I am reading data from our DB, creating features in several feature classes that exist in in-memory-workspaces, After building all of the feature classes, I want to center the map on the visible features. The original data is stored in decimal degrees but the map is displayed using Mercator projection with meter units.  I am able to project the original data to mercator without any issue. After all the data has been processed, I want to cener the map on the visible data.  I am attempting to do this by changing the central meridian from 0, to an appropriate value.  This seems to be where my problem is.  I am calculating the visible data by examing the extent of certain feature classes of my data. While running the application, I can reload some data and I want to repopulate some (but not all) of the feature classes and re-center on this data.  For performance reasons, I am attempting to leave some of the feature classes alone and not re-load / re-process this data. On the map, I am displaying the sample world shape file that came with the sample application data.  When I recenter the map, the shape file displays correctly, with the new central meridian.  My feature class data does not appear correct.  For example, if I change the central meridian to 180 degrees, the IDL of the shape file is displayed on the center of my map.  The coordinates that I display tracking the mouse movement are correct as well.  It is just the features related to the feature classes I created, that do not display correctly. Any ideas on what I am doing wrong? Thanks, Gary Hi Gary, The only thing I can think of is that the new feature layers don't have any coordinate system attached to them. If you check them in ArcMap, is their coordinate system "Unknown"? Melita ... View more

I created an extension to ArcMap (9.3.1) that employs a MapControl.  I recently upgraded to ArcGIS Desktop 10 and found to my dismay that the MapControl - which resides in ESRI.ArcGIS.Controls library - now requires an Engine license.  Setting aside my deep displeasure at this change in marketing strategy, my operative question is - how much do I have to spend to enable the MapControl so that I can upgrade my extension to run in ArcMap? I have had several conversations with ESRI's sales & marketing people and I'm still confused.  It appears that if I purchase an EDN subscription I will receive the ArcGIS Engine SDK and I should be able to convert my 9.3.1 extension to 10.0.   My problem is that the ESRI folks keep talking about "production level" licenses - assuming, I guess, that they think I am trying to create a stand alone app that will be available over the web.  But, my program is an extension to ArcMap and will only be used within that environment - therefore the user has a ArcGIS Desktop license and should not need any additional licensing run to run my extension with its embedded MapControl? I would appreciate hearing from anyone who has had to convert a MapControl bearing tool, command or extension from 9.3.x to 10 and what it took license wise to make it work. Thanks. I work at Esri in development; I freely admit that I do not understand all the licensing details. I did find this statement on the "How to create a MapControl application" page (http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#//000100000008000000😞 Running the application with an ArcGIS Desktop license Since the application is set by default, it requires an ArcGIS Engine Runtime license. However, if you want it to run on a machine with ArcGIS Desktop (ArcView, ArcEditor, or ArcInfo), or to check for and consume a floating ArcGIS Desktop license, update the LicenseControl to verify the other licenses as well. Do the following steps to run the application with a license: Right-click the LicenseControl and choose Properties. The Properties dialog box appears. Select the check boxes of the licenses you want to use when the application runs. The application always uses the lowest level of licensing that is available. For example, if a machine has an ArcGIS Engine and ArcEditor license, the application uses the ArcGIS Engine license. If there are ArcView and ArcInfo floating licenses available and both of those options were selected in the LicenseControl, the application uses the ArcView license. ----- When it says 'set by default,' it means that the VisualStudio Esri wizard automatically sets up licensing for Engine. Melita ... View more

I am working on an application where I want to change central meridian so that if the user is looking at data near the IDL, the data is in the center of the map and there isn't an issue where data must wrap around. I am adding the data to the feature class using standard central meridian = 0.  I then get the dataset extent's and then setting the map's spatialreference central merian to the midpoint of the extent. I have tried a variety of simple things such as setting the same modified spatial reference to the map layers.  I have removed and re-added the layers to the map. Is there a simple solution for this or do I need to retrieve each piece of geometry of each feature and somehow apply a transformation for the new central meridian? Also looking at the documentation for ILayer::putref_SpatialReference the Remarks say: This property is only used for map display, setting this property does not change the spatial reference of the layer's underlying data.  The ArcGIS framework uses this propert to pass the spatial reference from the map to the layer in order to support on-the-fly projection. To me, this sounds like exactly what I want.  I want to add the data to the feature class, change the map's spatial reference and see the data projected correctly. What am I missing? I load and display the world shape file (from the sample data) and this projects exactly like I would expect when changing the central meridian. I have tested this on 9.2 SP5 and 9.3.1 SP 2 on Solaris. Thanks, Gary Hi Gary, I can't understand what you're trying to do. It sounds to me like you're trying to change the coordinate system (spatial reference) of a layer/dataset without projecting it. I think you're on the right track that you just want to update the data frame's coordinate system. You want to set the central meridian of IMap::SpatialReference to the center of the area of interest. What coordinate systems will generally be used for the map and for your data layers? Melita ... View more

Hi,       Anyone please, can tell regarding the map projection. I have images of Moon, which is originally generated through Geographic lat/long, WGS 84 datum. I am working with ArcGIS 9.2. I worked applying the conversion factor of Moon ie 30.333 km is equal to 1 degree, in case of my applications. I try to reproject the image in Moon spheroid, it is done, I apply the flattening factor of the major and minor axis, I was very happy that I did it, when I check the projection it is showing moon spheroid, but in my application work the scale is not coming with respect to the Moon. And when I check the datum, it is showing unknown. So please tell me, is it not possible to project the images who were originally processed in Geographic Lat/long, WGS 84 datum to lunar projection (datum).       I ll be very thankful to you. with regards, guneshwar gune.ishwar@gmail.com I think you are running into a known issue or two in ArcMap/ArcGIS. Map projections and datum transformations for vector data are handled by the projection engine library (pe.dll). Rasters use both pe.dll and another library. The other library may not be able to save a moon-based coordinate system in that particular raster format. What format is the data? (tiff, jpeg, GRID) If you're trying to measure features in ArcMap, I believe that it was hard-coding an earth ellipsoid when calculating the measurements. You may want to check out PIGWAD (http://webgis.wr.usgs.gov/) which has some tools and forums specifically about non-earth data. Melita ... View more

Thanks for all the comments, feedback and help. I removed 'Class' from SpatialReferenceEnvironment and receive another COM exception: Retrieving the COM class factory for component with CLSID {7B5B7020-C4F6-11D1-BC92-0000F875BCCE} failed due to the following error: 80040111 ClassFactory cannot supply requested class (Exception from HRESULT: 0x80040111 (CLASS_E_CLASSNOTAVAILABLE)). I'm not a COM developer, so this is a little foreign to me. A bit more on the map I am using.  I was handed a MXD file that has Virginia's map data (Danville area), and in order to track across the US, ESRI's North America Street Map data was also added.  The map as previously indicated should be projected in the Virginia State Plane.  I did change by PCS to esriSRProjCSType.esriSRProjCS_NAD1983SPCS_VASouthFT, but as I mentioned, I get the error mentioned above. If I can get past that error, that would be great.  Is there something I have to do to put a wrapper around Arc Objects COM object? Thanks for the help. Bob Hi Bob, I tell people that I'm barely a programmer (I do enough for testing, but it sure isn't 'elegant') so I can't tell if you're using VB6 or VB.Net code, or maybe something else. For VB6, try    Dim pSpatRefFact As SpatialReferenceEnvironment    Set pSpatRefFact = New SpatialReferenceEnvironment and you also need to use    Set pGCS = pSpatRefFact.CreateGeographicCoordinateSystem(enum) For VB.Net, try         Dim t As Type = Type.GetTypeFromProgID("esriGeometry.SpatialReferenceEnvironment")         Dim obj As System.Object = Activator.CreateInstance(t)         Dim pSpatRefFact As ESRI.ArcGIS.Geometry.ISpatialReferenceFactory3 = obj Melita ... View more