Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1627 (N31E20 L=232)  1760    317        8 ( 879/1530-1072/1694) -15.5 E29-E11 N25-N33
      2  771 1578 (N22E31 L=221)   962    211        2 ( 715/1507- 846/1656)  85.6 E32-E29 N16-N29
      3 1482 1339 (N26W20 L=272)  2295    501        2 (1347/1126-1548/1502) -51.3 W08-W27 N14-N35
      4  150 1148 (S16E69 L=183)  3484    291        6 (  79/1029- 202/1277)  36.0 E89-E63 S23-S07
      5  846 1219 (N03E15 L=237)   190     89        2 ( 818/1205- 884/1249)   0.0 E16-E12 N01-N03
      6  399  987 (S20E39 L=213)   671    192        1 ( 333/ 927- 486/1030)  36.0 E44-E32 S25-S16
      7 1762  962 (N12W49 L=301)   391     79        1 (1745/ 919-1773/1010) -80.5 W47-W50 N09-N15
      8  833  856 (S18E07 L=245)   344    180        0 ( 802/ 796- 891/ 900)  71.6 E08-E04 S21-S14
      9  422  863 (S27E35 L=217)   481    120        0 ( 373/ 850- 488/ 880)  12.5 E39-E30 S29-S24
     10 1293  796 (S10W21 L=273)   346    125        0 (1234/ 765-1351/ 818)   7.1 W17-W25 S11-S08
     11  728  402 (S51E02 L=250)  1679    277        3 ( 622/ 317- 826/ 497) -16.9 E13-W10 S54-S46
     12  934  250 (S55W25 L=277)   317     22        0 ( 912/ 233- 953/ 268) -76.0 W22-W27 S57-S53
  

    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