Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1136 1633 (N41W26 L=338)   271     57        1 (1117/1605-1153/1657)  63.4 W24-W28 N38-N42
      2 1755 1432 (N11W63 L= 15)  1599    212        3 (1705/1366-1806/1521)  84.3 W58-W69 N07-N17
      3 1187 1491 (N31W23 L=335)   260     57        1 (1172/1460-1211/1516)  45.0 W20-W24 N29-N32
      4 1498 1211 (N07W33 L=345)   803    227        3 (1474/1096-1521/1318)  69.0 W29-W36 N00-N13
      5 1203 1261 (N17W16 L=328)   200     98        1 (1145/1241-1246/1272) -18.4 W12-W18 N14-N17
      6  628 1139 (N22E22 L=290)   217     77        0 ( 590/1120- 670/1152)   0.0 E25-E19 N20-N22
      7  138  988 (N19E68 L=245)  1286    283        2 (  95/ 935- 186/1041) -33.7 E72-E62 N15-N21
      8 1608  793 (S20W32 L=344)   602    138        1 (1561/ 730-1648/ 855)  32.0 W27-W35 S22-S17
      9  314  711 (N01E55 L=257)   480    134        2 ( 289/ 672- 344/ 748) -90.0 E57-E53 S01-N03
     10 1365  447 (S37W09 L=322)   275     52        0 (1343/ 426-1386/ 467) -76.0 W08-W11 S39-S35
  

    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