Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1331 1580 (N35W03 L=114)  2062    403        3 (1189/1399-1425/1716) -33.1 E11-W13 N24-N40
      2  560 1661 (N19E49 L= 62)   424    116        1 ( 513/1636- 596/1707)  63.4 E52-E47 N16-N22
      3  586 1565 (N13E42 L= 69)   564    140        0 ( 530/1535- 669/1607)  49.4 E46-E38 N10-N17
      4  258 1491 (N03E70 L= 41)   557     78        1 ( 240/1463- 286/1515)  45.0 E72-E67 N01-N04
      5  909 1477 (N16E19 L= 92)   206     66        1 ( 869/1464- 939/1487)  26.6 E21-E16 N14-N16
      6  988 1302 (N07E09 L=102)   241     69        0 ( 965/1276-1018/1341) -26.6 E11-E07 N06-N09
      7  251 1141 (S18E55 L= 56)   470     89        2 ( 231/1094- 272/1192)  90.0 E56-E52 S20-S14
      8 1316  605 (S22W30 L=141)   618    142        1 (1238/ 589-1389/ 627)  24.0 W24-W34 S24-S19
      9 1490  400 (S27W52 L=163)   489     43        1 (1463/ 382-1520/ 421)  21.8 W49-W54 S28-S25
     10 1650  371 (S22W68 L=179)  1183    176        2 (1594/ 334-1698/ 414)  49.4 W64-W71 S25-S18
  

    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