Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  480 1648 (N19E55 L=183)   961    130        2 ( 455/1581- 502/1707) -33.7 E63-E49 N14-N22
      2  366 1546 (N10E59 L=179)   624     87        1 ( 340/1497- 389/1587) -56.3 E61-E56 N07-N13
      3 1211 1350 (N16W03 L=241)   212     54        1 (1189/1320-1239/1371) -18.4 W01-W04 N13-N16
      4 1857 1283 (N32W56 L=294)   569     80        0 (1833/1242-1881/1317) -31.0 W52-W59 N29-N33
      5  248 1294 (S07E60 L=178)   627    101        0 ( 197/1273- 298/1311)   0.0 E66-E54 S07-S05
      6  207 1228 (S11E63 L=175)   385     82        0 ( 177/1205- 237/1246)   7.1 E66-E58 S12-S09
      7  185 1005 (S24E62 L=176)   786    121        0 ( 149/ 982- 232/1039)  16.7 E67-E56 S25-S22
      8  344  983 (S24E44 L=194)   390     82        0 ( 316/ 948- 379/1022) -26.6 E47-E40 S25-S21
      9  909  138 (S70W46 L=284)   541     68        1 ( 888/ 109- 931/ 167)  -5.0 W37-W62 S72-S66
  

    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