Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  632 1676 (N27E41 L=195)   738    175        0 ( 604/1603- 661/1758) -42.5 E48-E35 N21-N32
      2  330 1503 (N12E58 L=178)   372     71        1 ( 303/1486- 362/1522)   9.5 E61-E55 N10-N13
      3 1215 1375 (N17W05 L=241)   542    135        1 (1142/1349-1286/1392)  24.0 W00-W09 N14-N19
      4  293 1112 (S12E50 L=186)  1312    443        1 ( 173/ 956- 429/1258)  42.1 E62-E40 S21-S02
      5  159 1146 (S11E65 L=171)   715     83        2 ( 124/1118- 186/1177) -11.3 E71-E61 S12-S08
      6  698 1102 (S08E20 L=216)   255     86        1 ( 673/1057- 723/1138) -53.1 E22-E18 S09-S05
      7  143 1081 (S15E67 L=169)   922    121        4 ( 110/1038- 179/1127) -40.6 E72-E61 S17-S11
      8  225  746 (S35E59 L=177)  1394    158        3 ( 186/ 699- 273/ 802)  23.2 E67-E51 S37-S30
      9 1512  356 (S35W56 L=292)  1869    283        4 (1422/ 324-1618/ 417)  39.8 W47-W65 S38-S28
     10 1249  282 (S47W41 L=277)   399     60        1 (1209/ 268-1280/ 299)  31.0 W37-W43 S48-S45
     11 1430  229 (S45W65 L=301)   653     72        2 (1393/ 202-1461/ 252)  -3.6 W57-W73 S45-S43
  

    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