Posts

3D Toronto Images

Below are a few 3D Toronto images from a demonstration that I gave comparing Esri Arc Scene with FLY in PCI Geomatica. I generated the digital surface model (DSM) from some demo LIDAR all hits data that we had. The coverage area is for a small portion of downtown Toronto centered around Toronto City Hall.

Toronto City Hall 3D LIDAR image made with ESRI Arc Scene

Toronto City Hall 3D LIDAR image made with PCI Geomatica

LIDAR all hits point data of downtown Toronto

3D perspective view with LIDAR all hits point data of downtown Toronto

PCI EASI Script [PIX_to_SHP_batch_export.eas]

This simple EASI script used with the MODEL command in PCI will batch convert vector files stored in PCIDSK (PIX) format into Shape File (SHP) format with the exact same file name as the input files. This was originally created for PCI Geomatica v9.1 but was last tested and working with no problems in PCI Geomatica v10.0


!—————————————————————————-
!—————————————————————————-
!
! Batch Export PIX Vector to SHP Script
! [PIX_to_SHP_batch_export.eas]
!
! This script will export vector segments from PIX files located in a
! given directory into Shape files and place the new files into the same
! directory. The script assumes that no SHP files with the same names
! already exist and all the input files are setup the same with the
! vector segments stored the same.
!
!—————————————————————————-
!—————————————————————————-
! Define variables
!—————————————————————————-

!to store location of input & output files

local string in_files

!for the file format and extension types

local string type, ext

!file names

local string bn, fn

!to store directory listing of input files

local mstring dirlist

!to store vector segment number

local integer vec

local integer i

local $Z

!—————————————————————————-
! Clear the EASI window and then show the header information
!—————————————————————————-

PRINT @(1 ,1,CLREOS)

print “———————————————————————–”
print @reverse,” ‘Batch export PIX vector to SHP’ EASI Script “,@alloff
print “”
print “”
print “”
print “———————————————————————–”
print “”
print ” This script will export vectors from PIX files in a given directory ”
print ” into SHP format using the same file names as the input files. ”
print “”
print ” All input vector files are expected to be setup the same with all ”
print ” vectors to be exported, stored using the same segment number.”
print “”
print “———————————————————————–”
print “”

!—————————————————————————-
! Collect input from user
!—————————————————————————-

print “Enter the directory that contains the PIX files to export to SHP:”

input “>” in_files

print “Enter the vector segment number that the vector is stored in:”

input “>” vec

!——————————————————————-
! Get the contents of the directory
!——————————————————————-

type =”pix

dirlist = getdirectory(in_files)

let $Z = “\

for i = 1 to f$len(dirlist)

!——————————————————————-
! Extract parts of the filenames
!——————————————————————-

fn = in_files + $Z + dirlist[i]

ext = getfileextension(fn)

bn = getfilebasename(fn)

if (ext ~= type) then

print “”
print “Exporting:”, bn, “from PIX to SHP”
print “”

!—————————————————————-
! Set up the parameters and execute the FEXPORT command
!—————————————————————-

fili = in_files + $Z + dirlist[i]
filo = in_files + $Z + bn
dbiw =
dbic =
dbib =
dbvs = vec
dblut =
dbpct =
ftype =”SHP
foptions =

R Fexport

PRINT @(1 ,1,CLREOS)

endif

endfor

PRINT @(1 ,1,CLREOS)

print “——————————————————————-”
print “”
print @reverse,” ‘Batch export PIX vector to SHP’ EASI Script Finished “,@alloff
print “”
print “——————————————————————-”

return

!—————————————————————————-
!—————————————————————————-

Script to clip the unnecessary excess portions of an air photo mosaic

Below is a simple EASI script that when used with the MODEL command in PCI will clip the unnecessary excess portions of an air photo mosaic created from Ortho Engine to an irregular buffered shape around a set study area. This was originally created for PCI Geomatica v9.1 but was last tested and working with no problems in PCI Geomatica v10.0 – Just copy the code below into a blank text file and edit as needed.

!——————————————————————————————–

! Irregular_Polygon_Clip [bitmap_clip.eas]

! Ted MacKinnon – tmackinnon.com

! This simple script used with the MODEL command in PCI EASI

! will clip the unnecessary excess portions of an airphoto mosaic

! to an irregular buffered shape around the study area.

!

! The working file ‘working-file.pix’ has the existing

! mosaic image located in the first three channels,

! an existing irregular shaped polygon bitmap and

! three empty 8 bit channels.

!

! %%2 is the bitmap mask of the irregular shape

! %1, %2, %3 are the RGB existing image channels

! %4, %5, %6 will be the new modeled RGB image channels

! the RGB value of 255, 255, 255 will set the background to white

!

! Simply change the file name to use this script with a different

! file and ensure that the channels and segments are setup the same

!——————————————————————————————–

!——————————–

! Set up and run the model

!——————————–

MODEL ON “working-file.pix” OVER dbiw

if %%2 = 1 then

%4 = %1;

%5 = %2;

%6 = %3;

else

%4 = 255;

%5 = 255;

%6 = 255;

endif;

ENDMODEL

!———————————————-

! Export the resultant channels to a new file

!———————————————-

FILI = “working-file.pix

FILO = “mosaic_clipped.tif

DBIW =

DBIC = 4,5,6

DBIB =

DBVS =

DBLUT =

DBPCT =

FTYPE = “TIF

FOPTIONS = “”

RUN FEXPORT

!——————————————————————————————–

 

Color Shaded Relief Model of Vancouver Island

Color Shaded Relief Model of Vancouver Island, British Columbia

The above image is a Color shaded relief (CSR) model created using a DEM of Vancouver Island, British Columbia using PCI Geomatica software. The change in relief ranges from 0 m (blue) up to over 2200 m (red colors). Vancouver Island is a large island in the province of British Columbia off the Pacific coast of Canada, it is the largest island (32,134 square kilometers) on the western side of both North and South America.

This image has been published in various books, University Course material (such as UVic) and papers such as the Vancouver Sun & the Montreal Gazette

More examples of Color shaded relief (CSR) models …

Custom Pansharp Fused Image

Parc de la Gatineau - Pansharp Image

Above is a custom poster I created of Parc de la Gatineau in Quebec (just north of Ottawa) for Geocognition Consulting. It is a pansharpened image fused with a DEM, created using PCI Geomatica. The fused elevation data gives the image an extra 3D effect making the Park pop out at the viewer.

Enhanced LandSAT Images of Trout Lake


The following images are enhanced LandSAT images of Trout Lake, Nova Scotia (popular fishing place on the South Mountain – south of the town of Middleton) created using PCI Geomatica. The first image is a True Color Composite (TCC) showing bands 3,2 an 1. The second one is a False Color Composite (TCC) showing bands 7,4 and 3 and helping to clearly define the various clear cuts in the area. [Both images beloware actually low res versions of the original images]

 

Enhanced LandSAT True Color Composite Image of Trout Lake, Nova Scotia

 

Enhanced LandSAT False Color Composite Image of Trout Lake, Nova Scotia

Color Shaded Relief Models – Cap Pele, New Brunswick

Color Shaded Relief Models of Cap Pele, New Brunswick Color Shaded Relief Models of Cap Pele, New Brunswick

The two images above represent artificial three dimensional perspective views perspective views from different points of origin featuring color shaded relief models that were created from high resolution LIDAR digital surface models as part of a LIDAR flood modeling graduate research project.

The study area for the project consisted of the coastal Gulf Shore region of southeastern New Brunswick from Kouchibouguac National Park south to Jourimain Island (location of the Confederation Bridge). The coastal zone of New Brunswick is a picturesque fishing region that boasts several kilometers of sandy beaches with some of the warmest salt water temperatures north of Virginia.