Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1507 1758 (N45W66 L=218)   631     60        1 (1471/1738-1535/1778)   3.6 W58-W74 N44-N46
      2 1293 1733 (N48W38 L=190)   325     40        1 (1275/1704-1313/1756)  38.7 W35-W40 N46-N50
      3  660 1699 (N56E27 L=125)   479     68        0 ( 634/1667- 690/1725)  68.2 E29-E23 N53-N58
      4  424 1672 (N55E62 L= 90)  1049     96        1 ( 371/1648- 476/1691)  -1.7 E84-E51 N53-N56
      5  774 1384 (N31E13 L=139)   225     68        0 ( 749/1356- 793/1417)  45.0 E14-E11 N28-N32
      6  606 1357 (N30E26 L=126)  1083    212        4 ( 506/1326- 711/1377)  -3.4 E34-E17 N28-N32
      7  518 1113 (N16E32 L=120)   565    201        0 ( 479/1035- 553/1191)  66.0 E35-E29 N11-N20
      8  954  986 (N04E05 L=147)   978    402        0 ( 869/ 814-1001/1129) -75.3 E08-E02 S06-N13
      9 1658  951 (S07W41 L=193)   763    146        4 (1612/ 901-1733/1004) -41.2 W37-W47 S10-S03
     10  882  768 (S08E11 L=141)   272    114        0 ( 831/ 733- 937/ 785) -26.6 E14-E08 S10-S06
     11 1832  729 (S24W60 L=212)   906     92        2 (1813/ 686-1852/ 778) -60.3 W57-W63 S27-S20
     12 1763  662 (S28W53 L=205)   672     96        2 (1736/ 621-1787/ 714)  71.6 W50-W55 S30-S24
     13 1304  623 (S22W14 L=166)   929    175        4 (1248/ 558-1359/ 694)  48.8 W09-W17 S26-S18
     14  739  625 (S16E23 L=129)   312     52        1 ( 717/ 601- 762/ 658)  45.0 E24-E21 S17-S13
     15 1456  485 (S35W25 L=177)  3021    497        4 (1299/ 403-1607/ 626)  14.6 W12-W36 S40-S26
     16  839  471 (S27E19 L=132)   436    109        1 ( 806/ 450- 891/ 501) -31.0 E21-E15 S29-S24
     17  636  400 (S31E38 L=114)   661    101        1 ( 596/ 370- 675/ 430)  18.4 E41-E32 S33-S28
     18 1466  286 (S54W34 L=186)  1108    117        2 (1400/ 266-1521/ 310) -14.7 W24-W43 S57-S50
  

    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