Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  266 1485 (N11E71 L=206)  1565     83        4 ( 238/1437- 297/1519)  51.3 E79-E65 N07-N12
      2  306 1446 (N09E62 L=214)   577    105        1 ( 280/1423- 329/1494)  78.7 E65-E59 N07-N12
      3 1436 1404 (N30W19 L=296)   497    128        2 (1396/1342-1465/1442) -45.0 W15-W20 N26-N32
      4 1720 1034 (N14W46 L=323)   304     54        0 (1696/1001-1740/1063)  63.4 W44-W47 N11-N15
      5 1488 1032 (N08W28 L=305)   294     74        2 (1461/ 992-1509/1062) -26.6 W25-W29 N06-N09
      6  221  981 (S20E58 L=219)   453     47        0 ( 203/ 953- 240/1007)  63.4 E60-E55 S22-S18
      7 1777  960 (N10W53 L=329)   467    100        2 (1752/ 906-1794/ 998) -59.0 W50-W55 N07-N12
      8 1547  948 (N05W34 L=310)   470     69        2 (1514/ 922-1586/ 982)  31.0 W31-W36 N02-N06
      9  797  762 (S22E08 L=268)   224     44        1 ( 774/ 749- 826/ 784)  18.4 E09-E06 S23-S20
     10  534  634 (S36E26 L=251)   806    184        0 ( 469/ 597- 607/ 697) -14.0 E32-E19 S38-S33
     11 1178  635 (S22W18 L=295)   237    118        0 (1114/ 622-1228/ 650)  15.9 W14-W21 S22-S20
     12 1502  451 (S25W49 L=326)  1078    162        5 (1436/ 402-1566/ 500)  53.1 W45-W53 S29-S21
  

    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