Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1485 1732 (N31W59 L=156)   462     41        0 (1460/1713-1511/1747)   0.0 W56-W63 N29-N32
      2 1687 1359 (N05W51 L=147)   895    252        0 (1625/1238-1724/1445)  80.5 W47-W53 S01-N11
      3  730 1198 (N22E13 L= 84)   881    304        2 ( 599/1167- 866/1220) -23.6 E20-E04 N18-N25
      4 1662  802 (S25W34 L=131)  1683    545        1 (1549/ 594-1778/1053)  38.9 W22-W48 S34-S13
      5  978  990 (N04E03 L= 93)   231    119        0 ( 930/ 977-1032/1013) -31.0 E05-W00 N02-N06
      6  197  917 (N17E60 L= 36)  2293    297        8 ( 102/ 870- 337/ 996) -19.7 E74-E46 N12-N22
      7 1586  896 (S17W30 L=126)   428    133        0 (1547/ 840-1630/ 944)  29.7 W26-W33 S19-S15
      8  536  832 (N06E33 L= 63)   530    162        1 ( 501/ 779- 571/ 897)  50.2 E36-E30 N02-N09
      9 1832  759 (S34W52 L=148)   789    119        0 (1792/ 715-1870/ 820)  29.1 W46-W56 S35-S30
     10  393  768 (N05E45 L= 51)   350     75        0 ( 356/ 743- 434/ 794) -38.7 E47-E42 N02-N07
     11  815  710 (S07E20 L= 77)   299     93        0 ( 772/ 678- 861/ 741)   9.5 E22-E16 S08-S06
     12  393  604 (S04E51 L= 45)   810    149        2 ( 343/ 552- 443/ 669)  21.0 E57-E44 S06-S01
     13 1484  419 (S46W15 L=111)   652    116        2 (1428/ 398-1552/ 441) -32.5 W09-W21 S49-S41
     14 1166  320 (S42E13 L= 84)   388    114        1 (1121/ 287-1213/ 357)  18.4 E17-E08 S43-S40
     15 1034  229 (S45E28 L= 69)   329     46        1 (1003/ 215-1060/ 251) -11.3 E30-E25 S46-S42
     16  832  172 (S43E49 L= 48)   703     60        2 ( 796/ 151- 872/ 196) -15.9 E54-E44 S44-S41
  

    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