Downloads
View tile download instructions
Download from a spreadsheet of links (To download in bulk, save links in column C to a text file and use it as input to wget, WinWGet or Persepolis Download Manager)
Overview
In spring 2023, MassGIS acquired statewide, 15 cm resolution, 16-bit, 4-band (RGB-IR) digital orthophotos. Funding for this imagery was provided by MassDOT and the State 911 Department.
The data were delivered in late 2023 as 10,218 individual 1,500 m x 1,500 m tiles in a GeoTIFF format. The tile naming convention is based on the U.S. National Grid (USNG), taking the coordinates of the southwest corner of the tile. See the Index layer. To allow for easier distribution of this free imagery, MassGIS compressed the GeoTIFFS into the JPEG 2000 format, which retains the IR band. MassGIS also serves the data as a cached tile service.
The imagery’s horizontal accuracy specification was ASPRS Class 1. The data were developed based on a horizontal projection/datum of NAD 1983 (2011) UTM Zones 18N and 19N (meters).
This digital orthoimagery was created to provide easily accessible geospatial data which are readily available to enhance the capability of federal, state, and local emergency responders, as well as to plan for homeland security efforts.
These images can serve a variety of purposes, from general planning to field reference for spatial analysis, to a tool for data development and revision of vector maps. The imagery can also serve as a reference layer or basemap for myriad applications inside geographic information system (GIS) software and web-based maps.
Production
Aerial Acquisition
Planning, acquisition, processing, and creation of derivative products by MassGIS and NV5 Geospatial of Lexington, KY. The imagery’s horizontal accuracy specification was ASPRS Class 1.
Aerial imagery was acquired using an UltraCam Eagle M3 camera with a flight design that included a total of 214 flight lines to capture 11,158 photos. Aerial imagery was supplemented with the simultaneous acquisition of airborne GPS/IMU data, which captured the ground coordinate for the nadir point of each photograph.
Imagery was collected on 15 days, from 2/21 to 4/28, under conditions free from clouds and cloud shadows, smoke, fog, haze, light streaks, snow, ice on water bodies, flooding, excessive soil moisture, and foliage.
Within the UltraMap software suite used by NV5 Geospatial, raw acquired images were radiometrically and geometrically corrected using the camera’s calibration files and output as Level 2 images. The resulting radiometry was then manually edited to ensure each image has the appropriate tone and no pixels are clipped, and to blend each image with its neighbors.
Once radiometry had been edited, separate RGBI and Panchromatic images were blended together to form 10,218 single level 4 band TIFF images.
Photo position and orientation were calculated by linking the time of image capture, the corresponding aircraft position and attitude, and the smoothed best estimate of trajectory (SBET) data in POSPac software. Aerotriangulation was then performed using the Inpho software suite. Adjusted images were then draped upon a ground model and orthorectified. Individual orthorectified tiffs were blended together to remove seams and corrected for any remaining radiometric differences between images using Inpho’s OrthoVista.
Orthophotography Creation
Digital orthophoto frames were created by using existing public domain elevation data sets, which were in turn combined with Level 3 processed imagery, aerotriangulation data, and government supplied airborne topographic LiDAR bare earth datasets of various vintages. This orthorectification process was done in Hexagon Geomedia on the three data sources (processed imagery, aerotriangulation data and surface data). Manual Seamlines were drawn in ArcMap on every frame. Then, using the grid created with in-house software, a set of “base” mosaicked tiles were created in ImageStation OrthoPro using a bilinear interpolation method on the three data sources (rectified imagery, aerotriangulation data and surface data).
Post-processing and Quality Control
Radiometry was verified by visual inspection of the digital orthophoto. Slight systematic radiometric differences may exist between adjacent orthoimage files due primarily to differences in source image capture dates and sun angles along flight lines. These differences can be observed in an image's general lightness or darkness when it is compared to adjacent orthoimage file coverages.
All GeoTIFF data and image file sizes were validated using commercial GIS software to ensure proper loading before delivery to MassGIS.
Seamlines and tile edges were visually inspected. Seamline mismatches were corrected unless the overall displacement was less than one pixel. Orthoimages were visually inspected for completeness to ensure that no gaps or image misplacements existed within and between adjacent images.
In order to post-process the imagery data to meet task order specifications and ASPRS horizontal accuracy guidelines, the horizontal accuracy was tested by NV5 Geospatial using a total of 64 ground control points.
NV5 Geospatial created independent check points, none of which were used in the calibration or post-processing of the imagery data. The checkpoints were distributed throughout the project area and were surveyed using GPS techniques. The x and y positions of where the checkpoints fell on the orthoimagery were collected. These values were then compared to the surveyed control point x and y values.
Accuracy Statement for Zone 18:
Tested 0.334 meters (1.1 feet) RMSEx. The RMSEx of the imagery was calculated using 12 independent checkpoints.
Tested 0.192 meters (0.63 feet) RMSEy. The RMSEy of the imagery was calculated using 12 independent checkpoints.
Tested 0.385 meters (1.26 feet) RMSEr. The RMSEr of the imagery was calculated using 12 independent checkpoints.
Accuracy Statement for Zone 19:
Tested 0.192 meters (0.63 feet) RMSEx. The RMSEx of the imagery was calculated using 52 independent checkpoints.
Tested 0.175 meters (0.57 feet) RMSEy. The RMSEy of the imagery was calculated using 52 independent checkpoints.
Tested 0.26 meters (0.85 feet) RMSEr. The RMSEr of the imagery was calculated using 52 independent checkpoints.
Additionally, MassGIS collaborated with MassDOT’s Survey Section on an independent test of spatial accuracy. MassDOT provided surveyed locations for 80 photo identifiable point locations statewide. Example photo identifiable points include catch basin corners and corners of striping painted on tennis courts. Once the orthoimagery was delivered, MassGIS recorded coordinates from the imagery for each of the 80 photo-identifiable points. Coordinates for the photo identifiable points as determined from the imagery were compared to the corresponding surveyed coordinates and a root-mean-square error was calculated. This independent check confirmed that the imagery met the project’s horizontal accuracy specification.
JPEG 2000 Compression
MassGIS converted each GeoTIFF to an 8-bit JPEG 2000 (GMLJP2) format image using GeoExpress 10 software. These files maintain all four bands and were produced at a lossy compression ratio of 20:1, resulting in very minimal loss of image quality. Each image is approximately 19 MB, reduced from the 763 MB GeoTIFF, more suitable for distribution and download. MassGIS calculated statistics on the JPEG 2000 images in ArcGIS Pro. All of the imagery is available for public use and distribution.
When displaying the JPEG 2000 images, you can choose how to render them by assigning different band combinations to the red, green and blue channels. The default method is to assign bands 1, 2 and 3 to the red, green and blue channels, respectively, to display the natural color view. To display a color infrared view, assign bands 4 (near infrared), 1 and 2 to the red, green and blue channels, respectively. The difference in rendering is shown here:
2023 imagery basemap choices in MassMapper (top) and MassGIS’ ArcGIS Online maps and apps
Web Service
MassGIS created a tile service (cache), hosted at ArcGIS Online, from the JPEG 2000 files. This service type contains pre-rendered images to support fast visualization of large datasets. The cache uses the RGB bands and appears at scale levels 7 (1:4,622,324) to 20 (1:564). It is compatible not only with ArcGIS and web apps built with the ArcGIS APIs, but also QGIS and third-party apps that use OGC protocols such as WMTS.[Source] The cache appears as a basemap choice in MassMapper and all ArcGIS Online maps and apps hosted by MassGIS, including the Massachusetts Interactive Property Map.
In addition to the information below, we have created a step-by-step, illustrated guide to help you connect to the web service.
ArcGIS desktop software
ArcMap and ArcGIS Pro users can add data from ArcGIS Online and search for 'MassGIS 2023'.
QGIS
To add the service to QGIS, follow these instructions to create a new WMS/WMTS connection and for the URL enter:
https://tiles.arcgis.com/tiles/hGdibHYSPO59RG1h/arcgis/rest/services/orthos2023/MapServer/WMTS/1.0.0/WMTSCapabilities.xml
Use and Restrictions
No restrictions apply to these data. However, users should be aware that temporal changes may have occurred since this dataset was collected and that some parts of these data may no longer represent actual surface conditions. Users should not use these data for critical applications without a full awareness of their limitations.
Acknowledgement of MassGIS would be appreciated for products derived from these data.