Jump to content

File:Some sumerian city states with elevated sea level 1.png

Page contents not supported in other languages.
This is a file from the Wikimedia Commons
From Wikipedia, the free encyclopedia

Original file (2,570 × 1,965 pixels, file size: 2.71 MB, MIME type: image/png)

Summary

Description
English: Some Sumerian cities and towns
Date
Source Own work
Author Merikanto

This map is based on SRTM dem and knowlegde about cities of Sumer

GMT 6.3 code to make map


    1. map of sumer
    1. GDAL and GMT - Generic Mapping Tools Bash script
  5. 2024.07.17 0000.0003c


  1. !/bin/bash


  1. DEM0="./dem1/near_east_22.tiff"

DEM0="./dem1/near_east_1.tiff"


  1. DEM file

DEM1="./dem1/dem4.tif"

  1. Output file name

OUTPUT1="sumerian_cities_map_1.png"


paletti="paletti.cpt"

echo "-9000 100/100/200 00 100/200/255" > $paletti echo "00 80/200/100 10 100/200/100" >> $paletti echo "10 100/200/100 40 150/255/150" >> $paletti echo "40 150/255/150 100 255/255/150" >> $paletti echo "100 255/255/150 1000 255/200/80" >> $paletti echo "1000 255/200/80 2000 200/180/70" >> $paletti echo "2000 200/180/70 4000 185/130/55 " >> $paletti echo "4000 185/130/55 8000 185/130/55 " >> $paletti echo "B black" >> $paletti echo "F white" >> $paletti echo "N 128" >> $paletti

  1. Create a color palette for the land heights

gmt makecpt -Cetopo1 -T0/2000/1 > heights.cpt


    1. pre-process rasters
  1. gdal_calc.py -A $DEM0 --outfile=$DEM1 --calc="A-4" --NoDataValue=0 --overwrite
  2. gdal_calc.py -A $DEM1 --outfile=./dem1/sea_mask.tif --calc="A<=0" --NoDataValue=0 --overwrite
  1. gdal_polygonize.py ./dem1/sea_mask.tif -f "ESRI Shapefile" ./dem1/sea_mask.shp
  1. Create shaded relief from the DEM
  2. gmt grdgradient $DEM1 -Ne0.6 -A70 -G./dem1/region_dem_shade.nc


  1. echo "DEBUG"
  2. quit(-1)


  1. Start GMT session

gmt begin $OUTPUT1 png -C

  1. Plot the DEM with the land color palette and shaded relief

gmt grdimage -R43/49/30/34 -JM8i $DEM1 -Cheights.cpt -I./dem1/region_dem_shade.nc -B1


  1. gmt grdimage -R43/49/30/34 -JM8i $DEM1 -Cheights.cpt -B1

gmt grdcontour $DEM1 -R43/49/30/34 -JM8i -C+0 -W1.0p,blue -t50

  1. Draw the coastlines and borders without filling the sea

gmt coast -R43/49/30/34 -JM8i -W1/0.5p -Df

  1. Draw the rivers on top of everything

gmt coast -R43/49/30/34 -JM8i -Ia/0.25p,blue

  1. Draw the sea areas using the shapefile

gmt psxy ./dem1/sea_mask.shp -R43/49/30/34 -JM8i -Gblue


cat << EOF > zikkurats.txt 45.996111;30.815833;Eridu 45.637222;31.324167;Uruk 46.105126;30.961653;Ur 45.853611;31.285833;Larsa 45.233381;32.126444;Nippur 44.341689;32.391997;Borsippa 44.604722;32.540278;Kish 44.421111;32.5425;Babylon 44.252153;33.058829;Sippar 44.202222;33.353611;Dur-Kurigalzu 48.257778;32.190556;Susa 48.520833;32.008333;Chogha Zambil 43.152778;36.359444;Ninive 43.2625;35.456667;Assur 43.229444;36.509444;Dur-Sharrukin 42.449336;36.257086;Tell Rimah 43.328889;36.098056;Kalhu 43.270008;35.494685;Kar-Tukulti-Ninurta EOF

  1. Define city coordinates and names

cat << EOF > cities.txt 46.025556;30.97;Ubaid 45.996111;30.815833;Eridu 45.637222;31.324167;Uruk 44.779439;32.717781;Dzemdet Nasr 46.407222;31.411389;Lagash 44.604722;32.540278;Kish 45.233381;32.126444;Nippur 46.105126;30.961653;Ur 45.853611;31.285833;Larsa 44.252153;33.058829;Sippar 45.999722;31.379722;Bad-Tibira 45.509722;31.7775;Shuruppak 46.177583;31.562028;Girsu 45.8875;31.667222;Umma 46.495728;31.376817;Nigin 45.968889;31.946944;Adab 45.480556;31.838056;Kisurra 45.876667;31.743333;Zabalam 44.202222;33.353611;Dur-Kurigalzu 45.083333;32.266667;Abu Salabikh 45.660833;32.3125;Tell al-Wilayah 44.728333;33.484167;Eshnunna 44.555753;33.354675;Tutub 44.664722;32.781667;Tell Uqair 44.99;34.38;Tell Madhur 44.664722;32.781667;Urum 44.612861;32.760028;Kutha 44.341689;32.391997;Borsippa 44.421111;32.5425;Babylon 44.783333;32.092778;Marad 45.931389;33.123611;Der 47.147778;32.395833;Pashime 48.257778;32.190556;Susa 48.520833;32.008333;Chogha Zambil 48.940278;32.171111;Dehloran 48.546111;32.209167;Chogha Mish 52.428056;30.096667;Banesh 52.407778;30.011667;Anshan 40.89;34.549444;Mari 42.826389;33.645;Hit 40.569444;34.923333;Terqa 44.254917;35.369972;Nuzi 41.058644;36.667617;Tell Brak 43.197778;36.371389;Tell Arpachiya 45.44999;33.88333;Choga Mami 41.955833;36.815278;Tell Hamoukar 43.152778;36.359444;Ninive 43.2625;35.456667;Assur 40.997222;37.056944;Urkesh 41.505278;36.957222;Tell Leilan 39.025;36.870833;Harran 38.111375;35.986842;Emar 36.798;35.798;Ebla 48.243611;34.181389;Tepe Gijan 51.404722;33.968889;Tepe Sialk 43.229444;36.509444;Dur-Sharrukin 42.449336;36.257086;Tell Rimah 43.328889;36.098056;Kalhu 43.270008;35.494685;Kar-Tukulti-Ninurta EOF


cat << EOF > channel1.txt 44.216335227272737 32.875695197020214 44.350568181818183 32.804804524206197 44.582599431818146 32.531186066910536 44.728338068181785 32.404992055622998 44.816548295454496 32.275378782179224 45.006392045454554 32.138273382582348 45.201988636363645 32.076549200436666 45.3515625 31.982258220124343 45.453196022727241 31.839825403047456 45.593181818181826 31.721642790703573 45.620028409090914 31.519160239230526 45.598934659090872 31.343264975336034 45.66796874999995 31.208872123353601 45.733167613636333 31.12683137113526 45.823295454545459 31.030753567237866 45.876988636363642 30.899208093312126 45.957528409090919 30.820193902545316 46.137784090909093 30.774072239125214

  1. 46.496377840909055 30.814429904955798
  2. 46.615269886363599 30.804547961537182
  3. 46.789772727272734 30.83666056084871
  4. 46.899076704545408 30.873700212735312
  5. 46.985369318181831 30.928821621356441

EOF

cat << EOF > channel2.txt 44.48863636363636 33.178755356829228 44.795454545454554 32.949747606535105 44.891335227272734 32.86281019518978 44.849147727272737 32.795133228954761 44.956534090909102 32.70158983033182 45.232670454545463 32.514208715198968 45.635369318181823 32.420369203787672 45.738920454545465 32.28430329379534 45.71207386363637 32.141520867733156 45.926846590909101 31.952975772653783 45.980539772727283 31.806423567968839 45.980539772727283 31.702066969344013 46.003551136363647 31.502811816730752 45.923011363636377 31.290022424337156 46.16079545454545 30.995418531289268

  1. 46.367897727272727 30.850654999487432
  2. 46.536647727272729 30.765010026837757

EOF


cat << EOF > channel3.txt 45.934517045454555 31.718380440814656 46.064914772727278 31.626988054809161 46.191477272727276 31.591059307740046 46.352556818181824 31.499541789078343 46.429261363636378 31.404661263862025 46.544318181818184 31.339170463330959 46.598011363636374 31.211332247846421 EOF

gmt psxy channel1.txt -R43/49/30/34 -JM8i -W0.5,blue gmt psxy channel2.txt -R43/49/30/34 -JM8i -W0.5,blue gmt psxy channel3.txt -R43/49/30/34 -JM8i -W0.5,blue

gmt psxy cities.txt -R43/49/30/34 -JM8i -Sc0.3c -Gred -W0.5p,black -t20 gmt pstext cities.txt -R43/49/30/34 -JM8i -F+f12p,Helvetica,black+jLB -t20 -D0.2c/-0.2c

gmt psbasemap -R43/49/30/34 -JM8i -Lx1.5i/+0.5i+c12+w200k+l"Scale (km)"+f

  1. --LABEL_FONT_SIZE=12


  1. End GMT session

gmt end




R code do download SRTM map

library(elevatr) library(terra)

  1. extent1 <- extent(38, 42, 32, 35) # xmin, xmax, ymin, ymax
    1. manually specify corners of the bounding box of the US

bbox_usa <- data.frame(x = c(30,60), 

                      y = c(25, 45))
    1. turn into spatial, projected bounding box

sf_bbox_usa <- sf::st_as_sf(bbox_usa, coords = c("x", "y"), crs = 4326)

elev_usa <- get_elev_raster(locations = sf_bbox_usa, z = 8)

terra::plot(elev_usa) terra::writeRaster(elev_usa, "near_east.tiff", overwrite=TRUE)

Licensing

I, the copyright holder of this work, hereby publish it under the following license:
w:en:Creative Commons
attribution
This file is licensed under the Creative Commons Attribution 4.0 International license.
You are free:
  • to share – to copy, distribute and transmit the work
  • to remix – to adapt the work
Under the following conditions:
  • attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.

Captions

Some Sumerian cities and towns

Items portrayed in this file

depicts

inception

13 July 2024

media type

image/png

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current07:00, 16 July 2024Thumbnail for version as of 07:00, 16 July 20242,570 × 1,965 (2.71 MB)MerikantoSea, new code
11:13, 15 July 2024Thumbnail for version as of 11:13, 15 July 20241,952 × 1,510 (2.47 MB)MerikantoUpdate of layout
14:30, 13 July 2024Thumbnail for version as of 14:30, 13 July 20241,952 × 1,510 (2.5 MB)MerikantoUploaded own work with UploadWizard

The following page uses this file:

Metadata

This file contains additional information, probably added from the digital camera or scanner used to create or digitize it.

If the file has been modified from its original state, some details may not fully reflect the modified file.

Retrieved from "http://en.wiki.x.io/wiki/File:Some_sumerian_city_states_with_elevated_sea_level_1.png"