Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1198 1776 (N48W19 L=282)  1312    169        4 (1135/1726-1280/1840) -51.3 W14-W25 N43-N55
      2 1554 1665 (N38W48 L=310)  1714    308        1 (1394/1635-1680/1702)  -3.6 W33-W67 N35-N41
      3 1458 1528 (N28W33 L=295)  4614    723        2 (1133/1382-1740/1620) -19.2 W09-W54 N18-N35
      4 1506 1353 (N16W33 L=295)   491    159        0 (1485/1292-1540/1409)  81.9 W31-W35 N12-N19
      5  779 1352 (N18E15 L=248)   648    150        2 ( 713/1305- 842/1405) -33.7 E19-E10 N15-N21
      6 1847 1266 (N11W63 L=325)  1295    143        1 (1802/1202-1894/1310) -33.7 W57-W67 N07-N13
      7 1720 1251 (N10W48 L=311)   327     66        1 (1707/1214-1734/1291)  78.7 W46-W49 N07-N12
      8 1631 1257 (N10W41 L=303)   241     44        1 (1605/1236-1656/1282)  21.8 W38-W43 N08-N11
      9 1818 1213 (N08W58 L=320)   760     53        3 (1792/1177-1842/1247) -63.4 W54-W60 N05-N09
     10  939 1187 (N07E05 L=258)   266     74        0 ( 901/1155- 970/1219) -45.0 E06-E02 N05-N09
     11  370 1198 (N10E43 L=219)   289     67        0 ( 340/1172- 401/1218)  21.8 E45-E40 N09-N11
     12 1873  875 (S13W64 L=326)  2890    270        1 (1808/ 771-1934/1007)  54.5 W57-W72 S19-S04
     13 1576  793 (S18W36 L=299)   768    253        0 (1539/ 687-1641/ 902)  72.9 W33-W40 S24-S11
     14  385  767 (S15E45 L=218)   772    208        3 ( 295/ 742- 476/ 789)  -7.6 E52-E37 S16-S13
     15 1743  523 (S35W63 L=325)  2273    279        0 (1637/ 471-1829/ 628)  26.6 W49-W71 S39-S28
     16 1246  592 (S30W14 L=276)   258     71        0 (1217/ 561-1281/ 623)  36.9 W12-W16 S32-S28
     17 1451  465 (S39W33 L=295)   659    111        2 (1403/ 430-1494/ 507) -11.3 W27-W37 S42-S36
     18 1132  450 (S40W06 L=268)   838    179        4 (1050/ 417-1217/ 495)   9.5 W00-W12 S42-S36
  

    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