Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  832 1681 (N32E33 L=272)  1319    311        1 ( 739/1564- 930/1765) -19.8 E46-E21 N26-N35
      2  583 1622 (N22E47 L=258)   776    215        0 ( 529/1548- 612/1684) -33.7 E55-E42 N16-N25
      3 1404 1633 (N46W10 L=315)   497     92        1 (1369/1612-1449/1649)   7.1 W06-W14 N43-N47
      4  630 1526 (N17E38 L=267)   670    233        1 ( 565/1439- 685/1611)  85.6 E39-E36 N10-N23
      5  217 1285 (S06E65 L=240)   749     87        4 ( 199/1252- 244/1333) -18.4 E69-E60 S07-S02
      6 1788 1195 (N27W50 L=355)   826    111        1 (1754/1145-1814/1248) -45.0 W46-W53 N24-N29
      7 1871 1144 (N26W61 L=  6)   679    135        0 (1837/1092-1906/1200) -24.4 W56-W67 N24-N29
      8 1562  615 (S13W46 L=351)   340     91        1 (1525/ 584-1590/ 640)   9.5 W42-W48 S14-S11
      9 1301  561 (S23W30 L=335)   267     88        0 (1288/ 526-1320/ 612) -45.0 W27-W32 S24-S19
     10 1412  538 (S21W39 L=344)   254     71        0 (1396/ 510-1430/ 573) -78.7 W37-W39 S23-S18
  

    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