Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1023 1382 (N27W10 L=  1)   603    231        2 ( 897/1363-1114/1408) -19.7 W02-W16 N24-N30
      2 1524 1239 (N05W35 L= 26)   835    249        2 (1421/1169-1605/1316)   7.6 W26-W41 N02-N07
      3 1005 1248 (N20W05 L=356)   266    103        1 ( 949/1232-1052/1260)  -9.5 W01-W07 N18-N21
      4  421 1112 (N26E37 L=314)  1567    304        5 ( 305/1072- 548/1152) -19.3 E47-E27 N21-N30
      5  188  945 (N20E63 L=288)  1344    176        3 ( 114/ 908- 262/ 983) -17.5 E75-E54 N17-N23
      6 1508  946 (S10W26 L= 16)   228     70        1 (1475/ 920-1534/ 971)  14.0 W23-W27 S11-S09
      7 1177  789 (S10W02 L=353)   894    229        1 (1108/ 704-1236/ 886)  37.9 E03-W07 S13-S05
      8 1374  841 (S12W15 L=  6)   270     42        1 (1350/ 818-1404/ 866)   0.0 W13-W16 S13-S10
      9  364  765 (N07E49 L=302)   365    107        1 ( 331/ 720- 392/ 813) -80.5 E50-E46 N04-N10
     10  647  604 (S08E33 L=317)   461    151        1 ( 610/ 543- 685/ 676) -90.0 E34-E32 S11-S02
     11  925  561 (S17E18 L=333)   471    114        1 ( 886/ 514- 965/ 618)  29.7 E21-E14 S18-S14
     12  466  532 (S08E49 L=302)   487    112        1 ( 440/ 498- 504/ 589) -90.0 E52-E46 S11-S04
     13 1650  527 (S43W31 L= 22)   948    168        3 (1572/ 472-1711/ 579)   3.8 W22-W37 S45-S39
     14  762  494 (S17E30 L=320)   327    112        1 ( 739/ 448- 794/ 536)  29.7 E33-E26 S19-S14
     15  473  370 (S19E59 L=291)  1521    181        3 ( 448/ 288- 513/ 447)  30.1 E70-E50 S25-S14
     16  980  310 (S35E26 L=325)  1090    144        2 ( 906/ 203-1027/ 433)  69.1 E32-E22 S45-S24
     17 1311  217 (S57E09 L=342)  1574    201        4 (1218/ 184-1400/ 257) -12.5 E18-W01 S60-S52
     18  957  174 (S46E39 L=311)   549     42        0 ( 933/ 154- 983/ 193)   8.1 E43-E35 S48-S43
  

    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