Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1354 1760 (N45W51 L=352)  1633    189        4 (1293/1688-1420/1805)   0.0 W41-W62 N41-N47
      2  713 1732 (N61E13 L=288)   493     51        0 ( 689/1701- 736/1760)  71.6 E16-E10 N57-N63
      3  939 1561 (N41W05 L=306)   333     82        2 ( 910/1522- 966/1596) -80.5 W03-W07 N38-N44
      4 1209 1487 (N31W22 L=323)  1342    337        2 (1069/1440-1360/1531) -21.3 W12-W30 N25-N36
      5  334 1414 (N41E51 L=250)   756    125        2 ( 293/1367- 395/1453) -19.4 E59-E42 N37-N44
      6 1147 1122 (N10W09 L=310)   327    113        0 (1089/1099-1198/1140) -26.6 W05-W11 N08-N11
      7 1740  972 (S10W46 L=347)   543    116        0 (1724/ 911-1758/1023)  56.3 W44-W48 S13-S07
      8  228  954 (N13E58 L=243)   796    127        1 ( 211/ 890- 253/1009)  82.9 E60-E55 N08-N16
      9  643  855 (N02E26 L=275)   214    104        0 ( 604/ 831- 679/ 899) -51.3 E27-E23 N00-N05
     10 1262  787 (S11W10 L=311)   399    200        1 (1231/ 704-1323/ 831)  40.6 W07-W14 S16-S09
     11 1415  530 (S32W17 L=318)  1547    353        3 (1282/ 468-1526/ 658)  20.2 W07-W26 S34-S24
     12  762  206 (S45E45 L=256)  2440    266        3 ( 652/ 162- 865/ 246)   6.5 E69-E32 S48-S41
  

    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