Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1248 1812 (N55E14 L=244)   631     88        1 (1204/1797-1299/1827)  26.6 E18-E09 N51-N57
      2  869 1795 (N38E40 L=219)   840    134        1 ( 811/1770- 942/1819)  10.3 E44-E32 N36-N40
      3  293 1523 (N07E69 L=189)  1965    163        3 ( 239/1468- 357/1577)  52.1 E77-E64 N02-N11
      4 1322 1349 (N22W09 L=267)   540    159        1 (1282/1312-1379/1406) -23.2 W05-W12 N20-N25
      5  219 1294 (S08E65 L=193)   696     81        1 ( 178/1273- 253/1320)  18.4 E70-E60 S10-S06
      6 1855  889 (N12W61 L=320)   984    146        5 (1812/ 825-1883/ 945)  63.4 W58-W65 N07-N15
      7 1511  591 (S15W43 L=302)   287     73        0 (1479/ 562-1534/ 624)  90.0 W42-W44 S16-S12
      8 1338  566 (S21W33 L=292)   267     87        0 (1325/ 524-1355/ 618) -51.3 W30-W35 S23-S18
      9  744  554 (S40E04 L=254)   407    119        0 ( 717/ 511- 796/ 601) -29.7 E07-W01 S41-S37
     10 1394  509 (S23W40 L=298)   315     89        0 (1364/ 493-1433/ 527)  36.9 W36-W41 S23-S20
  

    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