Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1620 1504 (N38W43 L=131)  1441    253        2 (1529/1400-1694/1596) -58.4 W36-W48 N30-N44
      2  919 1513 (N32E12 L= 75)   224     42        2 ( 900/1492- 942/1536) -33.7 E13-E10 N30-N33
      3 1069 1411 (N26E00 L= 87)   280     51        0 (1045/1394-1101/1432) -45.0 E01-W01 N25-N27
      4  810  819 (S13E12 L= 76)   343    133        0 ( 791/ 760- 828/ 892) -69.4 E13-E10 S16-S08
      5 1617  716 (S13W43 L=131)   306     74        1 (1597/ 684-1647/ 737)  56.3 W41-W45 S14-S11
      6 1163  672 (S19W12 L= 99)   251     89        0 (1142/ 643-1186/ 717) -78.7 W10-W13 S20-S15
      7 1702  684 (S14W52 L=139)   466    116        3 (1670/ 645-1737/ 708)  45.0 W48-W55 S17-S12
      8 1737  613 (S19W58 L=145)  1110    172        1 (1689/ 554-1777/ 689)  68.2 W54-W61 S23-S13
      9  425  640 (S28E41 L= 47)   318     98        0 ( 392/ 620- 466/ 673)  -9.5 E42-E36 S29-S26
     10 1247  475 (S32W21 L=109)   241     45        2 (1225/ 454-1274/ 496)   0.0 W19-W23 S33-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