Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1706 (N46E50 L=139)  1853    376        6 ( 483/1624- 570/1774) -38.7 E59-E42 N39-N51
      2 1007 1659 (N45E05 L=184)   284     47        1 ( 980/1636-1035/1676)  11.3 E07-E02 N43-N46
      3  748 1642 (N42E27 L=162)   275     44        1 ( 723/1620- 764/1665)  90.0 E28-E24 N40-N44
      4 1115 1565 (N38W05 L=193)   289     82        0 (1089/1536-1142/1595) -45.0 W02-W06 N36-N40
      5  321 1392 (N22E55 L=134)  1073    205        2 ( 304/1327- 351/1477) -59.0 E58-E50 N17-N27
      6  248 1242 (N12E59 L=130)  2993    549        8 ( 180/1153- 310/1328) -39.3 E68-E52 N06-N16
      7  396 1252 (N14E45 L=144)   634    131        1 ( 361/1194- 441/1314)  69.4 E47-E41 N09-N17
      8  456 1248 (N14E40 L=149)   293    115        1 ( 429/1195- 492/1283)  71.6 E41-E37 N10-N16
      9  255  851 (S12E55 L=133)  1980    280        5 ( 219/ 783- 308/ 929)  59.0 E59-E50 S16-S06
     10 1773  831 (S08W55 L=244)   520     91        2 (1748/ 789-1792/ 863) -59.0 W51-W56 S10-S05
     11 1102  450 (S35W08 L=197)   768    141        2 (1027/ 425-1163/ 478) -10.3 W02-W13 S36-S32
     12  732  417 (S39E21 L=168)  1094    209        3 ( 636/ 386- 840/ 448)  12.5 E29-E11 S42-S36
  

    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