Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1758 (N56W12 L= 11)   355     36        0 (1259/1743-1314/1771)  15.9 W08-W15 N54-N58
      2  656 1339 (N15E28 L=331)   601    262        1 ( 594/1224- 691/1432) -47.5 E35-E23 N08-N20
      3  340 1339 (N11E53 L=306)   924    142        2 ( 325/1268- 363/1412) -60.9 E56-E49 N06-N15
      4 1765 1334 (N29W55 L= 54)   511     69        0 (1731/1303-1799/1366) -26.6 W50-W58 N27-N31
      5  906 1091 (N03E08 L=351)  1253    425        0 ( 836/ 911- 991/1272)  82.6 E10-E05 S08-N15
      6 1643 1082 (N12W40 L= 39)   348     83        1 (1598/1063-1672/1098)   0.0 W36-W43 N10-N12
      7 1578  963 (N04W36 L= 35)   240     53        0 (1549/ 944-1602/ 987)  36.9 W34-W38 N02-N05
      8 1174  635 (S21W15 L= 14)  2265    614        1 ( 929/ 524-1424/ 740)  34.6 W00-W29 S31-S11
      9  411  614 (S34E42 L=317)   346     91        0 ( 397/ 575- 433/ 665) -71.6 E43-E40 S36-S30
     10  620  593 (S32E23 L=336)   389     79        1 ( 571/ 579- 664/ 608)   8.1 E26-E19 S33-S30
  

    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