Home › Roadmap › GIS Fundamentals › Part 9

Rasterio & Xarray

📘 পাইথনে রাস্টার ডেটা পড়া ও বিশ্লেষণ।

🌐 Read and analyse rasters in Python.

🗺️ জিআইএস মৌলিক🎚 Advanced📍 Part 9 / 12🧩 3 সাব-টপিক
🗺️
1

ধারণা · Concept

কী এবং কেন
📘 বাংলা ব্যাখ্যা:
পাইথনে রাস্টার ডেটা পড়া ও বিশ্লেষণ।
🌐 English:
Read and analyse rasters in Python.

🗂 কোথায় পড়ে

৫. জিওএআই → জিআইএস মৌলিক → Rasterio & Xarray

🎯 এই শাখার লক্ষ্য

Handle CRS, vector/raster data, spatial SQL and web maps

🧰 টুলস

QGIS, GeoPandas, PostGIS, Leaflet

🪜

ধাপে ধাপে শিখুন · Step-by-step

কেন শিখবেন → কীভাবে শিখবেন → উদাহরণ → সাধারণ ভুল
💡 কেন গুরুত্বপূর্ণ?

স্যাটেলাইট, জলবায়ু ও উচ্চতার রাস্টার ডেটা Python-এ পড়া ও বিশ্লেষণ।

  1. ধাপ 1
    rasterio.open("img.tif") — src.read(1), src.crs, src.transform, src.nodata
  2. ধাপ 2
    mask দিয়ে পলিগনে ক্লিপ; window দিয়ে বড় ফাইলের অংশ
  3. ধাপ 3
    rasterstats.zonal_stats দিয়ে প্রতিটি পলিগনে গড়/সর্বোচ্চ
  4. ধাপ 4
    xarray/rioxarray: সময় × অক্ষাংশ × দ্রাঘিমাংশের বহুমাত্রিক ডেটা (NetCDF জলবায়ু ডেটা)

🧪 বাস্তব উদাহরণ

ERA5 তাপমাত্রা NetCDF → xarray → প্রতিটি জেলার মাসিক গড় তাপমাত্রা।

⚠️ সাধারণ ভুল

nodata মান (−৯৯৯৯) বাদ না দিয়ে গড় হিসাব।

🔑 মূল শব্দ (Key terms)

rasteriotransformnodatazonal statisticsxarrayNetCDF

3

সাব-টপিক ও বিস্তারিত · Sub-topics

টপিক → সাব-টপিক → সাব-সাব টপিক
9.1 Rasterio9.1.1 ব্যান্ড পড়া9.1.2 মেটাডাটা
9.2 উইন্ডো ও মাস্ক
9.3 Xarray9.3.1 বহুমাত্রিক সময়-সিরিজ রাস্টার
🎬

ভিডিও · Videos (বাংলা ও English)

ক্লিক করলে সাইটের ভেতরেই ফ্লোটিং প্লেয়ারে চলবে · ভিডিও MySQL-এ সংরক্ষিত
💻

হাতে-কলমে · Try it (GIS Fundamentals)

এডিট করে Run চাপুন · Ctrl+Enter
import numpy as np, matplotlib.pyplot as plt
# Buffer & point-in-polygon without GIS libs
poly=np.array([[90.3,23.7],[90.5,23.7],[90.5,23.9],[90.3,23.9]])  # Dhaka bbox (lon,lat)
pts=np.random.default_rng(1).uniform([90.2,23.6],[90.6,24.0],(60,2))
inside=(pts[:,0]>=90.3)&(pts[:,0]<=90.5)&(pts[:,1]>=23.7)&(pts[:,1]<=23.9)
print("Points inside polygon:",inside.sum(),"/",len(pts))
plt.fill(*poly.T,alpha=.25); plt.scatter(*pts.T,c=np.where(inside,"#059669","#ef4444"),s=18)
plt.gca().set_aspect(1); plt.title("Point-in-polygon (spatial join idea)"); plt.show()
❓

সেল্ফ-চেক · Self check

1. “ব্যান্ড পড়া, মেটাডাটা” — এটি কোন সাব-টপিকের বর্ণনা?

2. “বহুমাত্রিক সময়-সিরিজ রাস্টার” — এটি কোন সাব-টপিকের বর্ণনা?

🔗

সংযুক্ত কনসেপ্ট · Connected concepts

← আগের টপিকGeoPandasপরের টপিক →Web Mapping Leaflet Folium
🎬 Video