Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  374 1640 (N34E70 L=306)   767     61        3 ( 357/1608- 393/1675) -17.1 E78-E64 N31-N35
      2  583 1594 (N33E41 L=335)  1913    395        5 ( 449/1453- 682/1659) -18.4 E58-E28 N24-N37
      3  478 1450 (N22E44 L=332)   511    178        0 ( 429/1387- 512/1534)  84.8 E45-E42 N17-N28
      4  404 1482 (N23E52 L=324)   422     55        1 ( 387/1448- 417/1514) -63.4 E54-E49 N21-N25
      5 1042 1424 (N26E03 L= 13)   397    119        0 (1005/1376-1088/1450)  45.0 E05-W00 N22-N28
      6  901 1330 (N19E11 L=  5)   585    243        0 ( 786/1275- 991/1368) -12.1 E18-E04 N16-N20
      7  911 1069 (N03E07 L=  9)   208     50        1 ( 888/1050- 940/1089)  33.7 E08-E05 N01-N04
      8 1233 1027 (N03W13 L= 29)   599    160        1 (1166/ 991-1310/1054)  18.4 W08-W17 N00-N05
      9  761 1032 (S01E16 L=360)   224     84        0 ( 718/1014- 801/1053)   0.0 E18-E13 S01-N00
     10  144  929 (S13E70 L=306)  2154    155        5 ( 105/ 850- 171/ 996) -48.4 E77-E65 S18-S09
     11 1608  746 (S11W42 L= 58)  2471    548        3 (1436/ 610-1737/ 967)  53.5 W31-W49 S20-N03
     12  310  641 (S30E53 L=323)   896    138        3 ( 268/ 582- 365/ 683)  29.7 E60-E45 S34-S26
  

    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