Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  925 1376 (N28E04 L=157)   319     77        1 ( 888/1353- 967/1391)  -9.5 E06-W00 N26-N29
      2 1547 1294 (N17W38 L=199)   344    127        2 (1491/1279-1601/1314) -20.6 W34-W42 N15-N18
      3  726 1036 (N08E19 L=142)   386    106        1 ( 694/1006- 770/1071) -21.8 E20-E15 N05-N10
      4 1906 1019 (S05W70 L=231)  1126     79        1 (1870/ 999-1948/1043) -12.5 W64-W82 S07-S02
      5 1813  833 (S15W57 L=218)   471    101        1 (1778/ 798-1848/ 854) -18.4 W53-W62 S17-S13
      6  726  707 (S13E22 L=139)   485    162        1 ( 659/ 680- 777/ 746)   7.1 E25-E17 S14-S10
      7 1362  610 (S24W20 L=180)  1953    494        6 (1203/ 439-1524/ 743)  38.0 W08-W31 S35-S16
      8 1752  631 (S28W55 L=216)   430     56        0 (1730/ 605-1774/ 657)  45.0 W52-W57 S29-S26
      9  603  340 (S39E44 L=117)   564     85        1 ( 587/ 294- 626/ 378)  35.0 E48-E38 S42-S35
     10 1077  271 (S49E04 L=157)   530     88        0 (1045/ 251-1116/ 302) -18.4 E06-W00 S51-S45
  

    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