Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1646 1626 (N23W64 L=183)   637    111        1 (1626/1595-1663/1666)  45.0 W61-W66 N20-N24
      2 1485 1246 (N08W33 L=152)   287     65        0 (1467/1210-1501/1282)  45.0 W31-W34 N06-N10
      3 1801 1239 (S02W59 L=178)   486    136        0 (1758/1201-1843/1264) -35.5 W54-W62 S05-N00
      4  172 1205 (N35E61 L= 58)   732    138        1 ( 144/1174- 212/1224)   0.0 E63-E56 N32-N35
      5 1742  736 (S30W42 L=161)  1824    319        2 (1636/ 643-1839/ 834)  13.4 W31-W52 S33-S25
      6 1205  656 (S19W02 L=121)   248     91        1 (1157/ 641-1252/ 667) -31.0 E01-W04 S20-S17
      7 1004  503 (S23E15 L=105)   510    133        0 ( 944/ 473-1075/ 553) -45.0 E16-E10 S26-S18
      8  766  445 (S22E32 L= 88)   320     70        0 ( 748/ 410- 781/ 472)  68.2 E33-E29 S24-S19
      9  487  435 (S17E53 L= 67)   458     95        0 ( 452/ 398- 513/ 468) -78.7 E54-E51 S19-S14
     10  555  381 (S22E50 L= 69)  1630    379        3 ( 502/ 269- 602/ 451)  65.0 E57-E44 S31-S16
  

    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