Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  596 1625 (N19E44 L=182)   689    155        1 ( 554/1571- 638/1692) -24.8 E51-E38 N15-N23
      2  443 1623 (N16E57 L=169)   706    163        4 ( 403/1551- 482/1676) -77.5 E59-E54 N11-N20
      3  315 1527 (N08E64 L=162)  1808    192        7 ( 259/1454- 356/1603) -45.0 E71-E56 N03-N13
      4  510 1460 (N07E43 L=184)   281     70        1 ( 490/1429- 538/1484) -31.0 E44-E39 N05-N08
      5  374 1267 (S07E47 L=179)   337     78        0 ( 335/1244- 411/1288)   0.0 E50-E43 S07-S05
      6  322 1212 (S11E50 L=176)   340     54        0 ( 295/1187- 352/1232) -21.8 E53-E47 S11-S09
      7  896 1161 (S03E10 L=216)   189     38        0 ( 870/1144- 918/1176)   0.0 E11-E08 S03-S01
      8  618 1061 (S14E25 L=202)   345    110        1 ( 582/1012- 640/1111)  90.0 E26-E23 S17-S11
      9  440  983 (S22E36 L=190)   753    279        2 ( 346/ 913- 542/1075) -12.5 E45-E27 S24-S17
     10  218  890 (S31E58 L=168)   406     48        1 ( 200/ 859- 240/ 916) -63.4 E60-E55 S33-S29
     11  202  787 (S38E62 L=164)   500     38        0 ( 171/ 776- 234/ 801)   5.2 E68-E57 S38-S37
  

    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