Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1682 (N34E28 L=266)  1035    163        1 ( 847/1623- 973/1723) -13.0 E33-E20 N31-N36
      2  561 1653 (N24E51 L=242)   448     57        0 ( 543/1618- 580/1685) -45.0 E53-E49 N21-N25
      3  690 1614 (N24E39 L=255)   679    211        2 ( 615/1532- 745/1669) -20.6 E47-E31 N19-N26
      4 1500 1585 (N46W20 L=314)   345     88        0 (1473/1559-1535/1600)  18.4 W17-W23 N42-N47
      5  973 1387 (N17E12 L=282)   208     27        1 ( 953/1370- 997/1408)  45.0 E13-E10 N15-N17
      6  221 1203 (S10E62 L=232)  2516    198        9 ( 177/1114- 274/1296) -45.0 E69-E55 S15-S05
      7  143 1142 (S15E72 L=221)  1464    143        2 ( 116/1089- 178/1194)  37.9 E78-E67 S18-S11
      8 1886 1143 (N27W64 L=358)  1016     76        2 (1850/1099-1909/1185)  63.4 W58-W69 N23-N30
      9 1503 1096 (N13W27 L=321)   287    104        0 (1479/1046-1525/1141) -38.7 W24-W29 N10-N14
     10 1661  552 (S14W58 L=351)   735    110        2 (1632/ 503-1696/ 590) -68.2 W55-W61 S16-S11
     11 1423  487 (S24W42 L=336)   275     60        0 (1410/ 460-1436/ 529) -45.0 W39-W44 S25-S21
     12 1512  460 (S23W51 L=344)   447    106        2 (1480/ 416-1547/ 505) -81.9 W47-W52 S26-S19
     13 1056  223 (S54W37 L=331)   500     73        0 (1036/ 199-1081/ 258) -24.0 W33-W42 S55-S50
  

    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