Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1428 1564 (N24W40 L=358)   260     80        0 (1402/1532-1461/1585)  15.9 W36-W43 N22-N25
      2 1109 1496 (N28W17 L=334)   372     97        0 (1062/1467-1154/1521) -51.3 W14-W18 N24-N30
      3 1828 1370 (N04W65 L= 23)  1394    178        0 (1762/1316-1902/1407)  -2.4 W55-W79 N01-N04
      4 1479 1357 (N11W35 L=352)   476    163        1 (1396/1335-1543/1393) -33.7 W30-W39 N09-N15
      5 1931 1218 (S08W73 L= 31)  1482    246        2 (1907/1124-1955/1355)  69.0 W70-W78 S13-N00
      6 1811 1264 (S02W58 L= 16)   404    102        0 (1777/1224-1834/1303)  26.6 W53-W61 S03-N00
      7  610 1158 (N19E22 L=296)   230     81        0 ( 575/1136- 643/1187)  11.3 E24-E19 N17-N19
      8 1699 1019 (S13W42 L=359)   513    107        1 (1658/ 978-1735/1066)  15.9 W38-W45 S14-S10
      9 1505  759 (S23W24 L=341)  3107    836        4 (1148/ 577-1781/ 938)   7.0 E02-W47 S27-S18
     10  478  891 (N06E36 L=282)   348     80        1 ( 438/ 858- 508/ 914) -36.9 E38-E33 N03-N07
     11 1231  856 (S11W08 L=326)   301    157        0 (1159/ 832-1303/ 890)   0.0 W03-W12 S11-S09
     12  307  837 (N06E50 L=267)   795    118        2 ( 256/ 794- 367/ 870) -32.0 E54-E45 N03-N08
     13  229  791 (N05E59 L=258)   625     78        1 ( 199/ 758- 267/ 823)  12.5 E63-E54 N03-N06
     14  918  686 (S14E14 L=304)   328    112        1 ( 903/ 640- 933/ 746)  74.1 E15-E11 S17-S10
     15 1086  484 (S30E09 L=308)   605    177        1 (1043/ 407-1141/ 542)  35.0 E14-E04 S34-S26
     16  700  370 (S29E39 L=278)   344     59        1 ( 678/ 346- 730/ 394)   0.0 E42-E36 S30-S27
     17  560  299 (S31E56 L=261)   559     42        0 ( 530/ 277- 594/ 318)   0.0 E60-E51 S32-S29
     18  844  244 (S42E38 L=280)   374     66        0 ( 811/ 220- 875/ 272)  12.5 E42-E33 S43-S39
  

    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