Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  833 1899 (N69W01 L=328)  1051     99        0 ( 807/1853- 873/1931)  48.0 E03-W10 N62-N75
      2  670 1438 (N30E18 L=309)   892    230        3 ( 559/1401- 731/1479)  -8.7 E26-E13 N27-N33
      3 1176 1386 (N20W15 L=342)   718    289        0 (1064/1281-1243/1462) -61.7 W08-W17 N12-N25
      4 1808 1416 (N13W66 L= 33)  1222    140        3 (1751/1372-1860/1452) -17.1 W59-W73 N09-N14
      5  317 1251 (N23E47 L=280)   517    114        2 ( 300/1198- 336/1306)  80.5 E48-E44 N19-N26
      6 1195 1087 (N02W11 L=338)   307     68        1 (1173/1054-1220/1124)  45.0 W09-W12 N00-N03
      7 1549  566 (S36W32 L=359)  2890    483        6 (1356/ 442-1723/ 696)  17.9 W15-W47 S42-S28
      8  483  496 (S24E48 L=279)  3232    524        4 ( 388/ 299- 614/ 688)  76.0 E53-E42 S39-S11
      9  813  261 (S47E35 L=292)   397    112        0 ( 755/ 246- 866/ 275) -24.0 E39-E30 S48-S44
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y