Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1062 1509 (N32E03 L=287)   343     71        2 (1045/1470-1082/1551)  90.0 E04-E01 N29-N35
      2 1315 1197 (N15W17 L=307)   296     96        1 (1293/1159-1337/1247) -51.3 W14-W18 N12-N17
      3 1665 1117 (N13W43 L=333)   686    236        1 (1610/1007-1716/1206) -59.7 W38-W47 N06-N18
      4 1870 1063 (N12W65 L=355)  1659    222        1 (1850/ 959-1884/1171)  72.9 W62-W67 N05-N18
      5  139 1011 (S09E71 L=219)   721     79        0 ( 103/ 983- 175/1034)  -7.1 E81-E65 S10-S08
      6  815  887 (S09E12 L=279)   731    205        1 ( 763/ 808- 855/ 970)  78.7 E14-E09 S14-S04
      7  329  862 (S16E47 L=243)  1078    194        3 ( 261/ 781- 373/ 923) -71.6 E54-E43 S21-S12
      8  206  793 (S22E63 L=227)  2005    211        4 ( 165/ 715- 260/ 882) -70.0 E71-E56 S27-S16
      9  296  749 (S24E52 L=238)   856    103        4 ( 269/ 699- 316/ 805) -74.1 E54-E50 S27-S20
     10 1539  575 (S22W42 L=332)   539    170        0 (1478/ 545-1583/ 623)  20.6 W36-W45 S23-S18
     11 1261  402 (S37W27 L=317)   471    100        1 (1248/ 348-1272/ 453) -58.0 W24-W29 S41-S33
  

    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