Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  726 1837 (N36E52 L=105)  2010    328        3 ( 692/1773- 761/1908) -45.0 E60-E46 N30-N43
      2  344 1601 (N11E68 L= 89)  3402    488        5 ( 268/1538- 440/1659)  32.7 E76-E59 N06-N15
      3  929 1535 (N21E21 L=136)   622    118        1 ( 879/1488- 981/1567)  -6.3 E24-E15 N18-N22
      4  294  782 (S39E45 L=112)   364     56        0 ( 273/ 761- 316/ 808)  45.0 E47-E43 S40-S36
      5 1829  785 (N06W62 L=219)   481    151        0 (1784/ 770-1867/ 794)  16.7 W56-W66 N04-N07
      6 1236  726 (S16W21 L=178)   746    188        0 (1139/ 682-1332/ 784)  51.3 W16-W24 S20-S10
      7 1834  722 (N03W66 L=223)   487    108        0 (1799/ 710-1863/ 737)  20.6 W61-W70 N01-N04
      8 1564  584 (S13W47 L=204)   308     80        0 (1545/ 539-1595/ 612)  90.0 W46-W48 S16-S10
      9 1393  529 (S22W38 L=195)   668    298        0 (1306/ 470-1497/ 604)  -3.0 W29-W48 S22-S20
     10  861  446 (S44W10 L=167)  1016    245        2 ( 732/ 410- 982/ 488)  20.2 W00-W19 S48-S41
     11 1201  360 (S38W36 L=194)   399     80        0 (1154/ 343-1245/ 379)  51.3 W34-W38 S40-S35
     12 1471  294 (S31W59 L=216)   707    307        0 (1411/ 256-1525/ 337)   4.1 W51-W65 S32-S30
  

    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