Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1390 1677 (N35W49 L=321)  2327    407        4 (1237/1609-1579/1722) -35.5 W37-W59 N24-N40
      2  902 1373 (N30W01 L=273)   508    108        0 ( 852/1344- 957/1391) -15.9 E02-W05 N28-N31
      3 1609 1340 (N10W45 L=317)   841    182        2 (1539/1305-1666/1380) -21.8 W41-W51 N06-N13
      4 1226 1328 (N19W20 L=292)   196     79        1 (1193/1307-1268/1341) -11.3 W17-W22 N18-N20
      5 1783 1291 (N01W59 L=331)   582     62        1 (1758/1261-1809/1323) -90.0 W56-W61 N00-N03
      6 1827  957 (S19W54 L=326)   366     42        1 (1806/ 931-1842/ 980)  56.3 W51-W55 S20-S17
      7  218  677 (N02E69 L=203)  1011    137        2 ( 200/ 611- 248/ 732)  28.6 E74-E62 S02-N04
      8  532  690 (S02E38 L=234)   557    107        1 ( 469/ 661- 569/ 709) -33.7 E42-E35 S04-N00
      9  428  561 (S08E51 L=221)   501    128        1 ( 383/ 527- 478/ 603) -56.3 E53-E47 S10-S04
     10  744  476 (S19E32 L=241)   271    111        0 ( 706/ 438- 789/ 504)  14.0 E35-E27 S20-S17
     11  813  414 (S25E30 L=242)   323     63        0 ( 797/ 380- 831/ 446)  68.2 E31-E27 S27-S22
     12  875  264 (S38E35 L=237)  1321    195        2 ( 810/ 219- 952/ 337) -81.3 E37-E32 S43-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