Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1767 (N48W07 L=282)  1012    140        3 (1020/1729-1129/1825) -41.2 W03-W12 N44-N53
      2 1500 1659 (N37W42 L=317)  1827    436        3 (1279/1610-1670/1692)  -7.7 W23-W61 N33-N40
      3 1358 1487 (N25W24 L=299)  4321    984        3 ( 961/1254-1746/1606) -22.5 E02-W51 N09-N34
      4  853 1488 (N27E10 L=265)   222     87        1 ( 823/1470- 898/1515) -31.0 E11-E06 N25-N29
      5  983 1395 (N20E01 L=273)   280     74        0 ( 955/1361-1017/1425) -63.4 E03-W00 N18-N22
      6 1321 1359 (N17W20 L=294)   259     86        0 (1309/1317-1336/1404) -90.0 W19-W21 N14-N19
      7  586 1353 (N19E28 L=247)   373     91        0 ( 555/1335- 630/1389) -26.6 E30-E24 N18-N21
      8 1792 1314 (N13W57 L=331)  1406    210        9 (1729/1241-1853/1384) -38.7 W52-W62 N08-N17
      9 1559 1250 (N10W35 L=310)   325     75        1 (1546/1212-1578/1295)  78.7 W34-W36 N07-N12
     10  740 1160 (N07E17 L=258)   392    121        0 ( 689/1119- 784/1208) -39.8 E20-E14 N04-N09
     11 1754  895 (S12W50 L=324)  2167    305        5 (1679/ 761-1828/1041)  64.8 W44-W57 S20-S03
     12 1414  739 (S21W24 L=299)   984    371        3 (1365/ 619-1475/ 877)  86.4 W21-W28 S28-S12
     13 1877  785 (S18W67 L=342)   846     55        3 (1847/ 749-1900/ 818) -33.7 W62-W72 S20-S16
     14  298  688 (S19E55 L=219)   532    111        3 ( 267/ 630- 340/ 727)  71.6 E58-E50 S22-S16
     15 1089  591 (S29W03 L=277)   238    100        0 (1063/ 552-1121/ 632)  63.4 W00-W04 S32-S26
     16 1665  541 (S34W51 L=325)  1701    248        2 (1556/ 480-1750/ 626)  31.0 W41-W57 S38-S28
     17 1303  460 (S39W19 L=294)   449     96        1 (1270/ 433-1352/ 500) -33.7 W16-W23 S41-S36
     18  960  425 (S41E08 L=267)   476     94        1 ( 910/ 405- 995/ 457)  15.9 E11-E04 S42-S38
  

    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