Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  405 1626 (N13E62 L=259)   686    113        3 ( 371/1598- 439/1656)  18.4 E65-E59 N11-N15
      2  363 1564 (N09E63 L=259)   590    113        1 ( 340/1540- 394/1593) -14.0 E66-E58 N07-N10
      3  314 1216 (S13E51 L=270)  1138    211        1 ( 242/1155- 377/1288)  48.0 E56-E47 S17-S07
      4  274 1060 (S23E50 L=271)   865    195        0 ( 216/1013- 333/1111)  41.2 E55-E46 S25-S18
      5 1721  954 (N11W45 L=  6)   359     58        0 (1702/ 922-1740/ 987) -56.3 W42-W46 N09-N12
      6  330  806 (S37E42 L=279)   511    142        1 ( 291/ 762- 379/ 836)  45.0 E45-E37 S40-S33
      7 1809  680 (N00W65 L= 26)   519    114        1 (1790/ 650-1829/ 724) -18.4 W61-W68 S01-N01
      8  238  666 (S49E53 L=268)  1052    251        2 ( 185/ 649- 289/ 704)   9.5 E59-E46 S50-S45
      9 1610  613 (S10W49 L= 10)   355    111        1 (1572/ 596-1659/ 627)  23.2 W46-W53 S11-S08
     10 1696  562 (S10W60 L= 21)   533    120        5 (1673/ 528-1717/ 606) -53.1 W56-W62 S11-S07
     11 1387  356 (S31W48 L=  9)  1537    387        1 (1247/ 305-1492/ 435)   4.8 W34-W58 S33-S28
  

    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