Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  939 1473 (N34W09 L= 96)   255     63        1 ( 905/1454- 981/1485) -21.8 W06-W11 N32-N35
      2 1777 1358 (N01W60 L=148)   488    126        2 (1761/1312-1802/1419)  59.0 W57-W63 S01-N04
      3 1067 1282 (N19W10 L= 97)   301    126        1 (1039/1209-1099/1339)  40.6 W06-W13 N15-N21
      4 1685 1184 (S04W44 L=132)   821    137        2 (1649/1125-1720/1254)  82.9 W42-W46 S08-N00
      5 1790 1159 (S09W53 L=140)   418    153        0 (1771/1086-1814/1224)  49.4 W49-W57 S13-S06
      6 1881  968 (S23W57 L=144)  1103    241        2 (1835/ 886-1928/1058)  23.2 W50-W64 S26-S19
      7  450  991 (N19E36 L= 52)   274    109        0 ( 410/ 962- 506/1014) -39.8 E38-E32 N15-N20
      8  917  901 (N02E09 L= 78)   287     92        1 ( 872/ 865- 946/ 940)  21.8 E12-E06 N00-N04
      9 1404  756 (S19W14 L=102)   747    286        2 (1363/ 632-1438/ 910)  63.4 W09-W19 S26-S10
     10  459  650 (S01E45 L= 42)   308     71        3 ( 444/ 614- 479/ 684)  45.0 E46-E43 S03-N00
     11 1153  468 (S30E09 L= 79)   861    221        0 (1062/ 420-1214/ 562)  13.0 E15-E02 S32-S25
     12  447  335 (S21E64 L= 23)  1422    261        4 ( 407/ 290- 473/ 386)  35.0 E68-E58 S24-S17
  

    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