Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  644 1733 (N26E48 L=285)   741    118        1 ( 612/1696- 679/1775) -29.1 E52-E43 N23-N29
      2  459 1714 (N22E65 L=269)  1443    130        4 ( 405/1676- 511/1758)  26.6 E70-E59 N18-N25
      3  627 1614 (N18E42 L=292)   293     52        0 ( 615/1584- 640/1644) -56.3 E43-E39 N16-N20
      4 1024 1447 (N17E11 L=323)   299     53        0 (1008/1418-1043/1481) -45.0 E12-E09 N15-N18
      5 1710 1223 (N23W38 L= 12)  1321    314        2 (1660/1075-1767/1360) -47.1 W32-W45 N16-N31
      6 1832 1203 (N27W51 L= 25)   386     80        1 (1810/1162-1849/1237) -31.0 W48-W53 N24-N29
      7  903  968 (S13E05 L=328)   222     43        2 ( 879/ 942- 926/ 991)   0.0 E07-E03 S14-S11
      8  182  709 (S44E70 L=263)   720     83        0 ( 156/ 669- 210/ 744)  13.4 E83-E62 S46-S41
      9  967  604 (S32W08 L=342)   303     80        2 ( 935/ 585-1017/ 633)   8.1 W05-W12 S33-S31
     10  549  393 (S61E22 L=311)   403     53        0 ( 525/ 371- 569/ 426)  71.6 E26-E19 S63-S57
     11  708  328 (S61E00 L=333)   555     82        1 ( 672/ 296- 743/ 359) -11.3 E05-W05 S62-S57
  

    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