Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  823 1829 (N63E29 L=  7)   516     44        0 ( 794/1809- 848/1854)  45.0 E33-E25 N60-N66
      2  905 1817 (N62E16 L= 20)   424     48        0 ( 877/1802- 932/1835) -12.5 E20-E11 N60-N64
      3 1165 1730 (N53W14 L= 50)   325     71        1 (1136/1712-1196/1746)  -9.5 W11-W17 N51-N54
      4  567 1470 (N31E35 L=  1)   404    129        1 ( 530/1425- 606/1530) -37.9 E40-E31 N28-N35
      5 1812 1415 (N26W71 L=106)   608     52        1 (1792/1390-1840/1431)   4.8 W66-W78 N24-N27
      6 1206 1397 (N27W13 L= 49)   618    134        1 (1146/1356-1257/1424) -23.2 W08-W15 N24-N29
      7  876 1308 (N21E10 L= 26)   316    175        0 ( 829/1275- 922/1350) -26.6 E13-E06 N19-N23
      8 1754  855 (S08W52 L= 88)   553     43        2 (1729/ 829-1778/ 881) -18.4 W49-W55 S09-S06
      9 1443  758 (S13W28 L= 64)   734    315        0 (1391/ 699-1524/ 839) -77.5 W23-W33 S17-S08
     10  852  686 (S18E11 L= 25)   750    173        2 ( 777/ 646- 931/ 717)  20.0 E16-E05 S20-S15
     11  395  521 (S31E52 L=344)  2633    362        7 ( 315/ 429- 493/ 638) -54.5 E57-E44 S37-S22
     12 1426  578 (S25W29 L= 65)  1142    336        4 (1311/ 540-1548/ 638)   6.3 W20-W38 S27-S21
  

    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