Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1695 (N38E08 L=242)   375     51        2 (1110/1667-1169/1719)   0.0 E10-E05 N35-N39
      2  559 1664 (N24E48 L=202)   686    106        2 ( 520/1625- 591/1702)  81.9 E49-E45 N20-N27
      3 1209 1551 (N29W00 L=250)   568    127        1 (1169/1491-1245/1612) -49.4 E03-W03 N25-N32
      4 1737 1180 (N18W43 L=293)   273     37        0 (1711/1163-1761/1195)  14.0 W41-W45 N17-N19
      5  142  907 (S27E70 L=180)   758    112        3 ( 102/ 873- 196/ 938) -10.6 E77-E61 S28-S24
      6 1548  650 (S16W42 L=292)   257    118        0 (1498/ 608-1591/ 685)  63.4 W40-W44 S19-S13
      7  288  500 (S53E68 L=182)   675     78        0 ( 256/ 477- 329/ 513)   4.4 E83-E57 S54-S51
      8 1090  386 (S44W23 L=272)   624    148        1 (1053/ 334-1143/ 439) -15.3 W18-W29 S46-S40
      9 1450  387 (S34W50 L=300)  1023    223        1 (1366/ 361-1546/ 422)  36.0 W43-W55 S37-S29
  

    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