def proj_invsin_tile(iH, iV, fLin, fCol, iSize): """ Convert MODIS sinusoidal tile coordinates to latitude and longitude. Parameters: iH (int): Horizontal tile number. iV (int): Vertical tile number. fLin (float): Line (row) index within the tile. fCol (float): Column index within the tile. iSize (int): Dimension of one MODIS tile (e.g., 1200 for 1km, 2400 for 500m, 240 for 5km). Returns: tuple: (longitude, latitude) in degrees. """ y = (9 - iV - (fLin + 0.5) / iSize) / 9.0 fai = y * np.pi / 2.0 latitude = y * 90.0 x = (iH + (fCol + 0.5) / iSize - 18.0) / 18.0 lambuda = x / np.cos(fai) longitude = lambuda * 180.0 return longitude, latitude