Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1435 1545 (N36W11 L=112)  1886    380        4 (1302/1427-1541/1664) -20.0 E01-W21 N30-N41
      2  330 1577 (N10E68 L= 33)  1680    323        2 ( 273/1529- 385/1614)  -5.4 E80-E59 N06-N12
      3  710 1518 (N13E32 L= 69)   479    149        2 ( 640/1492- 794/1549)  36.9 E36-E28 N10-N16
      4  339 1473 (N03E59 L= 42)   431     73        1 ( 320/1436- 362/1512) -78.7 E60-E56 N00-N05
      5 1084 1347 (N12E05 L= 96)   715    274        0 (1004/1212-1159/1430) -32.7 E12-W02 N06-N16
      6  177 1202 (S15E67 L= 34)  1678    235        2 ( 126/1084- 210/1290)  61.7 E72-E64 S22-S09
      7 1431  541 (S21W41 L=141)   522    167        0 (1361/ 524-1500/ 565)  21.8 W34-W44 S23-S18
  

    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