Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  665 1752 (N62E04 L= 66)   521     89        0 ( 640/1708- 692/1777) -90.0 E07-W00 N57-N65
      2 1486 1702 (N29W56 L=126)  1561    242        1 (1391/1616-1550/1777) -85.2 W51-W61 N23-N35
      3 1585 1605 (N21W56 L=126)   364     36        0 (1563/1589-1614/1621)  -9.5 W53-W59 N19-N22
      4 1293 1508 (N24W30 L=100)   308     78        0 (1249/1491-1335/1522) -31.0 W27-W32 N22-N25
      5 1077 1417 (N25W14 L= 84)  1126    230        1 ( 965/1377-1175/1439) -26.6 W07-W20 N21-N29
      6 1383 1190 (N05W24 L= 94)   399    167        1 (1321/1165-1438/1239)  -7.1 W20-W28 N03-N07
      7 1828 1155 (S10W56 L=126)   377     55        2 (1810/1123-1848/1184)  45.0 W54-W58 S11-S08
      8 1746 1086 (S11W46 L=116)   356     64        2 (1716/1057-1778/1116) -51.3 W44-W48 S13-S08
      9 1870 1071 (S16W59 L=129)   484     86        1 (1845/1032-1899/1101)  14.0 W55-W63 S17-S14
     10  666 1035 (N14E20 L= 50)   224     81        0 ( 626/1020- 700/1057) -11.3 E22-E17 N12-N14
     11  924  974 (N04E07 L= 63)   270     55        0 ( 893/ 953- 945/ 999)  14.0 E08-E04 N03-N05
     12 1854  805 (S31W54 L=124)   591    115        0 (1826/ 760-1885/ 861)  35.5 W50-W58 S33-S28
     13  709  729 (S04E25 L= 45)  1094    337        1 ( 559/ 643- 807/ 836)  13.4 E37-E16 S06-N00
     14  228  715 (N05E63 L=  7)   665    105        0 ( 186/ 679- 282/ 750) -32.0 E67-E57 N01-N07
     15  358  505 (S10E60 L= 10)   540     62        2 ( 341/ 476- 380/ 545)  33.7 E62-E56 S11-S07
     16 1509  489 (S41W18 L= 88)   631    122        1 (1459/ 444-1550/ 534)  26.6 W13-W21 S43-S38
     17  418  397 (S18E62 L=  8)   610     51        1 ( 393/ 374- 438/ 422)  56.3 E65-E58 S19-S16
     18 1236  333 (S44E07 L= 63)   401    113        0 (1186/ 314-1294/ 347) -26.6 E10-E02 S45-S41
     19 1017  249 (S43E27 L= 43)   698    120        1 ( 955/ 230-1079/ 277) -10.3 E32-E21 S44-S40
  

    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