Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1018 1887 (N53E41 L=154)  1407    151        0 ( 958/1842-1066/1919)   4.4 E49-E30 N48-N56
      2  350 1567 (N11E66 L=129)  1096     84        2 ( 314/1535- 389/1611) -59.0 E71-E60 N09-N14
      3 1580 1479 (N39W26 L=221)   820    169        2 (1505/1452-1652/1498)  26.6 W19-W33 N34-N42
      4  389 1440 (N04E53 L=142)   794    129        0 ( 366/1371- 419/1482) -66.8 E56-E49 N00-N07
      5  179 1107 (S19E65 L=130)  1464    228        1 ( 134/1045- 210/1175)  56.3 E72-E61 S23-S14
      6  343  968 (S24E44 L=151)   686    153        1 ( 316/ 897- 368/1021) -60.3 E47-E42 S28-S21
      7  292  804 (S35E49 L=146)   631    134        1 ( 250/ 759- 332/ 870)  48.4 E54-E43 S39-S30
      8  776  760 (S26E08 L=187)   428    158        1 ( 713/ 729- 841/ 807)  -5.7 E12-E02 S26-S24
      9 1540  562 (S15W47 L=242)   314     62        1 (1514/ 539-1566/ 590)  63.4 W45-W48 S17-S13
     10 1331  512 (S25W35 L=230)   461     85        0 (1293/ 485-1371/ 547)  59.0 W33-W36 S26-S21
     11  639  270 (S68E06 L=189)   428     38        1 ( 623/ 250- 659/ 295) -38.7 E08-E02 S71-S65
  

    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