Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  966 1632 (N35E19 L=179)  2548    604        8 ( 731/1517-1192/1761)  -4.2 E42-W00 N30-N38
      2  384 1596 (N20E62 L=137)   953    334        1 ( 364/1501- 404/1681) -47.5 E68-E56 N13-N25
      3 1267 1503 (N33W05 L=204)   238     70        0 (1243/1479-1304/1520)  53.1 W03-W07 N30-N34
      4  593 1453 (N15E37 L=161)   314     74        0 ( 570/1423- 620/1490)  78.7 E37-E35 N12-N17
      5  359 1280 (N00E49 L=149)   402     83        1 ( 330/1257- 391/1313)  26.6 E51-E47 S01-N02
      6 1863  962 (N12W62 L=260)  1492    249        1 (1805/ 858-1896/1038) -36.9 W53-W67 N06-N15
      7  735  971 (S11E16 L=182)   224     59        0 ( 707/ 943- 763/ 999) -21.8 E18-E13 S11-S09
      8  890  919 (S11E06 L=193)   271     51        1 ( 856/ 891- 926/ 938)  36.9 E07-E03 S12-S08
      9 1652  687 (S08W49 L=248)   570    153        0 (1637/ 634-1678/ 752) -45.0 W46-W53 S11-S04
     10 1308  633 (S19W26 L=225)   256     43        1 (1284/ 616-1337/ 647)  33.7 W24-W28 S20-S17
     11  982  581 (S29W07 L=206)   320     96        2 ( 926/ 561-1032/ 598)   8.1 W03-W10 S30-S28
  

    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