Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1078 1700 (N40E01 L= 72)  1742    382        4 ( 918/1633-1233/1760)  -9.1 E14-W11 N35-N44
      2  714 1616 (N32E26 L= 46)   814    171        1 ( 671/1551- 753/1675) -52.1 E30-E21 N27-N36
      3 1275 1630 (N36W14 L= 87)   477     77        2 (1248/1593-1301/1662) -90.0 W12-W16 N33-N38
      4 1340 1594 (N34W19 L= 92)   460     80        1 (1315/1562-1372/1622)  38.7 W16-W21 N31-N36
      5 1510 1526 (N31W31 L=104)  1448    300        3 (1383/1465-1648/1573)  -2.7 W21-W42 N27-N33
      6 1213 1299 (N13W10 L= 82)   333    123        1 (1195/1225-1226/1362) -83.7 W08-W10 N08-N17
      7  573 1313 (N11E30 L= 42)   278     89        1 ( 534/1283- 609/1328) -18.4 E33-E27 N09-N11
      8 1799 1222 (N14W53 L=126)   490    105        0 (1766/1197-1844/1253) -23.2 W50-W58 N12-N15
      9  270 1191 (N03E53 L= 19)   305     69        0 ( 243/1161- 296/1218)  36.9 E55-E50 N01-N04
     10  257  587 (S34E59 L= 13)   641     89        1 ( 238/ 560- 283/ 618) -76.0 E61-E56 S35-S31
     11 1597  494 (S33W48 L=120)   337     61        0 (1568/ 480-1629/ 510)  11.3 W45-W50 S33-S31
  

    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