Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1177 1853 (N65E07 L= 81)  5313    421        7 ( 989/1772-1347/1932)  -3.2 E43-W17 N58-N72
      2 1132 1681 (N45E05 L= 84)  1399    255        3 (1015/1624-1216/1720)  -9.5 E15-W03 N41-N47
      3  542 1656 (N32E52 L= 37)  1356    301        3 ( 481/1546- 643/1708)  81.3 E57-E44 N24-N37
      4 1621 1570 (N46W44 L=133)   389     30        0 (1599/1551-1642/1589)  20.6 W40-W48 N44-N48
      5  344 1404 (N13E56 L= 33)  1654    357        7 ( 301/1277- 394/1542) -48.8 E64-E48 N05-N21
      6 1530 1187 (N17W31 L=119)  1045    293        1 (1439/1077-1605/1291) -37.6 W23-W37 N11-N22
      7 1484 1152 (N15W27 L=116)   237    105        0 (1435/1129-1526/1171)  26.6 W24-W30 N12-N16
      8 1855 1038 (N13W62 L=151)   848    176        0 (1831/ 958-1880/1103)  66.0 W59-W66 N08-N17
      9  380 1000 (S11E43 L= 46)   480    145        3 ( 319/ 948- 428/1040)  36.9 E47-E38 S14-S07
     10  158  866 (S23E69 L= 19)   801     80        2 ( 141/ 828- 176/ 896)  32.0 E74-E65 S25-S20
     11  269  783 (S26E54 L= 35)   391     96        0 ( 249/ 741- 291/ 827)  45.0 E57-E50 S29-S23
     12  633  720 (S25E21 L= 67)  1252    413        1 ( 486/ 688- 832/ 771)  15.6 E32-E07 S28-S19
     13 1136  656 (S21W13 L=102)   779    369        2 (1011/ 550-1273/ 760)  54.8 W07-W19 S29-S12
     14  382  667 (S32E43 L= 45)   322     87        0 ( 360/ 631- 418/ 688)  35.5 E46-E39 S35-S30
     15 1223  604 (S23W20 L=108)   299    138        0 (1183/ 548-1268/ 658) -36.9 W16-W24 S25-S19
     16  679  214 (S69E25 L= 64)   798     82        1 ( 653/ 181- 708/ 255)  74.7 E32-E20 S74-S63
  

    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