Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1632 (N35W40 L=169)   553    157        0 (1274/1617-1417/1645) -15.3 W34-W45 N32-N37
      2  734 1453 (N39E10 L=119)   656    144        4 ( 682/1411- 780/1487) -40.6 E14-E06 N35-N42
      3  484 1234 (N30E32 L= 97)  1108    246        1 ( 415/1153- 554/1326) -52.4 E37-E27 N23-N36
      4 1002 1255 (N21W04 L=133)   930    318        3 ( 864/1200-1167/1323) -40.9 E03-W12 N14-N27
      5 1307 1037 (N02W17 L=146)  1310    262        1 (1209/ 914-1368/1154) -79.4 W13-W19 S06-N10
      6  503 1029 (N17E33 L= 96)   309    112        0 ( 476/ 993- 542/1081) -63.4 E34-E30 N14-N20
      7 1845 1036 (S13W57 L=186)   769     84        2 (1819/1001-1866/1075)  63.4 W54-W60 S14-S10
      8 1301  716 (S17W10 L=139)   796    164        1 (1228/ 648-1336/ 768) -54.5 W05-W12 S22-S14
      9  301  518 (S12E67 L= 63)  1209    245        1 ( 281/ 468- 338/ 575)  32.5 E73-E61 S15-S07
     10 1039  357 (S35E16 L=113)   899    232        3 ( 987/ 295-1072/ 422)  36.0 E22-E11 S39-S31
     11  602  331 (S28E48 L= 81)  1095    189        2 ( 538/ 296- 645/ 380) -33.7 E53-E42 S31-S25
  

    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