Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1463 1431 (N27W34 L=301)  1275    277        3 (1399/1355-1561/1535) -65.2 W29-W40 N21-N34
      2  265 1430 (N32E67 L=200)   709     70        1 ( 235/1403- 293/1460) -11.3 E76-E61 N29-N33
      3  710 1153 (N14E20 L=247)   363    109        1 ( 682/1104- 757/1181)  51.3 E21-E16 N10-N15
      4  764  889 (S03E17 L=250)   270     95        0 ( 716/ 869- 798/ 917)  21.8 E19-E14 S03-S01
      5 1469  846 (S09W28 L=294)   587    136        2 (1444/ 783-1495/ 912)  66.8 W26-W29 S12-S05
      6 1745  832 (S12W51 L=317)   943    193        2 (1697/ 737-1770/ 893)  73.3 W47-W53 S18-S08
      7 1852  773 (S18W66 L=332)   524     68        0 (1821/ 757-1875/ 792)   0.0 W61-W70 S19-S16
      8 1886  741 (S21W75 L=342)  1045     68        2 (1868/ 706-1905/ 769)  36.9 W71-W80 S22-S19
      9  418  581 (S22E48 L=219)   686    139        1 ( 398/ 522- 444/ 657)  71.6 E50-E45 S26-S17
     10  285  569 (S24E65 L=202)   932    123        2 ( 253/ 525- 318/ 618)  37.9 E72-E57 S26-S19
     11  615  563 (S24E31 L=235)   305     91        1 ( 580/ 529- 644/ 602)  45.0 E34-E28 S25-S20
     12  698  517 (S27E26 L=241)   325     94        2 ( 662/ 485- 732/ 555)  35.5 E29-E22 S29-S24
  

    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