Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  575 1526 (N34E29 L=130)   320     77        0 ( 558/1491- 595/1568) -63.4 E30-E26 N31-N37
      2  878 1510 (N30E06 L=152)   215     72        1 ( 846/1487- 912/1537) -31.0 E08-E03 N28-N31
      3  354 1379 (N26E46 L=112)   544    133        1 ( 333/1329- 376/1462) -56.3 E50-E43 N22-N31
      4  643 1391 (N24E22 L=136)   519    122        1 ( 598/1354- 705/1430) -26.6 E25-E17 N21-N26
      5  236 1291 (N21E57 L=101)  1836    138        6 ( 176/1238- 281/1342) -26.6 E67-E51 N17-N25
      6 1376 1074 (S01W22 L=180)   452     54        1 (1345/1044-1402/1116)  71.6 W19-W23 S02-N01
      7 1865  971 (S10W61 L=220)   803    124        2 (1803/ 923-1903/1015)  -8.7 W54-W67 S13-S07
      8 1786  859 (S16W53 L=212)   507     89        2 (1755/ 826-1815/ 897)  38.7 W50-W56 S18-S14
      9 1849  858 (S17W61 L=219)   481     67        1 (1831/ 822-1865/ 894)  78.7 W58-W62 S19-S14
     10 1320  791 (S18W17 L=175)   330    113        1 (1277/ 751-1370/ 828)  45.0 W14-W20 S19-S15
     11  684  720 (S17E24 L=134)   247     59        0 ( 667/ 693- 718/ 744)  26.6 E25-E21 S18-S15
     12  261  700 (S15E59 L= 99)   431     62        0 ( 238/ 674- 290/ 723)   7.1 E63-E55 S16-S13
     13 1111  499 (S35W02 L=160)  2743    881        2 ( 847/ 366-1468/ 723)  14.0 E17-W27 S44-S22
     14 1719  656 (S29W50 L=209)   327     34        0 (1693/ 636-1742/ 676) -18.4 W47-W53 S30-S27
     15 1743  617 (S32W55 L=213)   346     61        0 (1707/ 605-1773/ 638)   8.1 W50-W58 S32-S30
     16  349  554 (S25E57 L=102)   592     86        0 ( 320/ 515- 388/ 596) -63.4 E58-E54 S27-S21
     17  734  548 (S29E24 L=134)   359     79        0 ( 700/ 524- 780/ 567)   8.1 E27-E20 S29-S27
     18 1382  249 (S60W32 L=191)  1111    126        2 (1332/ 227-1445/ 286) -14.0 W24-W44 S62-S55
  

    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