Here is a little something I cooked up, though it's not a perfect fit to your guidelines. I'll describe my methodology, as it could prove useful to others.
First, I went to the
USGS EarthExplorer site. I used to use seamless.usgs.gov, but that seems to have gone away
. To download data, you need to be logged in, but setting up a profile is free. So, why not?
Under the Search Criteria tab, I selected the Address/Place tab and entered "stewart island" and tapped Show. Click on "Stewart Island / Rakiura, Southland 9818, New Zealand". Zoom in a lot and tap Clear Coordinates. Click on the map to define the outline of the area you want to use. Now select the data sets tab. In the hierarchical list click on Digital Elevation, then SRTM. Under SRTM check SRTM Void Filled and SRTM Water Body Data. Then click on the Results tab. You should get four elevation data sets, and four water body data sets. Click on the footprints for each of these to insure that they cover the desired area., if they're where they should be click on the download button for each data set in the SRTM water body data. Then click on each of the SRTM void filled elevation data sets, in the resulting window, click the Download button next to TIFF, although it looks like Wilbur can handle DTED.
Now that you have the data, this next step will require QGIS. You can get QGIS free
here.
Next add the four raster layers to QGIS.From the main menu, select Raster>Miscellaneous>Merge…
Select the four input files, then browse a location to place the resulting merged elevation file. I named it merged_elevations.bt, and selected VTP .bt as the type. Now, I like to double click on the new layer in the legend area to open properties. In the Style tab, under load min/max values, I select Min/Max and Actual(slower), then hit Load and OK. This just makes the display look nicer.
Next, I'm going to clip the coverage area of the final raster to just the island. Start by clicking the New Shapefile Layer tool. This will create an empty new vector graphics object. Give it a name and set its type to polygon. We will fill that with the boundaries of our desired area. With your new shapefile selected, tap Toggle Editing in the Digitizing toolbar. This will allow you to define the area of the island's raster. Click the Add Feature tool in the Digitizing toolbar. Now click in the water areas around the island to create a new polygon encircling the island. As a guide, you can create a "landarea" raster using Raster Calculator under the raster menu. Simply set the expression to "merged_elevations > 0.0". Place that below the shapefile layer as a guide. Tap Toggle Editing again to finish your polygon.
Now pull down Raster>Extraction>Clipper… to actually clip the elevations. Set the input file to your merged_elevations, set the Clipping Mode to Mask Layer with the shapefile you just created as the Mask Layer. Name the resulting file to something like "clipped_elevations".
Open
Wilbur and open up clipped_elevations.bt.
Use the paintbrush with value set to something tiny like 0.001 and the Operation set to Maximum to paint over areas below zero within the island. These may be a fault in Wilbur, I'm not sure. I don't see them in QGIS.
Next Surface>Locate>Flip Vertically to get roughly the same arrangement as has been used previously in this thread.
We'll need to use an exponent operator to contract the mountain areas. Realistically, most continents would have a bit less highland, so that's what I'm going to do here. Select>From Terrain>Height Range… from Minimum: -1 to Maximum: 0. With that selection, Filter>Fill>Set Value… Set Value: -1. Otherwise the exponential operation tends to mess things up. Deselect. Next, Filter>Mathematical>Exponent… Set the exponent to 2 for the Land(above sea level), and 0.5 for the Sea(below sea level), 0 for the Sea Level. Preserve Height should be set to Absolute Low set to -1, Absolute High set to 791.
Use rectangle select to select an area very close all around the island. Now, Surface>Crop to Selection.
Finally, to place this where it belongs, sort of, Surface>Find Min/Max… Top: 60, Left: 80, Right: 140, Bottom: 0. File>Save As… set the type to Binary Terrain Surface(*.bt). Name it something like "squared_elevations".
Back to QGIS. While we used Wilbur to assign the new coordinates which were stored in the VTP .BT format, but the projection is hardwired as UTM. So, when the elevations file is loaded we need to reassign a WGS 84 latitude-longitude "projection"(Raster>Projections>Assign Projection… elevations.bt as the Input file, EPSG:4326 as the Desired SRS. This will be an equirectangular "latlon" coordinate system). Looking back on it, I think I can manage the projection described in this thread using georeferencing and probably a polar equidistant projection.
Finally, I'm going to use the Raster Calculator one more time. I'll use the expression: 4*squared_elevation + 2.375*unsquared_elevations to give me a range. That will give us our final elevations.
For now, I'm going to leave it at that. There's a lot more you can do, particularly when you've created a vector representation of the landmasses.