Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  681 1649 (N40E24 L=213)   904    234        1 ( 562/1629- 784/1669)  -5.7 E35-E15 N37-N42
      2  188 1083 (N05E60 L=177)  1308    159        2 ( 136/1000- 231/1130) -31.6 E68-E55 N00-N08
      3 1668 1047 (S03W42 L=280)   525     84        0 (1615/1017-1710/1069)  -8.1 W38-W45 S04-S01
      4  473  786 (S14E37 L=200)   303    117        0 ( 420/ 748- 523/ 823)  24.0 E42-E33 S16-S12
      5  229  701 (S18E63 L=175)  2054    319        7 ( 190/ 622- 273/ 781)  73.3 E69-E56 S22-S12
      6 1088  748 (S20W03 L=240)   279     53        0 (1067/ 712-1119/ 778)  63.4 W01-W04 S21-S17
      7 1508  692 (S24W32 L=269)   641    205        1 (1433/ 637-1574/ 769)  35.0 W26-W36 S27-S19
      8  555  514 (S32E38 L=200)   626    183        0 ( 538/ 415- 577/ 602)  50.2 E44-E33 S38-S26
      9  370  464 (S34E60 L=178)   726     91        0 ( 348/ 415- 394/ 506)  26.6 E67-E54 S37-S31
     10 1131  315 (S51W06 L=243)   263     39        0 (1116/ 292-1147/ 339)  63.4 W04-W07 S52-S48
     11  719  282 (S51E35 L=202)   371     46        0 ( 695/ 266- 743/ 303)  15.9 E38-E31 S52-S49
     12  873  263 (S54E20 L=217)   458     81        2 ( 839/ 230- 897/ 290) -31.0 E23-E18 S57-S51
  

    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