Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1463 (N13W50 L=128)  1201    190        3 (1534/1379-1662/1521) -79.7 W47-W52 N06-N17
      2  245 1460 (N52E64 L= 14)  1088     94        1 ( 196/1435- 292/1477)   0.0 E82-E54 N50-N53
      3  890 1384 (N29W02 L= 80)  1100    265        5 ( 767/1339-1023/1428) -37.6 E04-W09 N24-N34
      4 1173 1336 (N18W17 L= 95)   319    171        0 (1106/1320-1262/1359) -16.7 W13-W23 N17-N20
      5  259  902 (N16E55 L= 23)  1426    274        2 ( 232/ 794- 310/1029) -90.0 E57-E51 N09-N23
      6 1890  984 (S23W61 L=139)   539     50        3 (1867/ 956-1909/1012) -90.0 W58-W64 S24-S20
      7 1694  865 (S23W38 L=116)   535    121        0 (1664/ 817-1724/ 926)  39.8 W34-W41 S24-S19
      8 1047  568 (S21E11 L= 67)   219    112        0 ( 993/ 558-1102/ 579) -26.6 E14-E07 S22-S19
      9  530  317 (S25E60 L= 18)   916    131        2 ( 498/ 270- 553/ 386)  25.2 E70-E53 S27-S19
     10  790  290 (S33E41 L= 37)   348     54        1 ( 763/ 261- 809/ 316)  26.6 E43-E37 S34-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