Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1529 1531 (N25W49 L=199)  1695    261        2 (1425/1443-1609/1605) -66.8 W44-W52 N18-N32
      2  403 1434 (N44E44 L=106)   510     75        3 ( 374/1409- 446/1454) -16.7 E48-E38 N41-N45
      3 1165 1300 (N20W14 L=164)   283    132        0 (1125/1249-1216/1341) -60.3 W12-W16 N16-N23
      4 1812 1203 (S01W59 L=209)   413      9        0 (1791/1182-1831/1220)   0.0 W56-W62 S01-N00
      5 1575 1003 (S06W34 L=184)   481    116        2 (1521/ 970-1625/1034) -35.5 W30-W37 S08-S03
      6 1660  917 (S13W39 L=189)   396    103        0 (1637/ 866-1675/ 970)  71.6 W36-W40 S15-S09
      7  624  910 (N07E27 L=123)   223     34        0 ( 599/ 894- 651/ 929)  14.0 E28-E24 N06-N07
      8  282  703 (S01E60 L= 90)   927    101        2 ( 250/ 655- 309/ 745) -78.7 E62-E56 S04-N01
      9 1697  686 (S29W41 L=191)   946    137        1 (1644/ 643-1744/ 738)  12.5 W36-W45 S31-S26
     10  273  554 (S11E70 L= 80)   744     76        1 ( 246/ 522- 295/ 577) -63.4 E72-E66 S13-S09
     11 1207  516 (S28W01 L=151)   257     63        0 (1182/ 484-1227/ 545)  45.0 W00-W03 S30-S26
     12  458  508 (S16E51 L= 99)   415     71        0 ( 442/ 473- 472/ 543)  45.0 E52-E48 S17-S13
  

    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