Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1630 (N29W64 L=229)  1162    220        2 (1518/1584-1698/1674) -37.6 W58-W72 N23-N33
      2 1371 1461 (N25W33 L=198)  1455    292        2 (1252/1370-1449/1555) -71.6 W28-W36 N18-N34
      3 1236 1157 (N11W15 L=180)   687    323        1 (1176/1095-1329/1205)  -5.7 W11-W21 N06-N13
      4 1356  923 (S06W18 L=183)   451    129        1 (1299/ 903-1425/ 944) -35.5 W14-W21 S08-S03
      5 1604  708 (S25W33 L=197)  2215    499        2 (1448/ 553-1766/ 887)  25.8 W19-W48 S31-S17
      6 1464  830 (S13W23 L=188)   373    104        0 (1438/ 787-1481/ 886)  71.6 W21-W24 S15-S09
      7 1724  690 (S30W44 L=209)  1478    254        1 (1648/ 590-1804/ 781)  26.6 W37-W53 S35-S25
      8  424  542 (S13E52 L=113)   369     66        1 ( 396/ 528- 462/ 558)   0.0 E55-E48 S14-S12
      9  525  466 (S20E47 L=118)   378     72        1 ( 508/ 430- 547/ 489)  26.6 E50-E44 S22-S17
  

    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