Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1001 1911 (N62W52 L=199)   555     53        1 ( 986/1885-1018/1942)  10.8 W43-W64 N60-N64
      2  445 1466 (N47E32 L=116)   456     46        1 ( 421/1437- 473/1490)  56.3 E34-E28 N44-N48
      3 1048 1342 (N22W10 L=157)   228    101        0 (1000/1312-1089/1355)  -9.5 W06-W12 N20-N22
      4  376 1124 (N26E40 L=108)   855    246        1 ( 320/1051- 452/1170)   9.5 E45-E33 N22-N28
      5 1714  916 (S20W41 L=188)   700    179        1 (1675/ 853-1784/ 971)  14.0 W36-W48 S23-S18
      6 1632  687 (S32W31 L=178)   576    200        0 (1555/ 642-1714/ 748)   7.6 W23-W38 S32-S30
      7  868  577 (S17E21 L=126)   222     93        0 ( 824/ 548- 900/ 608)  18.4 E24-E18 S17-S15
      8  453  536 (S10E50 L= 97)   724    114        1 ( 430/ 482- 483/ 588)  29.1 E54-E45 S12-S06
      9 1412  490 (S38W10 L=158)   384    105        0 (1368/ 465-1469/ 522)   6.3 W06-W15 S38-S36
     10  911  314 (S35E30 L=117)   302     61        0 ( 894/ 286- 932/ 346)  26.6 E33-E27 S36-S32
     11  562  291 (S28E59 L= 88)  1564    150        2 ( 509/ 257- 622/ 340) -69.4 E63-E53 S31-S23
  

    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