Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  673 1601 (N36E33 L= 66)   353    112        0 ( 638/1558- 705/1638)  74.1 E34-E32 N32-N39
      2  300 1519 (N27E68 L= 31)  2569    222       10 ( 228/1447- 348/1592) -84.3 E79-E59 N21-N31
      3 1255 1419 (N28W13 L=112)   751    268        1 (1176/1316-1306/1528) -60.3 W07-W16 N21-N35
      4  989 1442 (N28E06 L= 94)  1157    217        5 ( 915/1362-1058/1507) -35.0 E11-W00 N22-N32
      5  357 1453 (N23E56 L= 44)   353     53        1 ( 339/1426- 377/1474)  71.6 E57-E53 N21-N24
      6  455 1362 (N18E43 L= 56)   275     75        1 ( 417/1337- 484/1380)   0.0 E46-E40 N16-N19
      7 1566 1302 (N23W36 L=135)   381     82        2 (1543/1259-1587/1340) -78.7 W34-W37 N19-N25
      8 1838 1152 (N15W61 L=160)   903     84        3 (1815/1105-1862/1194)  80.5 W58-W64 N11-N17
      9  174  908 (S12E66 L= 34)  1600    190        4 ( 138/ 843- 226/ 979) -90.0 E72-E58 S16-S07
     10  402  925 (S09E41 L= 58)   411    153        0 ( 342/ 874- 462/ 961)  26.6 E46-E36 S12-S06
     11 1212  674 (S19W15 L=114)  3219    746        8 ( 921/ 513-1435/ 913)  43.0 E03-W26 S31-S02
     12  740  799 (S14E16 L= 83)   250     88        0 ( 698/ 774- 779/ 813)   0.0 E19-E13 S15-S13
     13 1780  730 (S12W59 L=158)   611     85        1 (1752/ 698-1819/ 760)  33.7 W55-W62 S13-S09
     14  313  626 (S30E54 L= 45)  1640    218        2 ( 243/ 542- 380/ 685)  22.5 E74-E45 S36-S24
     15 1726  576 (S22W60 L=159)  2754    357        4 (1657/ 510-1795/ 656)  21.0 W51-W75 S26-S16
     16  612  431 (S41E29 L= 70)  1075    137        4 ( 561/ 386- 677/ 479)  31.6 E35-E22 S45-S36
  

    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