Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1802 (N37E40 L=102)  2411    387        7 ( 786/1708- 909/1891) -33.0 E53-E32 N29-N42
      2  426 1691 (N18E65 L= 77)   725    112        4 ( 397/1668- 456/1713)  53.1 E68-E62 N16-N20
      3  466 1650 (N16E58 L= 84)   478    110        2 ( 439/1621- 504/1677)  59.0 E60-E55 N13-N18
      4  446 1570 (N11E55 L= 88)  1684    290        4 ( 348/1495- 528/1646)  59.7 E60-E48 N04-N16
      5 1152 1422 (N20E04 L=138)   876    184        3 (1099/1361-1221/1470) -18.4 E08-W01 N16-N22
      6 1441  632 (S15W37 L=179)  1075    232        0 (1320/ 581-1545/ 691)  50.2 W31-W41 S21-S08
      7 1431  558 (S19W39 L=181)   379     78        2 (1395/ 530-1473/ 583)  59.0 W37-W41 S21-S16
      8 1730  514 (S11W67 L=209)   718     59        0 (1703/ 487-1756/ 543)  63.4 W63-W69 S12-S08
      9 1610  466 (S18W58 L=200)   368     43        1 (1589/ 446-1631/ 489) -63.4 W55-W59 S18-S16
     10 1569  416 (S22W57 L=200)   670    117        1 (1538/ 386-1592/ 450) -14.0 W52-W60 S22-S20
     11 1049  360 (S43W27 L=169)  1340    374        1 ( 891/ 321-1234/ 402)  19.8 W15-W40 S48-S37
  

    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